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REPORT ON SHOCK TUNNEL STUDIES OF SCRAM JET PHENOMENA 


NASA GRANT NAGW 674 - Supplement 5. 1989 


As in previous reports, this consists of a series of reports on specific 
project areas, with a brief general introduction commenting on each report. 
The introduction is structured by project areas, with the title of the 
relevant report stated under the project area heading. The reports 
themselves follow in the order of the project area headings. 

The commentary begins with a brief review of the program of work planned for 
1989. 
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Planned Program for 1989 


Hvnprsonle Combustion : Continue studies with a central Injection 

configuration, involving pitot and heat transfer surveys of the combustion 
wake, and the effects of pressure, Mach number and duct cross section. 

Experiments on hypersonic combustion with orifice injection from the walls 
were also requested. 

Continue development of a heated hydrogen experimental rig. 

Explore the effects of promoting mixing of a film cooling layer, in order to 
encourage combustion, and continue the development of skin friction gauges. 

Continue development of scaling studies experiments. 

rirap Measurement : Use a drag balance to measure drag on a cone, as a prelude 

to measurement of thrust. 

Expansion Tube Studies : Continue theoretical studies of expansion tube 

operation, seeking understanding of the operating limits imposed by flow 
disturbances, and continue studies of methods for predicting test section 
flow conditions. 
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SCRAM JET STUDIES 


I nvestigation of a Supersonic Combustion Layer R.T. Casey and R.J. Stalker 

This work Is a fundamental study of the wake formed by mixing of an hydrogen 
stream, issuing from the trailing edge of an Injection strut, with the 
surrounding highly supersonic air stream. A first attempt at measuring pitot 
profiles In the wake was made in 1988 but proved unsuccessful because 
unsatisfactory test flows were produced by the experimental apparatus. This 
problem has now been overcome, and measurements of pitot pressure and 
stagnation point heat transfer on small probes have been made In a flow at 
M = 4.2, and a stagnation enthalpy of 9.4 MJ/kg. 

Measurement of static pressure distributions on the wall of the duct 
Indicated that combustion was taking place, but this did not significantly 
affect the pitot pressure profiles. With an Injection slot width of 1.6 mm, 
measurements taken up to 154 mm downstream of Injection Indicated that the 
centreline pitot pressure changed very slowly with distance, Indicating slow 
mixing rates. No significant effects of combustion on the pitot pressure 
were observed, even when a silane - hydrogen mixture was used as fuel. 
Preliminary stagnation point heat transfer measurements are somewhat 
equivocal at this stage, indicating that this may perhaps, be an effective 
means of detecting combustion. 

Wall Injected Scramjet Experiments C.P. Brescianini 

Previous experiments had shown that, when hydrogen was injected in a film 
cooling mode along one wall of a combustion duct, very effective surface 
thermal protection was observed, but no significant downstream combustion 
occurred until equivalence ratios rose to almost 2. Numerical analysis 
indicated that this was because mixing between the hydrogen fuel and the air 
was suppressed by the presence of the wall. 

These experiments were aimed at promoting combustion by using ramp type 
mixing protrusions, mounted on the surface downstream of Injection. Although 
a number of different test conditions and mixer configurations were studied, 
no significant combustion effect was observed. 

Supersonic Combustion with Transverse, Circular, Wall Jets 

R.G. Morgan and R. Casey 

A series of experiments was performed in a constant area duct with fuel 
Injected through orifices in the wall of the duct. The centreline of the 
orifices was Inclined at 30° with respect to the stream direction, and the 
Injection orifices and duct cross sectional dimensions were chosen to 
maintain a 2: 1 aspect ratio for the part of the duct cross section fuelled by 
each orifice. The experiments were conducted at duct inlet flow Mach numbers 
of 4.2 and 5.5 approximately. 

It was found that substantial static pressure rises occurred In the duct when 
fuel was Injected, Indicating that mixing and combustion was taking place. 
The duct dimensions and operating conditions that supported combustion were 
found to be the same as those for a central injection strut. This suggests 
that, providing a configuration is chosen which promotes adequate mixing, the 
presence of combustion may not be sensitive to the injection configuration. 
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Dissociated Test Gas Effects on Scraajet Combustors 


R.C. Morgan 


Free stream freezing of oxygen constitutes a major disadvantage of reflected 
shock tunnels when operating at moderately high stagnation enthalpies. This 
is particularly important for combustion experiments, because the energy 
invested in oxygen dissociation represents a very substantial effective 
Increase in the calorific value of the fuel. Some indication of the 
seriousness of this effect, in practical terms, may be gained from 
experiments on the same model to be conducted at G.A.S.L. , using the 
"Hypulse" expansion tube, and at The University of Queensland, using T4. 

However, as a prelude to these experiments, the effect of oxygen depletion of 
the test gas is explored. The aim is to compensate for the increased 
effective calorific value associated with oxygen freezing by reducing the 
amount of oxygen in the test gas. In this investigation, one dimensional 
computations of premixed reacting flows are used to examine the effects of 
oxygen adjustment, as determined by a heat balance approach. The 
calculations suggest that simulation of the flow in a combustor is possible 
in the presence of free stream freezing by suitably adjusting the initial 
oxygen concentration in the test gas. 

Use of Silane as a Fuel Additive for Hypersonic Thrust Production 

R.G. Morgan 

The difficulties experienced in producing hypersonic combustion (reported in 
previous grant reports) are thought to be associated with the low pressures 
associated with hypersonic flow in the shock tunnel. One method of coping 
with this is to raise the tunnel operating pressures, and the effects of 
doing this are noted in later parts of this report. Another approach is to 
add silane to the hydrogen fuel, and this Is discussed in this section. 

It was found that silane allowed combustion at much lower pressures than with 
hydrogen alone, and may offer a means of extending combustion studies to the 
Mach numbers of 7 or 8 associated with very high velocity propulsion. 

The experiments were conducted under test conditions such that the nozzle 
stagnation pressure decayed by 15% during the test period, and the effect of 
pressure decay on the thrust measurements in the nozzle is discussed. A 
control volume approach is used, and it is argued that decay in test 
conditions has negligible influence on thrust increment measurements i.e. 
fuel-on minus fuel-off valves. A different approach is to use the Hypersonic 
Equivalence Principle in interpreting results in a decaying pressure state - 
this yields the same result. 

Pressure-Length Correlations in Supersonic Combustion 

P.A. Jacobs and R. J. Stalker 
This is the first test in what is intended to be a study of scaling effects 
in scramjet combustors. This is to be an investigation involving three 
models, of differing scales, and the present results were obtained in the 
largest of the three models. The results of this test show that the 
combustor duct is long enough to produce a significant combustion pressure 
rise (approximately 35%) at static pressures of approximately 30 kPa. This 
is low enough to permit much higher pressures with smaller models. Indicating 
that the range of pressures available can be expected to be sufficient to 
allow a meaningful study to be undertaken. 
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Hot Hydrogen Injection Technique for Shock Tunnels 


M. Vendt 


A safe, pulsed method of heating hydrogen for heated fuel Injection studies 
is required. A small gun tunnel Is being constructed for this purpose. It 
will be mounted In the dump tank of the shock tunnel and the hot, high 
pressure sample of hydrogen which Is created by the compression process will 
be fed directly to the fuel Injector assembly In the scramjet. 

Heat Release - Wave Interaction Phenomena In Hypersonic Flows 


N.R. Hard 


The management of the wave Interactions caused by combustion in a hypersonic 
duct Is expected to be a long term problem In high velocity scramjet 
propulsion. As a fundamental study of this effect, it Is planned to perform 
experiments with simple duct flows with dissociated nitrogen test gas The 
dissociated nitrogen is produced by the shock tunnel nozzle flow and In 
recombining, exhibits heat release characteristics which are not unlike a 
burning hydrogen/alr mixture. Observation of the waves produced by the 
recombination process, using differential Interferometry, Is to be used to 
Investigate thrust losses caused by wave mismatching. 


A Study of the Wave Drag In Hypersonic Scramjets 


Nesrin Aleman 


Previous analytical Investigations, using small perturbations theory, have 
shown that the effects of wave drag in a two dimensional combustion duct can 
be demonstrated by converting a Busemann Biplane to a Busemann Scramjet 
through the agency of heat addition. This study carries this approach Into 
the region where the wave Interactions are not linear, by simulating the flow 
numerically. 


Progress to date has Indicated that matching the duct to the wave patterns 
involves balancing a number of effects if It Is to be done accurately. 
However with the approximate matching which has been achieved so far, he 
variation of drag with Mach number for a non-linear Busemann Biplane has been 
studied, and is seen to exhibit maxima and minima which are 
qualitatively the same as In linearized theory. 

Parametric Study of Thrust Production in the Two Dimensional Scramjet 

G« A. Allen Jr • 

The process of thrust production In the expansion nozzle of a scramjet 
Involves interaction between the Mach number distribution In the combustion 
wake produced by the combustion chamber processes, and the expansion waves 
generated by the walls of the thrust nozzle. 

This process is investigated here using the method of characteristics In a 
numerical simulation. A simple two dimensional configuration Is considered, 
and the effects of varying parameters such as thrust surface angle and 
combustion wake Mach number distribution are considered. 


The Design of a Mass Spectrometer for use in Hypersonic Impulse Facilities 
“ K. Skinner 

In order to measure species concentrations in combustion flows, a mass 
spectrometer is being developed. This Is to use “time-of -flight" mass 
discrimination, and is expected to allow a number of scans of the mass 
spectrum at a point In the flow during each tunnel test. 
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Shock Tunnel Drag Measurement 


S.L. Tuttle and J.M. Simmons 


The initial work by Sanderson and Simmons is being extended to drag 
measurement on -a slender cone, with a 5° semi-vertex angle- This involves a 
cone approximately 400 mm long, and leads to much greater problems with 
signal to noise ratio than the Initial work. Some preliminary results are 
presented, indicating that accuracies achieved in the earlier work have not 
yet been achieved with this more difficult model. 

Development of a Skin Friction Gauge for use in an Impulse Facility 

Gabrielle Kelly 

This work is proceeding, mainly by meeting a series of difficulties and 
evolving methods of overcoming them. At this stage, measurements which 
appear to be skin friction signals are being obtained, and further refinement 
is expected to produce a working gauge. 

SHOCK TUNNEL STUDIES 


Hypervelocity Flow In Axlsymmetric Nozzles P.A. Jacobs and R.J. Stalker 


In developing a set of meaningful scaling experiments, it has been found 
necessary to design and manufacture shock tunnel nozzles for a range of Mach 
numbers and sizes. This has provided an opportunity to compare the quality 
of test flow produced by these nozzles. 

The Mach number range of the nozzles was from 4.0 to 10. The M * 4 nozzle 
produced a nearly uniform and parallel flow, according to design 
requirements. However, the M * 10 nozzle was found to be much more 
susceptible to boundary layer effects. In fact, when operating at enthalpies 
approaching 30 MJ/kg, with a nozzle stagnation pressure of only 20 MPa (both 
conditions which tend to produce thick nozzle boundary layers) it was found 
that steady flow could not be established during the test time. Whilst this 
situation could be rectified for this nozzle by raising the nozzle stagnation 
pressure, it suggests that it may be difficult to operate axisymmetrlc 
nozzles at significantly higher Mach numbers. 

Shock Tunnel Development (Supplementary Project) R.J. Stalker and R.G. Morgan 

Raising operating test section pressure levels can be done by raising the 
main diaphragm burst pressure and the driven gas volumetric compression ratio 
in a range of combinations. The trade off involved in these combinations is 
discussed. 

Hypersonic combustion of hydrogen and ethane is demonstrated. This is made 
possible by the increased operating pressure levels. 

Real Gas Effects in Hypervelocity Flows over an Inclined Cone R.M. Krek 


This study has three purposes. It is complementary to the study of drag 
forces, in that it reveals aspects of pressure measurement which need to be 
understood if pressure drag is to be calculated accurately. It is a study of 
the type of flow which is likely to characterize the forebody of a scramjet 
vehicle at high Mach number. It also constitutes a reasonably careful set of 
measurements, aimed at validating a 3-D c. f.d. code which, hopefully, may 
eventually be applied to scramjet flows. 
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EXPANSION TUBE STUDIES 


Investigation of Flow Characteristics in TQ Expansion Tube A. Neely 

This project is aimed at predicting the test section flows generated by an 
expansion tube. From the point of view of scramjet testing, the presence of 
dissociated oxygen is particularly important, and a ready means of estimating 
the conditions under which this occurs is needed. The work reported here is 
a record of an attempt to generate a theory which will readily predict test 
section compositions when dissociation fractions are low. It is being 
conducted in parallel with an experimental program, involving both air and 
argon (because of its theoretical simplicity) as test gas. 

Disturbances in the Driver Gas of a Shock Tube A. Pauli and R.J. Stalker 


This is a continuation of a study aimed at understanding the source of the 
disturbances on the pitot pressure which limit the range of operating 
conditions of expansion tubes. 

It was previously postulated that, over a substantial^^ange of potential 
operating conditions, these disturbances were associated^inclusions of driver 
gas in the test gas. In order to test this hypothesis, a differential laser 
interferometer was used to observe variations in the integral of test gas 
density across the tube as the flow passed a given station. The variations 
observed were much smaller than those expected, casting considerable doubt on 
the "bubbles" theory. 

Therefore a new approach was taken. This involved the recognition that the 
pitot pressure disturbances observed in expansion tube operation may have 
originated in the shock driver (possibly from disturbances generated at 
diaphragm rupture), before being transmitted into the shock tube test gas 
sample prior to its processing in the acceleration tube expansion. 

The acoustic disturbance modes which are compatible with the boundary 
conditions applying the shock tube are analysed, and attention is focused on 
longitudinal waves, and a first order lateral wave. The structure of the 
flow in the shock tube driver gas is analysed and it is argued that, under 
these conditions, the observed pitot pressure fluctuations are evidence of 
the dominance of the lateral wave. This is consistent with the failure to 
observe strong disturbances in the laser interferometer experiments. 

This work will continue, with a view to studying the way in which the lateral 
waves penetrate the driver gas-test gas interface, and traverse the 
acceleration tube expansion, into the test region. 
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INVESTIGATION OF A SUPERSONIC COMBUSTION LAYER 

by 

R. T. Casey and R. J. Stalker 


In 1969 a set of experiments were conducted in the T4 free piston shock 
tunnel which sought to investigate the development and mixing of a combustion 
region in a scramjet duct. The ambient Mach number and enthalpy were 4 ana 
9.4 MJ/kg respectively. A table of the other flow conditions is given in 
table 1. Hydrogen was injected into the flow from a 2-D slot type injector 
parallel to the free stream flow but at a slower velocity. The nominal 
equivalence ratio for the experiments was 1. A pitot rake and a heat transfer 
rake were used to measure the pitot pressure and heat transfer rate 
respectively. The rake was moved to various downstream positions in oraer to 
monitor the development of the pitot pressure and heat transfer within tne 
mixing region at various downstream positions. Vail static pressure 

measurements were also taken. A schematic diagram of the experimental set up 
is given in figure 2.1. 

For each downstream position the rake was moved to, a set of three shots 
were fired. For the first hydrogen was injected into an ambient 
airflow, (termed fuel on). For the second shot, hydrogen was injected into an 
ambient nitrogen flow (termed fuel into nitrogen) and for the third shot, no 
fuel was injected into an ambient airflow (termed fuel off). Comparison of 

the wall static pressure measurements for the fuel on and the fuei into 
nitrogen shots showed an increase in wall static pressure for the fuel on 

case thus confirming the presence of combustion. Comparison of the pitot 

pressure and heat transfer measurements between the fuel on and fuei into 
nitrogen shots served as a useful means of separating mixing effects and 

combustion enhanced mixing effects. 


Figures 2 and 3 each show the pitot pressure measurements taken for me 
fuei on, fuei off and fuel into nitrogen shots for the case of the rake being 
31 mm and 81mm downstream from the exit plane of the injector respectively. 

Several observations of these plots are as follows: 

l)The decay of centre line pitot pressure to its fiee 

stream value is very slow indicating that mixing is very weak, 

2 (Variations between the fuel on, fuel off and fuel into 

nitrogen profiles is very small for both the 31mm and the 81mm downstream 
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little effect on the mixing of 


cases. This suggests that combustion is having 
hydrogen and air streams 

3)A degree of asymmetry exists in the pitot pressure 

profile. 

Nitrogen was injected into quiescent air of ambient temperature and 
pressure of 40 kpa. The pitot rake was moved to a position 54mm downstream of 
the injector exit and pitot pressures were measured for this condition. T he 
pitot rake was then adjusted up half the pitch of the transducer spacing to 
increase the resolution of the pitot measurements. Again, nitrogen was 
injected into quiescent air. The results of the two shots is shown in figure 
4. From this figure, it can be seen, that the injector does not inject the 
fluid symmetrically about the centre line. This effect explains the asymmetry 
of the pitot pressure measurements ol figures ~ and j. 

Figure 5 gives a comparison of two pitot pressure profiles. The first 
was taken from a shot where hydrogen was injected into the Mach 4 air stream 
and the other a 205. silane 80S. hydrogen mixture was injected, it was expected 
that with the addition of siiane. the second case would produce improved 
combustion and this in turn wouid produce improved mixing. However, little 
difference can be seen between the hydrogen injection case and the 
siiane/ hydrogen injection case. 

The signals recorded from the heat transfer tiansducers were di a 
integrated. The value obtained from the fuel on shots and fuel into nitrogen 
shots were normalised against the value obtained from the fuel o.f shot. 
Since the fuel is injected at a cold temperature (approximately too ki tnen 
this normalised value of heat transfer is expected hi be less than unity wnen 
close to the fuel jet. When this normalized heal transfer value ,s h.gne. fo, 
the fuel on case than for the fuel off case, then it can be taken tnai the 
heat release of combustion accounts for this rise. Set of the ten neat 
transfer gauges used suffered constant malfunctions after each snot ana so no 
data was collected from these. Consequently heat transfer data Horn m,s set 
of experiments is sparse. Table 2 gives heat transfer data tor the case o 
the pitot rake being 1M«» downstream of the injector exit plane. The 
distances quoted at the top of the table are distances from the centre line 
of the duct. From this table , it can be seen that -ISmm away from the centre 
iin. the normalized heat transfer value for the fuel into nitrogen case is 
higher than for the fuel into air case. This is thought to he a boundary 
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layer effect. 


area. 

each 


It is intended for this year to conduct further experiments in this 
The mach number will be varied, up to Mach 5, and the resolution for 
pitot profile will be increased. 
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Table 2 Heat transfer values measured 154 mm downstream 
from the injector exit 
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Schematic diagram of experimental set up 
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WALL INJECTED SCRAM JET EXPERIMENTS 


C. P. Bresclanlnl 

A series of two-dimensional wall Injected scramjet experiments, 
performed in the T3 shock tunnel, have been reported previously by Morgan 
et. al. 1,2 Subsequent numerical simulation of some of these experiments, 
by Bresclanlnl and Morgan 3 , have indicated that combustion was being 
limited mainly by the fuel/air mixing rates. As a continuation of this 
earlier work a series of experiments were devised to be performed in the 
newer T4 shock tunnel at the University of Queensland. The aim of the 
experiments was to reproduce shock tunnel conditions similar to those 
experienced in T3. and also to see if fuel/ air mixing could be enhanced 
by the addition of "mixers" to the wall of the model. 

The experiments were based around the same wall injected scramjet 
model used in the previous T3 experiments. The scramjet was configured as 
a constant area duct, with a tangental wall injector, for all of the 
experiments. Three nominal stagnation enthalpy conditions, namely 8.8, 
6.6 and 3.8 MJ/kg, were chosen to match as closely as possible the 
stagnation enthalpy conditions used previously in T3. These three nominal 
enthalpies were taken from a table of standard T4 shock tunnel operating 
conditions. The actual experimental stagnation enthalpy was a function of 
the stagnation pressure at the chosen test time. 2 mm steel diaphragms 
were chosen for the compression tube/ shock tube Junction to keep overall 
pressures relatively low. This was desired so that less strain was placed 
upon the compression tube driver. 

The only shock tunnel 

nozzle available during the experiments was the nominally Mach 5 nozzle. 
To make a more direct comparison with the previous T3 experiments, which 
were performed at Mach 3.5, it was necessary to pass the flow through 
oblique shocks, created by 10* wedges, located upstream of the inlet to 
the scramjet model (See Fig. 1). The Mach number at the inlet to the 
model was estimated to be approximately 3.8, although the actual value 
was stagnation enthalpy dependant. 

The experimental results to date consist of approximately 28 shock 
tunnel runs in T4, using hydrogen, helium, silane and ethane, injected 
over a range of injection pressures. The experiments were supplemented by 
a further 27 shots using the 4 hump or single hump mixers shown in 
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Fig. 2. The 4 hump mixer Is only of a simple, preliminary design, with 
the Intention of enhancing the fuel/alr mixing rate. The mixer has a 
disadvantage in that the flow created by the mixer no longer remains 
two-dimensional', making numerical analysis of the results more difficult. 
Shock waves off the front of the humps will also produce pressure rises 
which may increase combustion, even If the mixing Is not improved. The 
single hump mixer was devised In an effort to distinguish between these 
two effects (i.e. pressure and mixing) as no significant mixing 
enhancement was expected with the single hump. The flow with the single 
hump mixer remains two-dimensional allowing easier analysis of data. 

Experimental data consists mainly of static pressure measurements 
along the top and lower walls of the model, and a few heat transfer 
measurements. The experimental data is currently being analyzed, however, 
an example of the preliminary results can be seen in figure 3. The figure 
shows static pressure along the lower surface of the model (i.e. the 
surface containing the wall Injector), normalized by the shock tube 
stagnation pressure, at an "absolute" time of 5 ms, or approximately 1.5 
ms after Incident shock reflection In the shock tube. The free stream 
test gas consists of air at a nominal stagnation enthalpy of 8.8 MJ/kg. 
The figure compares a run when hydrogen fuel was injected (RUN 709 , 
P = 200 kPa) with a similar run when no fuel was injected (RUN 712 , 
p =o kPa). There Is a noticeable, but small, difference between the 

O H 2 

fuel on and fuel off pressure traces for this low injection pressure 
case. Also shown on the figure is the result when the 4 hump mixer is 
added to the lower wall, downstream of the injector (RUN 759 , P 0 H2 ~ 200 
kPa). Unlike the result without the mixer, a significant pressure rise 
can be observed at the downstream end of the duct. This data will have to 
be looked at In the context of other experimental (and numerical) results 


before a definite conclusion can be drawn as to whether this is a genuine 
improved mixing effect, or whether it Is due to shock waves (or both). 
Unfortunately the single hump mixer was not tried at a similar hydrogen 
injection pressure in order to make a comparison with the above results. 

Experimental work on the wall injected scramjet configuration will 
continue into 1990. During the year it is hoped that a new mixer design 
will be tested, along with some perpendicular injection experiments. The 
T4 shock tunnel is now also equipped with a Mach 4 nozzle, which will do 
away with the need for the wedges placed in front of the inlet to the 


model. 
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Figure 1 


Scranjet Model In Shock Tunnel Test Section 



Figure 2(a) Four Hump Mixer Figure 2(b) Single Hunp Mixer 
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ABSTRACT 


A series of experiments arc reported using inclined circular wall jets in 
rectangular ac ramjet combustion chambers configured as constant area ducts. A 
Free Piston shock tunnel was used to create inlet flow Mach numbers between 4.2 
and 5.5, at enthalpies corresponding to flight speeds of between 3.5 and 5.3 km/s 
Measurements of duct static pressure levels indicated that significant mixing 
and combustion was achieved. Using hydrogen as a fuel, it was found that multi- 
point injection from two sides of the duct gave a larger, and more rapid pressure 
rise. 

In hypersonic flow it was found that combustion only occurred for the duct 
dimensions and operating conditions that support combustion with central 
Injection from a 2 dimensional strut. 
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1.0 INTRODUCTION 


To achieve efficient combustion In an air breathing combustor, mixing must take 
place on a scale varying from the molecular, to a characteristic duct dimension. 

In a two dimensional flow situation the mixing layer must propagate from the 
injection point to the duct walls, or to a plane of symmetry if multi layer 
injection is used, Figure 1. 

In this instance, the primary requirement is to allow adequate axial length in 
the combustion chamber for the mixing layer to speed across the duct. Because 
of the 2 dimensional nature of the flow, the transverse length scale is decoupled 
and, except for the edge effects, the inlet aspect ratio is not a consideration. 

However, in a flight configuration it may be necessary to depart from 2 
dimensional flows. Central, parallel injection has been shown to be effective 
in supersonic flow, (reference 1) . but there are likely to be structural problems 
with struts spanning the flow in real flight. Two dimensional injection has been 
used from a wall, with both parallel an transverse components, (reference 2), 
but the results were not very promising. In particular, as the combustion 
chamber mach number increases, the rate of spread of a two dimensional mixing 
layer is known to decrease, (reference 3). 

Mixing nay in principle be enhanced by using 3 dimensional configurations which 
allow the fuel jet to penetrate into the free stream flow. For instance, a 
circular Jet as shown in Figure 2 may have sufficient momentum normal to the flow 
direction to carry fuel right across the duct. 

However, the fact that injection is now from discrete stations, rather than as 
a continuous sheet means that the jet has an extra dimension to spread in, and 
that the mixing will be a more complicated process. 

For purposes of analysis it is convenient to treat the combustion chamber as a 
number of cells, each centered on a circular Jet. Apart from differing boundary 
conditions, i.e. walls as opposed to planes of symmetry, the cells are considered 
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to be Identical with respect to the »*Jor features of the flow. 

Intuitively, it would seen that the overall length scale of the cells, and the 
aspect ratio would be inportant parameters for this sort of injection. The 
length acale determines the total cross flow distance over which the fuel has 

to be transported, and the aspect ratio gives a measure of the relative 

le aulf Od v*' t**a tujo t/<vrs*'Wie d<mer\6>e>n$ 

propagation* Work performed at the NASA Langley Research Center, (reference 4), 
has indicated that a width to height ratio of 2:1 may be optimal for achieving 
good mixing throughout the cell. This shape is consistent with the assumption 
that the primary mechanism for large scale mixing is the formation of circular 
vortices on both sides of the fuel Jet. 

The experiments reported in this paper were designed to test the effectiveness 
of the injection concept in both supersonic and hypersonic flows. 

2.0 APPARATUS 

The facility used for the tests reported here was the Free Piston shock tunnel, 
T4 at the University of Queensland, (reference 5). The model used was a constant 
area duct, Figure 3 configured for single or two point injection from one, or 
both sides of the duct. The duct could be configured with an intake cross 
section with nominal dimensions 25 mm high by 50 mm wide, or 50 n square. 

From Figure 4 it can be seen that 3 basic configurations are available which 
preserve the 2:1 aspect ratio. 

Comparisons between configurations A and B, Figure 4 involve a change of length 
scale, and boundary conditions, as 2 walls are replaced by a plane of symmetry. 
Comparisons between configurations A and C involve only a change in boundary 
conditions, as one wall is replaced by a plane of symmetry. 

Using the same apparatus it is evidently possible to produce for comparison other 
Intake aspect ratios. Figure 5. Not all available configurations were tested. 
Injection was through sonic throats, inclined at 30* to the flow direction and 
flush with the wall surface. Figure 6. The fuel was supplied from a Ludweig tube 
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reservoir, (reference S), with equivalence ratio being controlled by the 
reservoir pressure. When compering single end Multipoint injection, i.e. 
configurstions A end B Figure 4, the total throat area was maintained constant. 
This was done ao that the priaary Jet parameters of dynamic pressure ratio, aach 
number and teaperature would be conserved at a given equivalence ratio. 

In calculating the fuel Jet flow properties, it was assuaed that the injectant 
expands isentropicaly froa the throat conditions to the free streaa static 
pressure level . 

3.0 RESULTS 

In table 1 a manifest of the runs performed in this series of tests is presented. 
The characteristic dimension chosen to categorize the cell sizes is the diagonal. 
Note that this refers not to the external duct geometry but to the size of the 
similar interval units which the injection pattern creates Figure 4. In table 
2 a summary of the test conditions is given. 

The first test was primarily to determine if combustion would occur at all, and 
configuration A, Figure 4 was used. The results are displayed for an enthalpy 
of 6 mj/kg in Figure 7, as a plot of static pressure against distance down the 
4 uc t s The static pressure levels are normalized by the shock tube stagnation 
pressure to compensate for minor variations between runs . An Indication of 
combustion was given by comparing the relative pressure rises produced by 
injecting fuel into air, and fuel into a flow of nitrogen at the same conditions. 

At this stage, a problem of interpreting the results of 3D flow data becomes 
apparent. In 2D flow a single static pressure monitoring point is adequate to 
fully define the wall static pressure at any axial location. 

To do the same thing in a flow where the pressure may vary in any of 3 dimensions 
requires a 2 dimensional array of pressure transducers for all surfaces of 
Interest. Such instrumentation was beyond the scope of the present series of 
experiments , and the pressure monitoring points were mainly located on the center 
line. This introduces the possibility that the pressure at the point chosen for 
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th . Miiuria*nt might be significantly different than the average wall static 
pressure at the axial location. This data ahould tharafora ba Intarpratad with 
care, especially when trying to compare the efficiency of different Injection 
configurations which produce different cross flow pressure profiles. 

gone Indication of the order of the transverse pressure changes which apply Is 
given by the 2 transducers mounted 335 am from the Intake. One of the- la on 
the centerline, and the other la offset. A difference In pressure was so.etl.es 
measured by these transducers, as seen for Instance In the Injection Into 
nitrogen data In Figure 7. Transverse pressure gradients In supersonic flow 
create a wave pattern, the overall effect of which Is to reduce the pressure 
difference. The data was obtained In a duct with a length to height ratio of 
30, and a length to width ratio of 15. This would allow for multiple reflection 
of waves within the duct, so that a local perturbation to the transverse pressure 
profile would not be maintained throughout the length of the duct. In Figure 
7 It Is seen that the Increased pressure measured when injecting fuel Into air, 
as opposed to nitrogen, Is sustained over a substantial length of the duct. This 
Is taken as evidence of a genuine co-bus t ion effect, rather than an accident of 
transducer location In a complicated flow field. 

Another Indication that transverse pressure gradients, in either direction, are 
not undue ly distorting the measurements is given in Figure 8. Pressure 
measurements were taken on the top and bottom surfaces simultaneously, and are 
represented by the symbols + and x, without lines. The axial development of 
pressure with distance is seen to be sensibly the same on both sides of the duct. 

Using the same configuration, the tests were repeated at two different atagnation 
enthalpies. The results, presented in Figure 9, show no significant temperature 
dependency. This is surprising, as the static temperature levels were varied 
between 1200 3075k, over which range a significant change in combustion heat 

release would be expected. The injection conditions for the 3 different 
enthalpies are included in table 1A, and are seen to be similar in Jet mach 
number and dynamic pressure ratio. 

The model was then re -configured as shown In Figure 4B, in order to investigate 
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th« effect. of varying th. longth icilt of the mixing call., but maintaining tha 
2:1 aspect ratio. In Figure 8 a comparison is shown at an enthalpy of 6I.JA* 
and an equivalence ratio of approximately 1 of the effect of the cell size. It 
is seen that a significantly increased pressure rise occurred with smaller cells. 
This is an intuitively sensible result, a. the large scale eddies have to decay 
over a smaller range before molecular .Ixing is complete. 

Fuel was again injected into nitrogen test gas to confirm that the pressure rise 
was combustion induced. In Figure 10. the fuel into nitrogen profile, for both 
configurations are shown. It is seen that without reactions, the presence for 
4 inert Jet. does not significantly after the duct static pressure levels. 

This configuration was also tested at 13.9 mj/k 8 . Figure 11 and it was again seen 
to give larger pressure rise than the single Jets. However, once again it was 
noted that the pressure profiles were not strongly temperature sensitive. (Figure 
12) . Combustion of a given fraction of the fuel would be expected to generate 
a larger pressure rise in the low enthalpy, cold test gas. The fact that this 
effect was not observed suggests that either more fuel was burned in the high 
enthalpy flow, or that for some reason the heat release was enhanced. 

The dissociated test gas at the higher enthalpies provides a possible mechanism 
for extra heat release due to the heat of formation of oxygen molecules. When 
combustion is from atomic rather than molecular oxygen the effective heat of 
combustion is increased. However, this would be in contradiction to previous 
results in a supersonic constant area duct using two dimensional injection, where 
the pressure rise was lower at the higher enthalpies. 

The model was then reconfigured a. a 50x50 - 2 duct with injection fro. single 
holes on the top and bottom surfaces, as shown in Figure 4C. The aspect ratio 
of 2:1 was preserved, the purpose of the tests was to investigate the effects 
of changing the boundary condition fro. a wall to an axis of symmetry. 
Unfortunately, the model was assembled with nozzles of a smaller area than was 
used for the 25*50 mm 2 duct test, in Table 1A. This meant that it was not 
possible to use the same jet conditions to match equivalence ratio, for the 2 
tests. It was chosen to keep the .me injection parameters and to reduce the 
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equilvalence ratio to nominally 0.5. Tha results of th# two conf lguratlona ara 
tharafora not diractly eoaparabla, but uaaful data was navarthalaas obtalnad. 


Tha runs performed in this configuration ara tabulated in Table 1C. Tha results 
arc summarized in Figure 13, plotted in tha fora of pressure difference between 
injection into air and injection into nitrogen. The saaller duct, with each 
injection cell enclosed by four walls was seen to give a lower pressure rise 
downstrean in the duct than the larger two cell configuration. However, the 
difference in the pressure rises is roughly i n proportion to the equivalence 
ratios, and it is tentatively concluded that the changing boundary condition did 
not leave a major lapact on the flow. 

Yhe »odel was then reconfigured with a single Jet as shown in Figure 5C, giving 
a square aspect ratio in a 50x50 an 2 duct. The runs for this geometry were done 
with the aach 4 nozzle as tabulated in Table ID. The equivalence ratio was set 
at 0.5, and the total nozzle area was unchanged so that a direct comparison with 
the results of Table 1C could be Bade . 

The results of these tests are shown in Figure 16, in the fora of the pressure 
differentials between injecting fuel into air and into nitrogen test gases. The 
shape of the surface pressure profiles is evidently different in the two cases , 
but the overall pressure levels were similar. On the basis of this single 
comparison, it would be premature to make any conclusions as to the relative 
efficiencies of the two aspects ratios. 


The next series of tests were performed using the hypersonic M5 contoured nozzle 
to evaluate the circular Jets in a hypersonic flow field. The model was 
configured as in Figure 4B, which had been found the most effective for the 
supersonic flow. The runs performed in this node are listed in Table IE. 

The results of these tests are shown in Figure 14. A substantial pressure rise 
was measured when injecting fuel into air, indicating that circular wall Jets 
may be an effective injection technique in a hypersonic duct. This is similar 
to the results reported in reference 7, where with central injection it was 
concluded that the compressive effects of the wall boundary layers were required 
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to initiate combustion. 


When the duct ««« subsequently enlarged to • 50x50 ..ctlon combustion could 
no longer be sustained. The dashed line In Figure 14 Indicate, the theoretical 
fuel off pre.aure level in the absence of boundary layer. . and the duct pre.eure 
Is seen to be mich higher than this . 

The circular Jet model was therefore .et up a. a 50x50 mm 2 duct to see if the 
sane effect reported in reference 7 would be observed. The configuration of 
Figure 4C was used, and the runs are listed in Table IF. 

The results are shown in Figure 15, and it is seen that no combustion was 
evident . 

4.0 CONCLUSIONS 

The results show that 13 wall injection through inclined circular Jets can be 
affective in supersonic and hypersonic combustion chambers. Most of the tests 
ware performed with mixing cells having a width to height aspect ratio of 2:1. 
Duct static pressure levels rose higher, and more rapidly when the transverse 
length scale of the mixing cells was reduced. 


Preliminary tests indicated that the mixing cell dimensions are the primary 
parameters governing mixing and combustion, when the injection conditions were 
maintained constant. Replacing wall boundary conditions with planes of symmetry 
did not appear to have a major influence on performance in supersonic flow 
although this test was not rigorous, as injection parameters were changed at the 
same time. This result, if confirmed, is significant in that combustion chamber, 
could be designed with multiple internal cells, centered on wall Jets, based on 
the results of scale testing of a single cell. A single comparison where the 
mixing cell aspect ratio was changed to 1:1 did not show a significant 

difference . 

In hypersonic flow the wall Jets were found to produce combustion where it had 
previously been observed with central injection. The range of conditions where 
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coabuatlon was auatalnabla In hypersonic flow was not axtandad by tha usa of tha 
wall Jata. 

axial praasura profiles wara found to ba largely Independent of duct static 
temperature , and this la not understood at present. 
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NOTATION 


H. 

P 

U 

H 

Me 

X 

UJ 

p 

p^ 


Stagnation Enthalpy 
Duct Static Praaaura 

Stagnation Preaaure Measured In Shock Tube 
Velocity of Free Stream 
Mach Number 

Convective Mach Number Based on Axial Velocity Component 
Distance from Intake 

Velocity of Injectant after Isentropic Expansion to Free Stream 
Pressure 

Free Stream Density 

Density of Injectant after Expansion 
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TABLE 1. RUM MANIFEST 


A. 50x25 ■»* Duct, Single 
Mach 4 Nozzle 

Injector 

Diameter 

7b, Configuration Fig. 

4A, 

Run Nunber 

H*(aj/kg) 

4 


MJ 

P^f/P u 2 

Me 

1380 

6.1 

0 

2150 

2.38 

0.29 

0.782 

1381 

6.1 

1 

2150 

2.38 

0.29 

0.782 

1382 

6 . 1(N2) 

- 

2150 

2.38 

0.29 

0.782 

1385 

6.1 

0.8 

20% SjH^/Hj alxture 


1383 

9.1 

1 

2072 

2.2 

0.32 

0.97 

1384 

13.9 

1 

1990 

2.04 

0.31 

1.25 


B. 50x25 am 2 

Duct, 4 Injectors Diameter 3.45 an, 

Configuration Fig. 

48, 

Mach 4 Nozzle 





Run Number 

Hs(mJAg) 4 


MJ 

^//pu 2 

Me 

1386 

6.1 1 

2150 

2.38 

0.29 

0.782 

1388 

6.1(N2) 

2150 

2.38 

0.29 

0.782 

1389 

13.9 

1960 

1.98 

0.301 

1.25 

1391 

13.9(N2) 

1960 

1.98 

0.301 

1.25 

1392 

13.9 0 





1383 

9.1 1 

2072 

2.2 

0.32 

0.97 

1384 

13.9 1 

1990 

2.04 

0.31 

1.25 


C. 50x50 an 2 Duct, 2 Injector* Diameter 5 an. Configuration Fig. 4C, 
Mach 4 Nozzle 


Run Number 

Hs (aj /kg) 

4 

UJ(a/.> MJ 

Pj Uj 2 /*u 2 

Me 

1399 

6.1 < 

0.5 

2210 2.53 

0.33 

0.76 

1400 

6.1(N2) 

- 

2130 2.33 

0.315 

0.79 

1401 

D. 50x50 am 2 

6.1 

Duct, 1 Injector 

0 

Diameter 7 am. Configuration Fig. 5C, 


Mach 4 Nozzle 
Run Number H* (aj /kg) 

4 

MJ 

pfif/p u 2 

Me 

1403 

6.1 

0.5 

2170 2.43 

0.32 

0.77 

1400 

6.1(N2) 

- 

2170 2.43 

0.32 

0.77 


1401 


6.1 


0 
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D. 50x50 b 1 Duct, 1 Injector Diaaeter 7 am, Configuration FI*. SC, 
Mach 4 Nozzle 


Run Number 

Ha(mJ/kg) 

4 

UJce/.> 

MJ 

PjOj a /pu 2 

Me 

1403 

6.1 

0.5 

2170 

2.43 

0.32 

0.77 

1400 

6.1(N2) 

- 

2170 

2.43 

0.32 

0.77 

1401 

6.1 

0 






E. 50x25 b * 2 Duct, 4 Injector* Diameter 3.45 bb, Configuration Fig. 4B, 
Mach 5 Nozzle 


Run Number 

Hs(mjAg) 

4 


MJ 

PjUjVpu 2 

Me 

1394 

8.4 

1 

2290 

2.74 

0.265 

1.14 

1395 

8.4(N2) 

- 

2290 

2.74 

0.265 

1.14 

1396 

8.4 

0 





F. 50x50 mm 2 Duct, 2 Injectors Diaaeter 5 mm, 
Mach 5 Nozzle 

Configuration Fig. 4C, 


Run Number 

Hs(mJ/kg) 

4 

UJ(./.> 

MJ 

PjUjVpu 2 

Me 

1397 

8.4 

0.5 

2290 

2.74 

0.265 

1.14 

1398 

8.4(N2) 

- 

2290 

2.74 

0.265 

1.14 

1396 

8.4 

0 






TABLE 2 TEST CONDITIONS 
A. Nominal Mach 4 Nozzle 






Fligiht Speed 
Km/* 


T(k> 

U («/») 

,<kg/m s ) M 

pu Kg/s/m 2 

3.5 

6.1 71 

1200 

3110 

0.2 

4.5 

622 

4.3 

9.1 97 

2025 

3520 

0.165 

4.2 

580 

5.3 

13.9 109 

3075 

4120 

0.1168 

4.35 

481 


B. Nominal Mach 5 Nozzle 

Flight Speed P(kpa) T(k) U (-/#) p(kg/m 9 ) M pu Kg/a/m 2 

Km/s 

4.1 8.4 20.5 1160 3700 .0609 5.54 225 
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wall or plane of symmetry 


freestream flow 



Fig 1, Two dimensional mixing layer. 



Fig 2 Transverse circular wall jets 

Shock tunnel nozzle 

injectors ^ 

constant area duct 

Fig 3 Constant area duct configuration. 
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a furl/W rwi 1382. M» *"y«9 connection 


aoo 300 400 9n 


000 1000 


Fiq 7 P/X 25*50nn duct, single injector 
supersonic Flow BM4, Hs 6nj/kg 


P/Pstg 


Kin MOL <n* taro* hoU on low wfocQiCTift;. 

♦ x !£r5n frwn *W writer 41 An 1* 
' Hi-taJ/ta* tn4, IS»J0" 


nil 


100 ION 


Fig 8 P/X single and double sided injection 
25*50 nn duct, supersonic flow 8 M4 


P/Pstg 


■ Run 1384, Hs lSnj/kg 
o Run 1383, Hs 9.1 nj/kg 
X Run 1381 Hs 8.1 nj/kg 



furl off thfory 


000 1000 


Fig 9,P/X 25*50 nn duct, single jet, M=4 
effect of stagnation enthalpy 




Fir 10, Hydrogen injection into nitrogen, 25«50 mn duct. 
Effect of injector configuration. 



Fig 11 Effect of Internal length scale, 25«50mm duct. 
P/Pstg 



Fig 12 Effect of enthalpy change, 25*50nn duct. 
Fuel Into air A point Injection. 
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Fig 13. Pressure differentials, injection into air 
ninus injection into nitrogen. 


P/Pstg 



Fig 14, Hypersonic flow, 25*50 nn duct, 4 point 
Injection, nltrogen/alr comparison. 


P/Pstg 



Fig 15, Hypersonic duct 50*50 mm, 2 point injection 
comparison of injection into air and nitrogen. 
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L “ AIR 



Fig 16, Effects of internal aspect ratio on pressure 
differential, large duct 50«50mm. 
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DISSOCIATED TEST CAS, EFFECTS ON SCRAMJET COMBUSTORS 


\ 


R.G. Morgan, Mechanical Engineering, The University of Queensland 


ABSTRACT 


The problem of the dissociated test gas produced by reflected shock 
tunnels is addressed as it relates to heat release in supersonic 
combustors. When hydrogen burns with dissociated oxygen atoms, the 
heat of formation of the oxygen adds to the normal enthalpy of 
combustion. It is proposed that this effect be compensated for by 
reducing the oxygen concentration of the test gas. A procedure for 
doing this is presented based on a simple heat balance approach. One 
dimensional computations of premixed reacting flow are used to support 
the calculations. It is suggested that direct simulation of the flow 
in a combustor is achieved if the oxygen adjustment is made. 


NOTATION 


a mass fraction of total oxygen in the form of 0 atoms 

•< ii ■ » " * “ " " NO molecules 

®H 

^ equivalence ratio 

*•= 4 > for * < 1. * 1 for * > 1 

h c heat of combustion of H2 from 02 (90 MJ/kg H2) 

hf heat of formation of 0 atoms from 02 (15.7 MJ/kg) 

h„ heat of formation of NO Molecules from N2 ♦ 02 (2.93 MJ/kg) 

h. heat release per of combustion products 

02 mass fraction 02 

0 mass fraction 0 

N2 mass fraction N2 

NO mass fraction NO 

X mass fraction total 0 content 

X, mass fraction total content in "ideal" air 

P static pressure 

u velocity 

p density 

T temperature 

M Mach No. 

H. stagnation enthalpy 
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INTRODUCTION 


Wind tunnel facilities have traditionally had a dual role as a tool for 
the Investigation of Aerodynamics phenomena. Firstly they may be used 
to study fundamental physical processes, using models which may not 
necessarily relate directly to flight configurations. Secondly, they 
may be used for the testing of models where the results may be used 
directly, with appropriate scaling laws, to Indicate the performance of 
flight hardware. Since the advent of high speed computers wind tunnels 
have increasingly been used for the valldif lcation of computer codes 
rather than for direct simulation. This application may be said to 
combine the two roles mentioned above, with the code being validated 
using contrived models, but subsequently being applied to real flight 
configurations. 


In recent times some doubt has been expressed as to the value and 
usefulness of the direct simulation role of high speed wind tunnels. 
This particularly relates to hypersonic combustion where the 
dissociated test gas produced by ground facilities leads to changes In 
chemical kinetics and an Increase in the effective enthalpy of 
combustion. It is the contention of this paper that these problems can 
be adequately addressed, and that direct simulation is still necessary 
given the current status of numerical and analytical understanding. 


The test flow in a reflected shock tunnel Is created by means of a 
steady expansion from the stagnation region in the shock tube. The gas 
In the stagnation region is approximately In a condition of equilibrium 
composition due to the high density levels and the low velocity. At 
the stagnation enthalpies corresponding to high flight speeds, the 
stagnation temperatures are such that significant dissociation occurs 
to the nitrogen and oxygen molecules. As the flow starts to travel 
down the expansion nozzle, the temperature drops and the equilibrium 
composition shifts towards a reduced dissociated content. However, due 
to the high particle velocities and lower density in the nozzle, 
reaction rates may in some circumstances be inadequate to maintain 
equilibrium composition. The test gas at the nozzle exit will 
therefore have a frozen content of dissociated radicals created In the 
stagnation region. The higher the stagnation enthalpy, the higher will 
be the concentration of dissociated gas. 


The presence of the dissociated components may effect the flow in 
several ways. Some of the stagnation enthalpy will be stored as the 
atomic heat of formation, and will not be converted directly to thermal 
or kinetic energy. It is important, when comparing the results of 
shock tunnel tests with real flight, or with other tunnels which do 
not produce dissociated air, to match the total enthalpy Including 
thermal, chemical and kinetic components. If the flow is subsequently 
processed by a strong shock, then the Increased pressure and density 
levels lead to rapid reaction rates and the establishment of near 
equilibrium conditions, with the recovery of the heat of formation. 


41 



When the flow. Is subsequently expanded around the body It will be 
Insensitive to' the gas condition upstream of the shock. In these 
circumstances very little difference will be observed between shock 
tunnel tests and tests on an Ideal wind tunnel producing undissociated 

test gas. 


For the purposes of aerodynamic testing. Mach number Independence 
applies for flows behind a strong shock, and matching of the total 
enthalpy has been shown to give good simulation of flight conditions^ 
However for shock tunnel testing of scramjets the situation Is slightly 
different. The flow is never processed by a strong shock, and non 
equilibrium gas composition persists throughout the duct. 


The presence of free stream oxygen atoms may have a significant effect 
on combustion and it has to be examined closely. There are three main 
mechanisms by which the free stream oxygen atoms may be seen to enhance 
combustion. and influence comparisons made between shock tunnel 
experiments and real flight. 


Combustion heat release . When combustion is from the atomic 

oxygen form, the heat of formation is added to the combustion 
heat release. This Increases the net energy release from 
approximately 90 MJ to 215 MJ per kg of hydrogen burned 
This Is clearly a substantial Increase, and relatively small 
oxygen dissociation fractions may be significant. 


Ignition delay times . Forming of oxygen radicals Is a 

necessary step In the combustion of hydrogen, and the pre 
dissociated shock tunnel flow may be expected to react faster 
than ordinary air. 


(ili) Mixing Effects . In some circumstances a scramjet may be 

considered to be a pure diffusion flame, where reaction rates 
are so fast that heat release is mixing controlled. Oxygen 
molecules have different diffusion rates than the atoms, and 
the local chemical composition will effect the rate of beat 
release. Furthermore, the macroscopic development of the 
fuel air mixing layer Is dependent on the local flow 
properties such as Reynolds number, this may also be 
influenced by the dissociated oxygen content. 


For a given static pressure, dissociated oxygen increases the partial 
pressure of the oxygen components, effectively enhancing mixing and 
combustion. 
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In this report only the first mechanism will be addressed, that of 

-s 

SKEtr. \t SK2, ‘S-TSUrTu^ «3 ,"*£ 0 £. 
in th ;? t 2 Sf^iJSiFS’srss: ssiStSL - .now 

tunnel ~P.rl.ents end tests us In, 

undissociated air. 

However, the heat release per unit mass of fuel ls stl “ h \ g . he ^ 
despite the reduced free stream oxygen content, and this should be 
taken Into consideration when comparing specific Impulse measurements. 


ANALYSIS 

The comparison of -lnd tunnel tests with flight ®“ er *} 

Involves scaling of a number of flow parameters. When dealing with 
hieh energy shock tunnels, an additional uncertainty arises due to the 
difference* in gas composition. The scop, of this paper is to address 
tie difference between a wind tunnel using dissociated oxygen and 
nitrogen and an Idealised wind tunnel using air with normal 
atmospheric composition. The more generalised problem of c °“P ari ^^ 
Idealised wind tunnel and free flight Is well understood, and is not 

considered here. 

Analysis of the reacting nozzle flow using a non-equilibrium chemical 

kinetics code. Ref. 1 Indicates there are 4 sl « n J fl ^“ t 

as far as hydrogen combustion Is concerned, namely N2, OZ, nu, u. 


Combustion to water from the 3 oxygen bearing constituents proceeds 
with heat release as Indicated below 


Cl) 

(11) 

( 111 ) 


H2 + 02 ■= H20 + 120.6 Mj/kg H2 

H2 + 0* H20 ♦ 246 Mj/kg H2 

H2 + NO=H20 + H N2 + 164.6 Mj/kg H2 


Note 

schemes, 


These equations are not Intended to represent reaction 
but are used to calculate energy release between end products. 


Combustion of undissoclated air Involves only reaction from ‘ 
molecules, and may be seen to give significantly less ^ . ^ 

the two other main combustion paths. The approach of this paper Is to 
apply a correction to the free stream oxygen concentration In order to 
match the heat release. In other words the shock tunnel tests would be 
performed with a reduced oxygen concentration Intended to reproduce the 
results of combustion In atmospheric air with the same flow properties. 

The primary measured experimental parameter In scramjet testing Is duct 
static pressure level, and corrections made to the oxygen levels would 
be calculated to produce identical pressure/distance profiles in the 
two test facilities. To match this exactly a numerical simulation is 
required which fully evaluates the changes in combustion species 
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composition produced by changing the free stream oxygen level. 
However, the approach of this paper Is to do a first order correction 
based on the heat release of the three reaction mechanisms Identified 
above. Numerlfcal simulation Is then performed to assess the validity 
of this approach. A series of experiments Is envisaged for 1990 to 
evaluate the concept. The analysis of this preliminary paper Is based 
on ID premixed fuel/alr flow. 

The heat release of 120.6 mj/kg for reaction (1), combustion from 
oxygen molecules. Is based on combustion to water alone. The 
equilibrium composition of a hydrogen flame contains dissociated 
radicals, *rtri the full heat release Is not achieved. A heat release of 
90 mj/kg H2 Is more representative of that developed In a real 
combustion chamber, and this value is used here. 

For a kilogram of test gas the number of moles of each species will be 


of atomic oxygen 
of NO 

( 1-tt-OjJ of molecular oxygen 


jl - X |l ♦ a nj|] j fg of »°lecular nitrogen J 

This will require + if 0 } ^ moles of molecular 

hydrogen to establish the equivalence ratio 4- 


as follows: 

Xa 

16 

16 
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The Idealised reaction scheme for the 4 component air mixture Is 
therefore as given below: 


101 • INO] ♦ jSj (1-a-o.) 1021* |l-x(l«o„ 55 ]} 55 [N2I 

‘{I * x ‘ 1 Ti s) * IS” } ♦ ,H2) 

■ {ll * T?" * 16 55 IN21 

t t t t t 


> 


from 0 from NO from 02 from NO original free H2 

(add h f ) (add h n ) (h c ) 

. jL (l-«-«.)(l-*'HO2).jg»(l-**HNO]*g(lVH01*55(f-*‘) IH2) 

t t t t 

residual 02 residual NO residual 0 residual N2 

+ heat release 2 h c ♦ Xa h f * Xa,, h„ ygj^* 

t t T 

heat release from heat release heat release 

converting converting 

H2, 02 -» H20 0 •* 02 NO -» N2, 02 

where [ ] denotes a mol of a species. 

The function is used to allow for equivalence ratios other than 1 
where there Is either unburned fuel of oxidant. 


For 


# < 1 . # 
# > 1 . - 1 


The coefficients on the RHS of Eq. (4) may evidently be simplified 
slightly, but they are presented here In full In order to clearly 
Illustrate the origin of the end products. For equivalence ratios of 
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less than 1 It Is assumed that the 0, 02 and HO are consumed in equal 
proportions. 


Lowering the oxygen mass fraction reduces the heat release, as may be 
seen from Inspection of Eq. 4. However, the question arises as to the 
flow parameter against which the heat release should be conserved when 
comparing different test facilities. With respect to the fuel content, 
energy release is a function only of a and and cannot be matched by 
changing the total oxygen mass fraction. Comparisons of measured 
specific impulse therefore will still have to be adjusted to correct. 
For the effective change In fuel calorific value even if oxygen 
depleted test gas Is used. 

Ideally the heat added due to combustion at a given equivalence ratio 
should release the same amount of heat when normalised by the 
stagnation enthalpy. When comparing results at a given stagnation 
enthalpy, this Is equivalent to matching heat release per unit mass of 
air flow. However, reducing the oxygen content changes the composition 
of the combustion products, and therefore conserving heat release on a 
unit mass basis will lead to different pressure rises In the two cases. 

For the purposes of this study It was chosen to conserve heat release 
per unit mol of combustion product. To a first order this should lead 
to the same pressure rise, and should create the same wave pattern and 
flow field for the 2 cases. This Is considered to be of fundamental 
importance for scramjet ducts, which can only be analysed correctly In 
the light of a full wave capturing treatment. Ref. 8. 


The heat release per mol of combustion products is given by 


r h c (i ♦ 


8a h f 


+ 15 «» 


H (2 - a4>* - ^ + a) + hk 


28X 


(5) 


From this the modified oxygen mass fraction for the shock tube test gas 
may be computed to be: 


X = 32 

28 { B( * * SiJ - *-«<'-*•>} 


( 6 ) 


Where A - 24 >~ 4 >* - 4/28 and B * 1 ♦ * 15 ^ ^ 


Note that this mass fraction includes the total oxygen content in all 
its forms. It is sensitive to 0 and NO concentrations, which are 
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functions of flow enthalpy, and also equivalence ratio. The test gas 
mixture Bust therefore be specifically targeted to precise test 
conditions, which will require careful planning for any experimental 
program. 

It is not possible to conserve all the flow parameters when comparing 
the 'real* and 'ideal' facilities. For Instance, the molecular weights 
of the two test gases will be different, and therefore if flow velocity 
and temperature are matched then the Mach number will have to be 
different. 

The approach of this paper is to identify certain ‘primary flow 
parameters' which it is thought important to conserve, and small 
discrepancies will be required in the other parameters in order to 
achieve this. 


The primary flow parameters are as listed below. 


(1) Total Enthalpy This is an important parameter as it relates 

directly to flight speed. Furthermore, it may be calculated with 
reasonable accuracy from shock speed and pressure measurements. 
Energy stored in heats of formation for shock tunnel flow means 
that the balance between kinetic, thermal and chemical storage 
must change. 

(2) Pitot Pressure. This relates directly to the free stream for 

scram Jet flight path evaluation and high altitude performance. 
Can easily be measured or calculated in most situations. 

(3) Static Temperature. This is not normally conserved in 

hypersonic aerodynamic testing. Mach number Independence applies 
behind strong shocks, and the free stream temperature is allowed 
to rise so that the total enthalpy can be achieved at a lower Mach 
number. However, when the combustion chamber alone is modelled, 
the flow Mach number can be matched quite closely and real static 
temperature levels can be used. Temperature is so important in 
combustion processes that it is preferred to conserve temperature 
at the expense of static pressure. 

The total enthalpy restriction means that once a static temperature and 
gas composition have been selected, then the velocity is prescribed. 
The pitot pressure restriction then determines the density, and hence 
the static pressure level. Static temperature and pressure can both be 
matched, but only at the expense of pitot pressure. 

When analysing premixed combustion processes which are meant to 
represent the combustion of fuel and air jets which were initially 
separated, it would seem appropriate to calculate an equivalent 
premixed flow conserving energy, momentum and mass fluxes. However, it 
has been found in the past, that better agreement with experiment is 
obtained by using the flow properties of the air, rather than the 
mixture in the ID analysis. This may be because in a reacting/mixing 
process combustion is expected to start in a region with a local 
equivalence ratio of ~ 0.2, and the flow properties here would be 
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closer to the free stream values. 

In tMs analysis the flow properties were calculated using the air 
fluxes only, and the hydrogen was added as a diluent at the calculate 
temperature and pressure. At this stage It Is not l -port ant. « 
romDarlson Is between analytical and numerical procedures, but It will 
have to Ee lo^edat more* closely when experimental comparisons are 

made. 

The effect of the changes in oxygen level were assessed using a ID 
Semlcal kinetics code. Ref. 3. and a 28 reaction hydrogen co»busUon 
. In table 4 The flow was constrained to follow a 

=£s 

composition was being approached. 


Results 


Reoresentatlve shock tunnel test conditions were chosen at two 
enthalpies for the evaluation, having total dissociated oxygen mass 
factions (including 0 and NO) of 0.35 and 0.63. The flow properties 
and gas composition are presented In table 1 and were calculated 
according to the procedure outlined In Ref. 2. These are established 
operating conditions using air, and were obtained using a contoured 
nozzle with an area ratio of 109. 

When eodlf led conditions were postulated » lth . ‘ 

the values of « and a_ were assumed to remain constant for ease oi 
computation. To use this analysis In conjunction with experiments ^ 
flow conditions would have to be recalculated using experimental data 

from the new test gas. 

At the 2 enthalpies the analysis was performed for equivalence ratios 
of 0.5, 1 and 2 both for the shock tube and for ideal gas flows. 

The modified flow conditions were calculated as described In the 
^vlous JoVtlon snd ,r. dlspl.yod In tsblos 2 jnd 3 Hydrogen «s 
then added to the air flow to give the target equivalence ratio, 
maintaining free stream pressure, temperature and velocity. 

The predicted pressure-distance profiles are shown In Figs, 
the 13 mj/kg condition. The results show that the shock tunnel with 
reduced oxygen concentration gives very good agreement with the ideal 
tunnel, except for the ^ * 0.5 case. 

Inspection of the species concentration axial profiles Indicated that 
the^oxygen atoms were being preferentially consumed In favour of the 
molecules for equivalence ratios of less than 1. The oxygen adjustment 
tad Sen calculated assuming that the 0. 02 and NO would be conned in 
equal proportions when the hydrogen content was not enough for full 
combustion. This leads to underestimation of the heat release an 
overestimation of the number of mols of combustion products. The latter 
arises because although the same amount of water Is produced, more mols 
of 1 0 are required to produce It than would be needed for combustion 
from the 02 molecules. The net result Is overprediction of the 
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The heat release was then recomputed assuming the 0 and HO burn, 
preferentially to the 02 molecules. The results of this analysis Is 
shown as the diagonal crosses In Fig. 1. and Is In reasonably good 
agreement with' the Ideal flow predictions. 

The results for the 18. 6 mj/kg enthalpy case are presented in Fig. 4—6. 
For the + * 0.5 case, Fig. 4, the oxygen content was calculated 
assuming combustion was only from the monatomic 0 form because the « 
value was bigger than 0.5. The results are shown to give reasonable 
agreement with the Ideal case, but It Is not as close as for the lower 
enthalpy. 

One reason for this Is possibly the lower pressure levels required In 
the Ideal flow In order to conserve the pu product. The temperature 
levels at this enthalpy are very high, ~ 3000 K after combustion, and 
heat release will be reduced due to dissociation of combustion 
products. Equilibrium composition is highly sensitive to static 
pressure level, and to check this the Ideal gas computations were 
repeated with the pressure levels matched to the corresponding shock 
tunnel cases. The results of this are shown as the diagonal crosses In 
Figs. 4-6, and they do show a reduction In the discrepancy between 
modified shock tunnel and Ideal flow prof lies. This correction Is 
achieved at the sacrifice of the pu parameter by about 15%. 

At an equivalence ratio of 1, Fig. 5, there Is significant disagreement 
between the 2 flows, as shown by the open and closed squares. 
Stoichiometric mixtures produce the hottest combustion temperatures, 
and are the most subject to dissociation. The effects of dissociation 
have been Included In the analytical treatment by using a reduced value 
for the heat of combustion from molecular oxygen, h c . For most of the 
cases considered this has proven to give a good first order estimate of 
the eventual heat release. It would appear, however, that at the 
hottest case considered # » 1, H. * 18.6 mj/kg, the assumption Is 
beginning to break down. The effective value of h,. Is being reduced, 
and with the heats of formation remaining constant the effect of the 
dissociated oxygen content becomes more significant. Intuitively, and 
by inspection of Eq. 6, It can be seen that the oxygen levels must be 
reduced even further at the conditions of highest heat release In order 
to match the ideal flow. At this point, the analysis has passed beyond 
the scope of a simple heat balance calculation and a trial and error 
approach must be made, using the ID reaction code, to find the level of 
free stream oxygen concentration which produces the required result. 
The calculated gas composition would be used when experiments were 
required for that condition. 

PROPOSALS FOR FUTURE WORK 

This paper Justifies the use of a heat balance approach for tailoring 
the oxygen concentration of the test gas to match the combustion heat 
release achieved In different test facilities. Premixed one 
dimensional flow Is assumed, and preliminary experiments should 
preferably approximate this condition as closely as possible. 

Shock induced ignition rigs, as used In Refs. 4 and 5 come close to 
achieving this. The fuel is Injected into a hypersonic flow at a low 

pressure where it can mix, but cannot burn. Ignition Is subsequently 
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Initiated by the passage of an oblique shock through the partially 

premixed layer. Experiments using the same apparatus, and adjusting 
oxygen level should give some indication as to the validity of the 
concept. 

As mentioned In the Introduction, the heat release parameter Is only 
one way In which the oxygen composition influences combustion. A more 
complete 2 dimensional approach Is required to fully analyse the mixing 
and combustion processes. It is intended that this will be done 
firstly by Beans of an existing code. Ref. 6. This will be followed by 
a series of experiments to see If the macroscopic features of the 

supersonic flow are unduly Influenced by the changes. That Is, It Is 
Important that the same basic flow pattern Is being produced in both 
cases. 

Conclusive proof of the validity of the technique can only come through 
comparison of the same experiment In different facilities. The 

expansion tube at GASL, Ref. 7, and the reflected shock tunnel T4 of 

The University of Queensland, Ref. 2 provide the potential for doing 
this. The expansion tube can produce the same enthalpy as the shock 
tunnel, but with a lower level of free stream dissociation. Parallel 
experiments are planned for the 2 facilities in 1990 to confirm that 
the same pressure/distance profiles can be produced In both cases. 


CONCLUSIONS 

Based on a 1 dimensional analysis It Is shown that the effects of free 
stream dissociation In shock tunnels can be accounted for in scramjet 
combustors. A simple heat balance approach gives adequate correction 
for most of the conditions considered. At the highest enthalpy 
examined, at an equivalence ratio of 1, the heat balance approach 

was found to be Inadequate, and a numerical analysis would be needed to 
match pressure profiles. Pressure rise downstream In the duct was used 
as the criteria to evaluate the effectiveness of the corrections, which 
consist of Baking adjustments to the oxygen content of the test gas. 
The effect of changing oxygen concentration on ignition and reaction 

delays, mixing and other flow properties was not considered. 
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Intake normalised 
static pressure 


T4 reduced oxygen content 
o preferential radical 
combustion 

_ T4 Simula tKxv normal 
oxygen content 



X<cm> 


Fig 1 Hs=12.98 MJ/kg, phl=0.5, p/x. 

O * -* 

COnTflT* 


Intake normalised 
static pressure 



Fig Z Hs*12.98 MJ/kg, pN«t p/x. 


Intake normalised 
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jarif l Representative Shock Tunnel Test Condltlgns fgr A ll- 


n nu 

Flight Sneed 5.1 Km/S 

Flight Speed 6.1 Kl 

H, (mj/kg) 

13 

18.6 

T (K) 

1874 

2520 

pu 2 (Kpa ) 

823 

847 

p (Kpa) 

23.2 

24.3 

p (Kg/m 3 ) 

0.0412 

0. 03026 

u (m/s) 

4467 

5237 

M 

5.25 

5. 18 

Gas ComDOSltlon 

(mass fractions) 


a 

0.203 

0.532 

«* 

0.145 

0. 0988 

a ♦ «• 

0.348 

0.63 

02 

0. 145 

0.08192 

0 

0.0451 

0. 1183 

N2 

0. 7493 

0.758 

NO 

0.0603 

0.412 
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TABLE 2 Flow Condition s Before Mixing. 



H. ■ 13 MJ/Kg; p3 - 823 Kpa. 1 

‘ ■ 1874 K 



I PEAL FLOW 

Unmodified 

Air 

Reduced Oxygen 

Preferential 





*0’ burn 


4 « 0.5 x 

0.2223 

0. 1636 

0.128 

0.2223 

o 2 

0. 145 

0. 1067 

0. 0835 

0.2223 

0 

0.0451 

0. 0332 

0.026 

0 

N 2 

0. 7493 

0.8157 

0.8558 

0.7777 

NO 

0.0603 

0.0444 

0.0347 

0 

P (Kpa) 

23.2 

22.7 

22.3 

20.5 

u (m/s) 

4467 

4514 

4548 

4662 

p (Kg/m 3 ) 

0.0412 

0. 04041 

0. 03982 

0. 03789 

M 

5.25 

5.26 

5.3 

5.55 

4 - 1 x 

0.2223 

0. 1573 


0.2223 

o 2 

0.145 

0. 1027 


0.2223 

0 

0.0451 

0.0319 


0 

N 2 

0. 7493 

0. 8227 


0.7777 

NO 

0.0603 

0.0427 


0 

P (Kpa) 

23.2 

22.6 


20.5 

u (m/s) 

4467 

4519 


4662 

p (Kg/m 3 ) 

0.0412 

0. 04032 


0. 03789 

M 

5.25 

5.26 


5.55 

4 = 2 x 

0.2223 

0. 1434 


0.2223 

o 2 

0. 145 

0.0935 


0.2223 

0 

0.0451 

0. 02909 


0 

n 2 

0.7493 

0.8385 


0.7777 

NO 

0.0603 

0.0389 


0 

P (Kpa) 

23.2 

22.6 


20.5 

u (m/s) 

4467 

4519 


4662 

p (Kg/m 3 ) 

0. 0412 

0.04032 


0. 03789 

M 

5.25 

5.26 


5.55 
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TABLE 3 Flow Condi tl ona Before Mixing 


PRIMARY FLOW PARAMETERS : H. ■ 18.6 KJ/Kg; pu 2 - 847 Kpa. T » 2520 K 



IPEALJUPW 

Unmodified 

Air 

Reduced Oxygen 

Preferential 
*0’ burn 


4 ■ 0. 5 x 

0.2223 

0. 1238 

0.0911 

0.2223 

Oa 

0.08192 

0.0456 

0. 03358 

0.2223 

0 

0. 1183 

0.0659 

0. 0485 

0 

N 2 

0.758 

0. 8655 

0.901 

0.7777 

NO 

0.0412 

0. 02296 

0.0169 

0 

P (Kpa) 

24.3 

23.6 

22.9 

20.4 

u (m/s) 

5237 

5284 

5339 

5489 

p (Kg/m 3 ) 

0. 03026 

0. 03032 

0. 0297 

0. 0281 

M 

5. 18 

5.22 

5.31 

5.66 

# - 1 X 

0.2223 

0.1155 


0. 2223 

Oa 

0.08192 

0.0426 


0. 2223 

0 

0.1183 

0.06151 


0 

*2 

0.758 

0.8745 


0.7777 

NO 

0.0412 

0.0214 


0 

P (Kpa) 

24.3 

23.5 


20.4 

u (m/s) 

5237 

5298 


5489 

p (Kg/m 3 ) 

0. 03026 

0.03016 


0. 02819 

M 

5. 18 

5.24 


5.66 

# • 2 x 

0. 2223 

0.0996 


0.2223 

Oa 

0.08192 

0. 03672 


0.2223 

o 

0. 1183 

0.053 


0 

*2 

0.758 

0. 8918 


0.7777 

NO 

0.0412 

0.0185 


0 

P (Kpa) 

24.3 

23.1 


20.4 

u (m/s) 

5237 

5324 


5489 

p (Kg/m 3 ) 

0. 03026 

, 0.02987 


0.02819 

M 

5. 18 

5.28 


5.66 
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ORIGINAL PAQ£ ie 

OF POOR QtMLrTY 



Kpa, T • 1874 K 




0.2223 

0.145 

0.0451 

0.7493 

0.0603 

23.2 

4467 

0.0412 


0. 1636 
0.1067 
0.0332 
0.8157 
0.0444 
22.7 
4514 
0.04041 


0. 128 
0.0835 
0.026 
0.8558 
0.0347 
22.3 
4548 
0.03982 


0.2223 

0.2223 

0 

0.7777 

0 

20.5 

4662 

0.03789 









TARi r a Flow Condition! Before Mixing 


PRIMARY FLOW PARAMETERS : H, - 18.6 KJ/Kg; p u 2 - 847 Kpa. T - 2520 K 




IPEAlr FLOS 

Unmodified Air 

Reduced Oxygen 

Preferential 
•O’ burn 


f « 0.5 x 

0.2223 

0.1238 

0.0911 

0.2223 

o 2 

0.08192 

0.0456 

0.03358 

0.2223 

0 

0.1183 

0.0659 

0. 0485 

0 

N 2 

0.758 

0. 8655 

0.901 

0. 7777 

NO 

0. 0412 

0. 02296 

0. 0169 

0 

P (Kpa) 

24.3 

23.6 

22.9 

20.4 

u (m/s) 

5237 

5284 

5339 

5489 

p (Kg/m 3 ) 

0. 03026 

0.03032 

0. 0297 

0.0281 

M 

5.18 

5.22 

5.31 

5.66 

♦ - 1 * 

0.2223 

0.1155 


0.2223 

o 2 

0. 08192 

0.0426 


0.2223 

0 

0. 1183 

0.06151 ' 


0 

n 2 

0.758 

0. 8745 


0.7777 

NO 

0. 0412 

0.0214 


0 

P (Kpa) 

24.3 

23.5 


20.4 

u (m/s) 

5237 

5298 


5489 

p (Kg/m 3 ) 

0. 03026 

0. 03016 


0. 02819 

M 

5.18 

5.24 


5.66 

4 B 2 x 

0.2223 

0.0996 


0.2223 

° 2 

0.08192 

0. 03672 


0.2223 

0 

0. 1183 

0.053 


0 

N 2 

0.758 

0.8918 


0.7777 

NO 

0.0412 

0.0185 


0 

P (Kpa) 

24.3 

23.1 


20.4 

u (m/s) 

5237 

5324 


5489 

p (Kg/m 3 ) 

0. 0302< 

> 0. 02987 


0. 02819 

M 

5.18 

5.28 


5.66 
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THE USE OF SILANE AS A FUEL ADDITIVE FOR HYPERSONIC THRUST PRODUCTION 


R.G. Morgan, Department of Mechanical Engineering, 
The University of Queensland 


ABSTRACT. 


Scramjet propulsive ducts operate by converting some of the flight free 
stream kinetic energy to heat In order to initiate combustion. As the 
flight speed Increases, the available kinetic energy which can be used 
for heating also increases. However, the required temperature range for 
efficient combustion does not change much with flight speed. 
Temperatures are limited at the lower end by Ignition and reaction 
characteristics. and at the upper end by combustion product 

dissociation. It is realistic to expect combustion chamber intake 

temperatures to be within the range of 800 to 1800 K throughout the 
whole flight envelope. This means that at the higher flight Mach 
numbers, the combustion chamber Mach number must be in the hypersonic 
regime in order to avoid overheating the flow. The static pressure 
levels associated with the hypersonic flow will be low, and this will 
adversely effect combustion and heat release characteristics. This 
paper reports on the use of silane as an additive to enhance combustion 
in shock tunnel tests of a hypersonic scramjet. Silane, which is 
spontaneously combustible at room temperature, was mixed with hydrogen 
at 20 */. by volume. The mixture has been used before in a constant area 
duct configuration, and for this work a 2 dimensional nozzle was added 
to produce thrust . 


The silane was found to produce vigorous combustion, and significant 
thrust, but when specific impulse was calculated, on a per unit mass 
flow rate basis, it did not perform better than pure hydrogen over the 
same enthalpy range but at a higher static pressure. This is attributed 
to the high molecular weight of the silane molecules, and the associated 
high mass flow rate of fuel required to burn the oxygen. At the static 
pressure levels for which the tests were performed, pure hydrogen did 
not burn at all within the confines of the duct. In conclusion, the 
silane was useful for performing hypersonic combustion studies under 
test conditions where hydrogen combustion was not possible. 
Additionally, the tests indicate that combustion can be enhanced, but 
that the specific impulse will not exceed that of hydrogen in regions 
where the hydrogen itself can be made to burn well. 
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1. INTRODUCTION. 


Experiments In hypersonic combustion of hydrogen have been performed in 
The University of Queensland shock tunnel, T4 U) . Whilst the results 
indicate that some mixing and combustion does occur. ‘ 2 >. the performance 
when a thrust nozzle is attached does not approach that achieved in a 
supersonic duct in the shock tunnel T3, <3) . at the Australian National 
University. The reason for this is not fully understood at present. 

The hypersonic flow is created at the expense of free stream static 
pressure, and pressure scaling tests in a supersonic duct in T3 ( * have 
shown combustion to be very pressure sensitive. 


The supersonic duct typically operates at static pressure levels 
~ 150 kpa, but when the pressure level is reduced below 50 kpa 
combustion is inhibited. The hypersonic duct has an intake pressure 
level of ~ 20 kpa, and if a direct comparison with the supersonic tests 
applies, it might be expected to be adversely effected by the low 
pressures. 

There is some evidence to suggest that low static pressure levels are 
the primary cause of the poor combustion. When the scramjet was 
configured as a 51 x 25 mm section constant area duct, boundary layer 
growth on the internal walls caused significant pressure rise, (2) and 
Fig. 1. When hydrogen was injected at this condition it burned fairly 
readily. Fig. 2. However, the boundary layer compressive effects would 
have reduced the free stream Mach number to ~ 4.5, and it could not 
truly be called hypersonic combustion. 
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When the duct section was increased to 50 x 51 w the compressive effect 
of the boundary lawyer was reduced, the Fuel Off pressure levels were 
lower and the flame was extinguished. Fig. 3. When the pressure levels 
In the square duct were subsequently Increased by means of a compression 
wedge. Fig. 4, combustion effects were again evident In the duct. 

Fig. 5. 

Both the compression wedge, and compression waves off the developing 
boundary layer. Increase the static temperature as well as static 
pressure levels. However the static temperature levels may be held 
constant by adjusting the stagnation enthalpy. This leaves only Mach 
number and static pressure level as the parameters which are different 
between conditions where combustion does, and does not occur. 

The pressure may Influence combustion through lt*s effect on density, in 
two ways. Firstly, reaction rates will be enhanced, to a degree 
depending on the order of the reaction mechanisms involved, due to 
Increased collision rates. In a flame where mixing is expected to occur 
rapidly, this will be the major heat release controlling factor at a 
given free stream static temperature. 

Secondly, the Increased Reynolds number may increase the rate of 
entrainment between the fuel and air layers, which, if reaction rates 
are fast will Increase the overall rate of heat release. The T3 
supersonic duct tests were performed In the free stream unit Reynolds 
number range of 1 to 2 10 7 m' 1 . whilst the T4 hypersonic runs were In 
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the range 0.5 to 1 *10 7 ■*». This In itself nay be enough to explain 
the increased combustion in the supersonic flow. This effect could be 
partially checked by raising the pressure levels of the T4 tunnel, and 
this will be done when the facility Is capable of doing so. However, 
any change would be coupled with Increased reaction rates and may be 

hard to Interpret. 

The hypersonic conditions where combustion occurs have, as well as 
higher static pressure levels, lower Mach numbers than for test cases 
where no combustion Is observed. The effect of Mach number on the rate 
of development of a mixing layer Is not fully understood. <». however It 

may be a significant factor. 

The purpose of these experiments was to partially decouple the effects 
of chemical kinetics by the use of a combustion enhancing additive In 
the fuel. Silane at 207. by volume was mixed with hydrogen and injected 
at equivalence ratios of up to 1. based on the consumption of all fuel 
and oxygen components. The resulting pressure traces Indicated 
substantial combustion and heat release. However, this cannot be solely 
attributed to increased reaction rates, as the silane mixture has a 
molecular weight of 8. and has to be injected at 2.5 times the mass flow 
rate of hydrogen for the same equivalence ratio. This leads to a slower 
jet. with a large momentum deficit which may enhance the mixing. 

When a nozzle was added to the combustion chamber, significant thrust 
was measured for the fuel on condition. The thrust levels when 
normalised against reservoir stagnation pressure were constant over the 
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usable test Use of the shock tunnel. Indicating that steady flow and 
combustion had been achieved by the us. of the silane additive. This 
enabled thrust and specific Impulse measurements to be mad. In a 
hypersonic flow where current T4 pressure Hal tat Ions preclude 

combustion of hydrogen. 


2. EXPERIMENTS . 

The initial experiments were performed In a constant area duct 
configuration to confirm that significant combustion could be achieved, 
and to determine at what axial location the thrust nozzle should start. 
* sample set of Fuel On and Fuel off pressure profiles are shown In 
Fig. 6 for a nominal stagnation enthalpy of 6.8 MJ/kg. 


In the light of these results the scramjet was re-conflgured for thrust 
production as shown In Fig. 7. The thrust surface diverged at an angle 
of 150 fro » the incoming flow, and started at a distance of 225 mm 
downstream of Injection. Injection was fro. a two dimensional central 
injector with a throat of 1.6 ■» height. The Injector spanned the full 
width of the duct. Fuel mass flow rate was controlled by adjusting the 
stagnation pressure of the fuel reservoir, which was maintained at room 
temperature. Changes In equivalence r.tlo therefore -ere coupled to 
changes In fuel Mach number, velocity and temperature, which had an 
unknown effect on the rate of fuel air mixing. 


The thrust surface on the duct was instrumented with quartz transient 
pressure transducers, and thrust was calculated assuming 2 -dimensional 
flow with linear interpolation between data points. 
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3. RESULTS. 


Shock tunnels are Inherently unsteady facilities, and care must be taken 
when trying to draw conclusions based on steady combustion assumptions 
from transient data. The first step when analysing shock tunnel data 
must be to confirm that steady flow has been established. As the 
stagnation pressure decays during the test time of the tunnel, steady 
Mach number will only be maintained If the pressure ratio across the 
nozzle attached to the shock tube Is constant. 

In Fig. 8 the static pressure/tlme history Is shown for a typical fuel 
on shot for a transducer on the thrust surface. It should be noted that 
for the duration of the steady flow period, although the stagnation 
pressure changes by ± 15*/.. the normalized static pressure level stays 
sensibly constant. Similarly, the integrated thrust when normalised In 
the gan>«» manner Is nearly constant, Indicating that the macroscopic flow 
field Is steady. 

If therefore the nozzle exit conditions are calculated on the assumption 
of a steady expansion from the stagnation region, then the scram jet 
tests are conducted In an intake flow field with the following 
approximate variation in properties from the nominal conditions: 


Flow duration 1 ms 

Stagnation pressure variation ± 1554 
Stagnation enthalpy variation ± 454 
Static pressure variation ± 1514 
Static temperature variation ± 414 
Mach number constant 

Velocity variation ± 214 

Density variation ± 10% 

Momentum flux variation ± 15% 

Equivalence ratio variation t 13% 
Reynolds number variation ± 10% 
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[These quantities estimated at a stagnation enthalpy at 12.5 MJ/kg from 
run 874]. 

To assess the significance of this unsteadiness on thrust production. It 
Is necessary to estimate the decay of momentum within the control volume 
during the test period. The total momentum stored within the control 
volume is 0.135 kg. ms’* for the 12. 5 MJ/kg case, and the decay rate 
during the 1 ms of test time therefore represents a force of ~ 35 N. 

It can be seen from the thrust balance Illustrated schematically In 
Fig. 11 for a 2 -dimensional duct of unit width that this momentum decay 
will detract directly from the measured thrust on the thrust surface. 
Moreover, this value of 35 N Is a significant fraction of the fuel on 
steady thrust levels of ~ 50 N. 

However, when comparing Fuel On and Fuel Off thrust levels it is 
estimated that the difference between the two due to unsteady momentum 
effects reduces to a second order error. This is because, due to the 
lack of upstream Influence in supersonic flow, the intake time history 
of pressure and momentum flux Is the same for both Fuel On and Fuel Off 
cases. The new thrust and specific Impulse calculated by subtracting 
Fuel Off from Fuel On values will therefore be primarily due to the 
presence of the fuel jet as governed by steady flow and steady 

combustion effects. 

Also, the Internal momentum decay can be considered to detract directly 
from the thrust as indicated in Fig. 11 only when comparing steady and 
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unsteady flows with the same Incoming and outgoing momentum flux and 
pressure integrals. In practice with the onset of unsteady flow the 
momentum decay term will be partially met by changes to the incoming and 
outgoing integrals. Take for example the case of a constant area duct 
In unsteady flow. No thrust is produced, and a rise in downstream 
pressure balances the momentum decay. The more the nozzle departs from 
a constant area duct, the more sensitive the thrust production will be 
to unsteady effects. 

In Fig. 9 scatter bars are shown illustrating how the normalised 
pressure/distance profiles change over the nominal test period. The 
steadiness of the profiles Indicates that the macroscopic features of 
the flow, such as wave propagation and thrust production are constant. 
In Fig. 10 the normalised integrated thrust/tlme profiles are shown for 
Fuel On and Fuel Off conditions. The constant levels of the two traces 
adds credibility to the assumption that the unsteady thrust decrement is 
the same for Fuel On and Fuel Off runs. 

In Figs. 12-17 the pressure profiles are shown for 205C silane injection 
in the enthalpy range 4.46 to 17.1 MJ/kg. For all of the conditions 
tested except 4.46 Mj/kg strong combustion affects were observed in 
terms of increased pressure on the thrust surface. At the lower 
enthalpy the static temperature was ~ 500 K which appears to be too low 
for silane combustion at this pressure level. 

In Fig. 18 the effects of injecting helium into a flow at 6.81 Mj/kg are 
shown. It is seen to produce negligible pressure rise over the Fuel Off 
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case. 


This suggests that the pressure rises observed when Injecting the 
silane mixture are due to combustion rather than the intrusive presence 
of an inert Jet. The molecular weight of helium was 4, compared with 8 
for the silane mixture, but when coupled with the evidence of the low 
enthalpy, no pressure rise test. Fig. 12. it seems certain that the 
silane mixture was burning at the high enthalpies. 

The measured net thrust with silane injection was high, but when the 
specific impulse was computed. Fig. 19. it did not compare favourably 
with the hydrogen results from (2) . shown in Fig. 20. This is 
attributed to the high molecular weight and mass flow requirements of 
the heavy silane molecules. However, the pressure level of these tests 
was lower than for <2) , and hydrogen would not burn at all under these 

conditions. 

4 . CONCLUSIONS . 

The use of silane additive at 20% induced vigorous combustion at 
conditions where hydrogen was not able to burn. This tends to suggest 
that kinetics, rather than mixing are limiting heat release in the 
hypersonic duct. However, this contradicts to some extent the findings 
of (6> whose numerical study suggests that mixing will be the limiting 
factor. A possible reason for enhanced mixing with silane might be the 
lower velocity, high density fuel jet required for stoichiometric 
combustion. 

Alternatively, rapid vogorous combustion might induce turbulence which 
enables the fuel to mix over a wide spread region of the duct. 
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An estimate of the extent of the mixing layer may be made by Inspection 
of the fuel on pressure traces, such as Fig. 9 and Figs. 13-17. As 
discussed In Ref. (9), thrust Is transmitted to the thrust surface by 
two separate mechanisms. These mechanisms are (1) direct propagation of 
waves from the combustion zone to the walls, and l2) the Interaction of 
the corner expansion fan with the Mach number gradients In the flow. 
Waves from mechanism 1 are evident upstream of the thrust surface, 
whilst the onset of mechanism 2 Is clearly shown from the pressure rise 
downstream of the corner. By following waves back upstream from the 
start of the thrust surface pressure rise. It Is possible to locate 
where the leading edge of the expansion fan first starts producing 
compression waves. This Identifies the point where perturbations first 
occur In the free stream Mach number profile, which may be taken to 
represent the limits of the combustion zone. 

The results of these calculations are shown in Fig. 21. The Mach number 
before the expansion fan Is estimated by assuming that the freestream 
flow upstream of the corner Is compressed isentropically from the 
nominal Intake conditions to the locally measured static pressure. The 
compression waves are assumed to propagate at the local Mach angle. Run 
874. equivalence ratio of 1 at a stagnation enthalpy of 12.5 Mj/kg was 
used for the calculations. The first compression waves are noticed 
Incident on the thrust surface at a distance of approximately 70 mm from 
the corner. 

The results of this treatment Indicate that the combustion zone extends 
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to within approximately 3. 5 am of the wall. Given that 
boundary layer at the wetted length of - 400 mm will be of the order of 
2 thl. Implies that the mixing layer has effectively propagated 

across the whole duct. If this result can later he confirmed by optical 
techniques, and by a more detailed wav. analysis of the pressure 
profiles, then the prospects for hypersonic combustion are good. 


Silane with vigorous combustion produces comparable performance to 
hydrogen with partial combustion, as Indicated by small pressure rises 
associated with Injection. This suggests that complete combustion of 
hydrogen, whether Induced by Increased pressure levels or though mixing 
enhancement techniques, can potentially perform very well In 
hypersonic duct. Alternatively, the us. of a lower concentration of 
silane, as tested by "> might be used to Improve specific Impulse. 


Ecramjet performance In a hypersonic duct does not yet compare with that 
of combustion In a supersonic duct at higher pressure levels. However. 
It is hoped that increased pressure simulation capability will be able 


to improve hypersonic performance. 
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TABLE 1 RUN MANIFEST IN 15° THRUST MODE 


Stagnation 
Enthalpy (H.) 
MJ/kg 

Run 

No. 

Injection 

Pressure 

kPa 

Equivalence Ratio 

4.46 

875 

325 S,H 4 mix 

0. 99 

4.46 

883 

Fuel Off 

0 

6.81 

882 

0 

0 

6.81 

871 

0 

0 

6.81 

876 

202 S t H 4 mix 

0.66 

6.81 

870 

393 S,H 4 mix 

1.28 

6.81 

884 

380 H e 


8.5 

880 

0 

0 

8.5 

879 

217 S s H 4 mix 

0.82 

8.5 

873 

373 S,H 4 mix 

1.41 

8.5 

885 

473 H e 

— 

10 

868 

0 

0 

10 

877 

209 S t H 4 mix 

0. 88 

10 

867 

410 S(H 4 mix 

1.73 

10 

886 

382 H e 

— 

12.5 

881 

0 

0 

12.5 

874 

373 S t H 4 mix 

1.77 

12.5 

887 

366 H e 


17. 1 

889 

0 

0 

17. 1 

888 

350 S,H 4 mix 

2. 1 

17. 1 

891 

399 H e 



(Injector calibration 6.6 10“ 5 Kg s _1 kg -1 


TABLE 2 - NOMINAL TEST CONDITIONS 


H s 

Static 

Static 

density 

Mach 

Velocity 

Intake 

MJ/kg 

Temp. K 

pressure 

Kg/m 3 

No. 

m/s 

Oxygen 



kPa 




mass flow 

4.46 

529 

7.54 

0.0493 

6.09 

2820 

.0695 

6.81 

888 

9.72 

0.0379 

5.73 

3430 

.065 

8.5 

1120 

9.8 

0. 0303 

5.5 


.056 

10 

1350 

10. 1 

0.025 

5.43 


.05 

12.5 

1570 

10 

0.0205 

5.43 

4350 

.0446 

17. 1 

1980 

8.98 

0.014 

5.44 

5000 

.035 


72 














LIST OF FIGURES 


1 . Fuel Off pressure profiles for 25 x 51 mm 2 constant area duct. 
H. * 8.4 MJ/kg. 

2 H Fuel On, Fuel Off and helium injection comparisons for 
2§ x 51 mm 2 constant area duct, H. = 8. 4 Mj/kg. 

3 . Fuel On and Fuel Off pressure profiles for hydrogen injection in 
49, 5 x 51 mm 2 constant area duct, H, 8.39 MJ/kg. 


4 . Schematic of square duct with 5° wedge. 

5 . Hydrogen Fuel On and Fuel Off pressure profiles for 49.5 x 51 mm 2 
duct with 5° wedge. H, = 8.4 MJ/kg. 

6. Fuel Off and Fuel On profiles for S,H 4 20% mixture in 49.5 x 51 mm 2 
constant area duct. Hj = 6.81 Mj/kg. 


7. Schematic of 49.5 x 51 mm 2 duct with 15° thrust surface. 

8. Pressure and thrust/time profiles for silane mixture injection into 
49.5 x 51 mm 2 duct with thrust surface. H. = 12.5 Mj/kg. 


9. 


10 . 


Scatter bars for pressure profiles during steady 
mixture in 49.5 x 51 mm 2 duct with 
H. * 12.5 Mj/kg. 


flow period. SjH 4 
thrust surface. 


Normalised thrust comparisons between Fuel Off and S t H 4 mixture 
Fuel On for 49.5 x 51 mm 2 duct with thrust surface. 

H. * 12.5 MJ/kg. 


11. Thrust balance for unsteady flow. 

12-17P/X profiles for fuel on (silane mixture) and Fuel Off conditions 
in 49.5 x 51 mm 2 duct with thrust ramp, enthalpies as indicated 

12. H. * 4.46 Mj/kg 

13. H, = 6.81 Mj/kg 

14. H, = 8.5 Mj/kg 

15. H. = 10 Mj/kg 

16. H, = 12.5 Mj/kg 

17. H. = 17.1 Mj/kg 

18. H, = 6.81 Mj/kg, comparison of Fuel Off and helium injection. 

19. Specific impulse/stagnation enthalpy. S,H 4 mixture, 49.5 x 51 mm 2 
duct with thrust surface. Intake pressure = 10 kPa. 

20. Specific impulse/stagnation enthalpy for hyrogen injection into 
supersonic and hypersonic flow. 25 x 51 mm 2 duct with intake 
pressure 20 kPa. 

21. Compression wave tracking to find extent of mixing layer. 
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PRESSURE-LENGTH CORRELATIONS IN SUPERSONIC COMBUSTION 
P.A. JACOBS & R.J. STALKER 


An experiment 1 has been proposed for study of the pressure sensitivity 
of supersonic hydrogen combustion in a constant-area combustor. The 
specific aim of the experiment was to investigate the theoretical 
pressure- length correlations 2 by conducting trials in three 
geometrically similar models. The pressure in the models was to vary 
over an order of magnitude but the pressure by model - length product 
was to be held constant. 

Here we report that the largest of the three models (and in some ways 
the most difficult model) has been constructed and the initial flow 
measurements have been made. The correlation P 1 - 7 l = const and the 
length of the combustor (1.2m) determined the smallest pressure below 
which we would not get combustion in the largest of the three models. 
The duct pressure in these runs was ~ 29 kPa. 

The overall flow configuration is shown in figure 1. The shock tunnel 
is fitted with a nozzle that expands the gas to approximately M = 8 
(design area ratio = 670). The test gas is recompressed to M * 3.4 
through a pair of wedges (and their associated shocks) and is then 
allowed to enter the model. The first 113 mm of the duct expands the 
gas flow to M * 4.4 and at the end of this expansion the hydrogen fuel 
is injected from the rear of a centrally located strut. The hydrogen 
is injected parallel with the air flow and forms a wake-like mixing 
region downstream of the injector. The estimated flow conditions at 
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the Injector exit plane are given in table 1. 


The wall pressure, as a function of distance downstream of the 
injector, is shown in figure 3 for hydrogen injection into air, 

hydrogen injection into nitrogen and air flow without hydrogen 
injection. From the general trend of Increasing pressure when hydrogen 
is injected into the air flow and constant pressure otherwise, it 
appears that the hydrogen does burn at these conditions. 

(H « 9.1 MJ/kg, P duct a 29 kPa ♦ * 0.76). Hence we have something 

interesting to measure. 

The next step is to try to match the CFD codes to these pressure data 
and then use the codes to compute burning lengths and ignition lengths. 
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TABLE 1 ; Nomihal conditions In large scramjet. 
Shock tube tube data taken from shot 1410. 
diaphragm thickness = 5mm 
fill pressure =150 kPa 
fill temperature = 300K 


Shock speed = 3. 145 km/sec 

steady reservoir pressure ~ 60 MPa = P, 

ESTC (conditions after expanding reservoir to 60 MPa) 
stagnation enthalpy H s * 9. 13 MJ/kg 
temperature T B * 5747 K 

NENZF calculation of conditions at the end of the M=8 nozzle, (stop 
expansion when P p ,i ot /P. * 0.004 , assume P pnot « 0.92 pu 2 ) 
effective length of nozzle = 9.89 cm 
M = 7.12 y = 1.36 u = 3962 ms -1 

P = 0.01660 kg/m 3 P.tatic “ 3-795 kPa T, tatlc = 790 K 

flow through shocks (15° wedges) and then an isentropic expansion 
through the initial part of the duct (area ratio * 2.653) to the exit 
plane of the injector. 


Hdeet 85 4- 37 


P.tatic * 28.8 kPa 
T.tatic * 1802 K 


disassociation fraction 
through shocks) 

«o a 


(assuming the composition remains 


[0] + [NO] 

[0] + [N0J = 2 [0 2 1 


0.207 


frozen 
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Table 1 (continued): Nominal conditions in large scramjet. 


Hass flow of air 

after shocks p = 0. 1571 kg/m 3 , u = 3420 ms' 1 

intake size = 0. 014 m x 0. 100 m 

m alr = 0.7522 kg/sec 

ntQ 2 = 0. 232 M alr = 0. 1745 kg/sec 

Mass flow of hydrogen 

H 2 manifold pressure “ 250 kPa (shot 1420) 

Injector calibration constant = 6.60e-5 kg/sec /kPa 
mH 2 = 0. 0165 kg/sec 

m H 2 

Equivalence Ratio ^ = - — 

“02 


— = 0.756 
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HOT HYDROGEN INJECTION TECHNIQUE FOR SHOCK TUNNELS 
by Michael Wendt 

The effect of fuel temperature on combustion In a supersonically and 
hypersonically combusting ramjet Is currently being studied In the T4 shock 
tunnel. It is intended to not only study the temperature effect, but the 
corresponding density distribution and velocity distribution effects as well. 

A device is currently being manufactured to produce hot hydrogen for 
injection into the ramjet model. Many difficulties were encountered in 
its development, including: 

1) The ramjet model instrumentation calibration was very temperature 
sensitive and hence heat sources had to be kept well away from the test area, 
ii) To minimize losses only a short length of piping was allowable between 

the heat source and injector nozzle. 

ill) Many standard components had difficulty withstanding the required 
stagnation temperature ( 1200 K ) at the injection pressure, 
iv) Only a small amount of fuel was to be used to minimize the build 
up of unburnt hydrogen inside the Dump Tank. This corresponded to a 
maximum supply time of 20 - 40 milliseconds. 

After evaluating a number of solutions it was decided to use a Gun Tunnel. 

The gas is heated by a combination of free piston compression and entropy 
rise through shock waves. It has the advantages of: 

i) The hydrogen is only hot for a very short duration « 30 msecs) and 
hence the temperature of the vessel walls is only raised slightly. 

ii) The danger of having hot fuel in a vessel for an extended 
period of time is reduced. 

ii) No heat is transferred to the ramjet model. 

iii) Dry hydrogen is produced. 
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The layout of the gun tunnel Is shown in the diagram over leaf. It 
consists of two tubes separated by a bursting diaphragm. The high pressure 
reservoir contains nitrogen and provides the driving pressure for the 
piston. The compression tube initially contains hydrogen at low 
pressure. A light nylon piston initially sits in the compression tube Just 
in front of the primary diaphragm. The diaphragm is mechanically burst 
allowing the high pressure nitrogen to expand and drive the piston down the 
compression tube. The light piston is fast enough to cause shock waves ahead 
of it in the hydrogen gas giving an entropy increase (the entropy rise is 
small, however, due to the high speed of sound in H 2 ). The pressure in 
front of the piston increases causing the piston to decelerate near the high 
pressure end. The pressure increase also bursts the secondary diaphragm. 

The piston fluctuation dampens rapidly as pressure on both sides of the 
piston approach equality. Steady injection occurs with a near constant 
pressure on the rear of the piston. 

The tunnel design has the following notable features: 

1) The tunnel is positioned within the Dump Tank of T4 to increase the 
safety in event of explosion and to minimize the length of tubing between 
the hot hydrogen reservoir and ramjet model injector. 

ii) The entire assembly is free to reciprocate linearly to allow the 
dynamics of the device to lower the stress induced in its structure. 

iii) Bursting diaphragms are used to provide simple, low 
maintenance, reliable flow control capable of containing the high 
temperatures and pressures. 

Construction of the tunnel is due to be completed in early 1990 with the first 
ramjet experiments due in the middle of the year. 


87 




Hydrogen Gun Tunnel Layout 



HEAT RELEASE - WAVE INTERACTION PHENOMENA 
IN HYPERSONIC FLOWS 


AN OPTICAL STUDY 


Annual Report, February 1990 
by N.R. Ward 


INTRODUCTION 

The study of combustion waves in a scram jet duct is 
made very difficult due to the imperfect nature of the 
injection and mixing of fuel into the flow and the inherent 
disturbance to the flow which this causes to produce the 
desired heat release. This can be avoided by bringing about 
this heat release in another way. The flow produced by a 
free piston shock tunnel using a nitrogen test gas is highly 
dissociated but basically 'frozen' in this state due to it 
having been expanded to a temperature too low for a 
significant rate of recombination. When, however, it has its 
temperature and pressure raised passing through the system 
of shocks at the entry to the scramjet duct recombination 
progresses rapidly to bring the flow to its equilibrium 
level of dissociation. This recombination results in heat 
release, simulating quite closely the heat release during 
combustion. Thus by studying the waves generated in this way 
in a scr am jet model in nitrogen flow in the T4 shock tunnel 
it is hoped that much can be learned about 'ideal' 
combustion and the combustion waves it generates. 

It was decided to study these waves using optical 
techniques in order to keep disturbance of the flow to a 
minimum and collect as much information at each shot of the 
shock tunnel as possible. The optical system decided upon 
was the Differential or Schlieren Interferometer. 

Since the analysis of these flows was to be done using 
optical methods development of a short exposure shutter and 
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associated triggering system was, in 1988, carried out. 
Combined with a polaroid camera this was used to take 
photographs of test flow over models. Some of the 
information gathered was then used to check experimental T4 
test conditions against those computationally predicted, 
showing good agreement. A computer program was then 
developed to predict post-shock recombination profiles for 
different experimental conditions, allowing for detailed 
design of test conditions and model shapes for future tests 
in the T4 shock tunnel. 


REPORT 

in preparation for experiments in the T4 shock tunnel 
it was neccessary to design and build a model of an 
open-sided scramjet duct in which generated 'combustion' 
waves could be studied freely using optical methods. It was 
decided to begin with a simple model consisting of a pair of 
wedges as an intake followed by a pair of plates as a duct, 
with simply an open path to exit. Predicted intake shock 
wave angles were calculated to allow mapping of these shocks 
and consequent design of the model components such that 
these shocks are allowed to escape the duct through gaps in 
the walls rather than having them reflect and propagate the 
length of the duct, thus confusing the wave pattern. The 
resultant model design is shown in Figure 1. This model was 

then built. 

The most critical task to be completed before the 
actual experiments with the T4 shock tunnel could begin was 
to be the design, assembly and alignment of the optical 
system. The extent of the difficulties would, however, 
exc6Gd all expectations* 

A CCD camera system capable of capturing images with 
exposure times of as low as 100 us, which had been developed 
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in this department , was need with only a *T 
probleme. Details of this systsm can be found 

<5> ' h Laser Electronics ruby laser system with He-He 
, . ,_ ser was installed on the roof of an enclosed 

alignment laser ^ on the te . t .action -dance 

° p tics b ox mou diagram of the optics 

r ^ tt^am was lead down into the optics box 

box. The wa . decided that due to the short 

using an °Ptic^ ^ ^ recorded by the optlcs . hou ld be 

teS Dossible to instantaneous, thus it was decided 

t: :rairrele b ctronic O-switch to bring the ruby laser 
pulse time down from about 1 ms to less than 50 ms. 

unfortunately the Pockells Cell component of 

..... ■"“> - 
of T4 experiments. 


FIGURE 2. SKETCH OF OPTICS BOX FROM THE SIDE 



\ opemi^s 

optics BO* 

6CARD 
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The next task was to design and build mounts for the 
optical components of the actual optical system. Once the 
components were mounted it was time to develop a procedure 
for aligning them so as to give a proper interference 
pattern and a sharp image of an object in the test section, 
and hence a good differential inter ferogram. The design that 
had been settled upon was that of a double pass differential 
interferometer, which made the alignment more difficult. The 
layout for this optical system is shown in Figure 3. 

It was found that most of the alignment was easiest 
carried out using an extra He-Ne laser in place of the 
output end of the fibre optic inside the optics box, as this 
gave a bright and non-diverging colimated beam which was 
easy to trace through the system. Also, in the course of 
alignment a number of the optical components were found to 
be of insufficiently high quality or of inappropriate 
specifications to be able to obtain an image. These were 
replaced in due course and eventually a systematic procedure 
for aligning the system was developed. This can be 
summarised as follows: 

(i) Centre alignment laser to give a well directed and 
horizontal beam. 

(ii) Work through the system of mirrors ( with other 
components removed ) adjusting each to pass the beam on to 
the centre of the next. 

(iii) Use a screen with a hole for the outgoing beam in 
it to adjust the final mirror so that the returning beam 
follows the same path as the outgoing beam. 

(iv) Position each lens in turn in the beam centred so 
as to not alter the path of the outgoing beam. 

(v) Adjust lens positions in direction of beam so that 
the focus of the beam before the concave mirror coincides 
for the incoming and outgoing beams. Readjust the final 
jiirror if neccessary so that incoming and outgoing beams 
still coincide. 
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(vi) Place polaroid and wollaston prism in the path of 
the beam near the focus of the concave mirror. 

Some adjustment of some components ( such as the cubic 
beamsplitter ) has also been found neccessary to eliminate 
stray reflections which disturb the image. This procedure 
has been found successful and good inter ferograms ( using 
the direct He-Ne beam as a source ) have been obtained. 

Currently final alignment is being carried out with the 
laser source travelling through the fibre optic. Presently 
the original multi-mode fibre optic is being replaced with 
single-mode fibre optic to eliminate severe distortion in 
the illumination pattern of the light emitted from the 
fibre. 


FUTURE 

It is expected that the differential interferometer 
will be fully operational shortly. Tests to determine how 
severely the vibrations during the operation of the shock 
tunnel affect the alignment of the optical system will then 
be carried out in order to minimise this disturbance. 

Experiments with the previously described model in the 
T4 shock tunnel will follow, using the differential 
interferometer to obtain information on the simulated 
combustion waves caused by nitrogen recombination aswell as 
other wave phenomena in the scram jet duct. 
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A STUDY OF THE WAVE DRAG IH HYPERSONIC SCRAM JETS 


Nesrln Altman 

The Mechanical Engineering Deptartment 
The University of Queensland, St. Lucia, 4067 


ABSTRACT 

A model two-dimensional flow In a convergent-divergent duct Is 
analysed. It Is found that wave drag varies with the Mach number In a 
“saw tooth" manner with increasing amplitude and period. The drag 
value reaches Its minimum when symmetrical wave distribution can be 
obtained In linear theory within the system, the maximum drag value 
coincides with Mach numbers larger than that In which the shock wave 
emanating at the leading edge of one plane completely misses the other 
plane. It is also found that the thrust varies linearly when 
combustion is simulated by adding heat without addition of mass or 
change in the compressibility factor. 

1.0 INTRODUCTION 

Wave phenomena and management of waves play an important role in 
hypersonic scramjets. In order to achieve beneficial Interaction 
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^ r f.f.*e it Is necessary to 
.. and the thrust surfaces. It 

between th result of the wave 

, . the drag produced within the system as a 

understand the drag P ^ ^ ^ flw 

tlon This understanding becomes m 
dlstrlbutl . values and Inefficient 

, tMn the duct take hypersonic values. 

Mach numbers within tne au 

flow aan.ge.ent results In Increased dr as- 

th . results obtained ,n this stud, -Ub that of linearised 
TO a validate tb. nuaerlcal aetbod used In tb. -l,sU. an 

Z LU -sentl,. tb. teaperature dependence effects 

USC d the viscosity effects are Included 

r the specific heat constants, and the vlscos 
of the spec li Kiaa-dlns" to achieve 

, tb. stud. The concept of "boundary Lifer bleeding 

— r: 

at the throat, without the addition 

v , . normally is associated with combustion, 
compressibility factor which normally 

The flow is assumed to be two-dimensional. 


•7 The Busemann Biplane . 

A tuo dlaenslonal convergent-divergent duct, for^ by pUCl ^ 

hln double -edge aerofoils opposite to each other. L caUed^ 
„ blplW — »• - " “ ^ 



Figure 1. Busemann biplane. 
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literature and the author's 


earlier report (1989). 


The drag of the busemann biplane according to the two dimensional 
linear theory as a function of the Mach number Is given In Figure 2. 



Figure 2. Drag of the Busemann biplane as a function of the Mach 
Mach number obtained by two dimensional linear theory. 


The Busemann biplane Is referred to as the -Buseaann Scranjef 
(Supersonic Coabustlon Ramjet) when fuel Is Injected Into the system. 
The combustion of the fuel at supersonic velocities and at low 
densities Is very difficult to achieve. To obtain good mixing, and 
hence efficient burning of the fuel. Is a study In Its own and will not 
be discussed here. However, the allowable t.aperature within the 
scramjet should be limited to 2000 K; since, the fuel will not burn at 
higher temperatures. 
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In tb e context or ttl. study. the coition I. r.pr.s.nt.0 by hoot 
addition .t th. throat. addition or .... or change In the 

coapressiblllty r.ctor which 1. — U» — **- 


3. Jhe Numerical Code 

, wrltten u based on th. two-dlaenslonal unsteady Euler 
The code written i 

. -thod The geometry adopted to represent 
equations and the Godunov s method. 

VI ni ana Due to the symmetry 
, , tha t of the Busemann bl plane. uu 

the scramjet Is that oi 

v,, f of the duct Is considered, 
only the upper half ot tne a 



Figure 3. 


■me adopted geo^try or th. scr»jet 


The upstrea- (region 1) and downstreaa (region I) regions are added to 
the Buseaann biplane conTlguratlon. Tb. regions 2 and 3 are teraed th. 
-co.pr.ss ion- and -expansion- regions (Figure 3). The angle between 
the x-axls and the coapresslon (drag) surface Is tahen to be .dual to 
th . angle between th. x-axls and the expansion (thrust, surface, d, 

cc«.*l«««‘ i ^^° ptl *“ thratCOndiUOnS ‘ n 1 "' £C " fl2<4 £ “ r ' 1 ' 

Der lection angles or f. 2«. »*. »* “» “ 

. u ..... the design condition was taXen to be at the upstrea. Mach 

nuaber o f e Th. results ror deration angle or 2* haw. been 
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presented In the earlier report (1989). 


4. Results 

The results obtained using the nuKrlcl cod. ere presented above In 
Figure «. these results are In close agreeaent vlth the linearized 
theory in the prediction of the drag within the syste. (Figure 2). The 
wav. drag varies with the Mach nu.be r as predicted fro. the linear 

theory. 

The code falls to give zero drag at the design Mach number. the 
non-zero result Is .or. accurate since In the region of Interaction 
between the expansion fan and the surface, the pressure drops 
isentroplcally. reducing the Mount of thrust which would have been 
otherwise developed The percentage of the expansion surface which 
intercepts the expansion fan Increases In a non-llnear fashion as shown 

in Figure 5. 

Surface contouring nay be employed to counteract the effect of the 
expansion fan on the surface pressures. This procedure ^stores the 
pressure values to the value prior to the expansion fan; however, this 
procedure reduces the projected thrust surface, hence the available 

thrust (Figure 6). 


For six degree deflection angle the ratio of the calculated biplane 
drag coefficient to the nonoplane drag coefficient at the design 
condition increases from 0.18613 to 0.29896 due to decrease in the 
projected thrust surface when the surface contouring is employed. 
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coefficient as a function of 
Design Mach number Is 6. 


12 




Figure 5 The percentage of the surface which Intercepts the 

expansion fan as a function of the deflection angle. 

If the pressures in the affected portion of the surface are restored to 

+t -tkx.4 | . t s the surface without change 

the value prior^the expansion fan Intercepts tne sun a . 

of the surface geometry^ the drag reduces to zero. 



expansion 

region 

axis of syaaetry 


The effect of surface contouring on the projected 
thrust surface. (A, A*ProJected thrust area) 
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The pressure . contours snd the pressure distribution .long the surface 
for verlous deflection angles at the design Hech nunber 6. ere given In 

Figures 7 and 8. 




the system and interacts further with a. skeclc *+- «f 

<Aoc+. The expansion fan generated at the shoulder aeets the axis 
of symmetry at AB. The reflected expansion fan therefore starts at 
point A and interacts with the original expansion fan. The effect of 
this is seen as the curved contour lines in Figure 7. 



Figure 8. Pressure distribution along the surface for 

various deflection angles. 
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A ,l.ll.r problem exist, et Hach number, mt which the reflected .hock 
weve interact, with the expansion wave. At these Hach numbers, the 
reflected shock wav. curve, and meet, the .urf.ce earlier than 1. 
predicted fro. the linear theory (Figure <|>. thl. result. In Increased 
thrust In the expansion region. Instead of the ratio of the biplane 
drag to that of monoplane drag Increasing to the value of 1.0 along a 
straight line, a more gradual Increase Is obtained as shown In 

Figure 9. 



Pressure distribution along the surface. 



Figure 9. 


( 5 ) Pressure contours within the systea. 

The Interaction of the reflected shock wave 
with the expansion wave froa the shoulder. 


4.1 Shock Sh earing 

The smearing of the discontinuities, such as shocks, presents a problem 
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The 


in numerical methods when high accuracies are required, 

calculation , of the wave drag and complex geometries usually require 
such discontinuities to be represented quite precisely. Godunov's 
first order method tends to be diffusive at weaker shocks in two 
dimensional flows. However, it does capture the shock well, and the 
calculated pressure, density and temperatures agree with the values 
obtained using the oblique shock equations within 2X in the range of 
deflection angles studied (1° - 15°). 

The Busemann scramjet used in this study is a long and very narrow 

biplane. In order to show pressure contours clearly they are plotted 

on a geometry which is stretched in the vertical direction by a factor 
of 2. This tends to spread the shock waves twice as much in this 

direction hence falsely accentuates the smearing. 


The width of the shock waves seen in these pressure contours is not 
only due to the smearing but also due to the finite dimensions of the 
cells. Shock waves will pass through several cells on each row as 
shown In Figure 10. 



Figure 10: 


Surface cells which the incident and reflected 


shock waves pass through at the design condition. 


\ 
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Figure 11 gives the pressure distribution .long the axis of sy-.tr, 
for the various deflection angles at the design condition. The snail 
hunp. A. seen In these pressure profiles 1. the result of this 
averaging process. The shoe, wave reflect, fro. the surface just prior 
t„ the shoulder causing the -11 Increase, hence the hunp In the 
pressure distribution along the axis of the sy-etry. This effect Is 
cancelled -1th the expansion wave e-natlng fro. the shoulder. 

At the design Mach number and at Mach numbers less than 
design value (at -h.ch a sy-trlcal -,v. ^distribution obtain* ,n 
linear theory) there Is - over-expansion region starting at the 
shoulder (Figure ... The flow conditions recover and reach the correct 
values rapidly. This effect again Is due to the finite d,nenslons of 
the cells -Ploy* to represent the flow and the first order nature of 
the code, and prov* to be very difficult to overcon. without post 

processing of the results. 

The advantage of the G*unoVs neth* Is Its accurate pr*lctlon of the 
fluid properties In conplex flows. Fro. the ratio of the pr*lct* 
pressure to the upstrea. pressure and upstrea. Mach neber the shoe* 
angles can be calcul.t* accurately. This allows the use of post 
processors to ell.lnat. the nu.erlc.1 diffusion of the discontinuities 
fro. the results. However, writing a post processor which deals with 
the regions -her. the Interaction of waves ta*. place. 1. difficult due 
to the complexity of the logic Involved. 
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Figure 11: 


The pressure distribution along the axis of 
symmetry at the design mach number, 6. 


4.2 A Primitive Combustion Model 

A primitive attempt has been made to study the fuel combustion In the 
Busemann scramjet by adding heat at the throat (Figure 12). This Is an 
Idealized situation, In which It Is assumed that the fuel Is added and 
completely burned in the rectangular combustion chamber located between 
the combustion and the expansion regions. In this approximation, the 
increase In mass and the change In the compressibility factor which Is 
normally associated with combustion, are ignored. 
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Figure 12. 


Heat addition at the throat 


of the Busemann scram Jet. 


„ .... — - “ — - ““ 

varies linearly with the amount of heat added. The 

, that the fuel at the hypersonic speeds does 

experimental work shows that, the lue 

J Ic reaulred to examine various 
not burn efficiently Further study Is required 

d numerically determine the most efficient ones, 
mixing functions and numerically 

W K c to be carried out on how to achieve 
Once this Is achieved, work ha 

mixing In the physical sense. 



Figure 12. 


Variation of thrust as a function of added heat 
at the design Mach number. 6. 


N 


no 



Figure 13: 


Variation of thrust as a function of added 
heat at H = 8. 



Figure 14: 


The gap to chord ratio as a function of the 
deflection angle. 


7.0 Conclusion 

The results obtained with the code, agree closely with the results 
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derived fro. the two-dl.en.lon.1 linear theory. In the linear theory 
the width of the expansion fan. hence the lsentroplc pressu 
the surface, as -ell as the non-linear effects of expansion and shock 
wive interactions are Ignored. The nt-erlca. analysis Includes these, 
therefor, the numerical results at the Mach hunters, where the 
symmetrical wave distributions and .onoplane drag obtained In the 
linear theory, are better representations of the physical solutions. 

*1 the other Mach numbers. the result, of the cod. and the linear 

theory agree within 0. 154. 

The non-linear effects increase rapidly at first with the Increasing 
deflection angle. This Increase becomes .ore gradual for the larger 

values. 

ror deciding the value of the deflection angle to be used, one has to 
consider the overall problem. As the deflection angle Increases the 
gap to chord ratio (Figure .4) Increases. The scramjet has to be 
located at the botto. of the spacecraft for aerodynanlc considerations. 

This allows the oor. dense and high pressure air behind the shock wave 

r xv- rraft to be used as the upstream 
originating at the leading edge of the craft 

air entering the sera. Jet. Hence, the gap to chord ratio should be 

such that only the air behind the shock can enter the scranjet. The 

other crucial consideration Is that, the te.peratur«s should be U.lted 

to 2000° in order to facilitate combustion. For this reason. 

i r a® - 6° are thought to be suitable in this 

deflection angles of 4 - 6 are mou* 

appl icat Ion. 

The co.bu.tlon model used In this study Is not an adequate one by any 
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The aim of this study Is to determine. If possible, a mixing 
function which will allow the fuel to nix at the hypersonic speeds. 
Once this is achieved, physical means of achieving such nixing have to 

be sought. 

It is possible to obtain more accurate results by using a second order 
method. However, this procedure is expensive in computer memory since 
it requires four times the memory of the first order method and may 
become prohibitive if the real gas effects are to be Included without 
the use of supercomputers. Supercomputers seem to be outside the 
budget of many universities at present. This may force more 
experienced researches into developing post processors which require 
better understanding of the physical problem under study, but does not 
require the use of the larger and more expensive computers. 
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Parametric Study on Thrust Production in the Two Dimensional Scram jet 


by 

Gary A. Allen, Jr. 


1. Introduction 

There are many possible scramjet geometries with different thrust production 
processes. In this paper we will examine one of the simplest geometries 
which Is the two dimensional open duct scramjet. The particular thrust 
production process which we will study occurs where the burnt and quenched 
fuel-air mixture has passed over the diverging section of the scramjet. A 
phenomena of particular Interest in this paper Is the Influence upon thrust 
by the variation of diverging section angle. 

Our theoretical model for the thrust production process was the steady, 
compressible, two dimensional Euler equations, assuming a perfect gas. The 
Euler equations were solved through a computer program using a numerically 
integrated form of the nonhomen tropic method of characteristics. The 
advantages in using this simple theoretical model were in the resultant 
computer program being small, extremely fast and usable on a low cost 
personal computer. 

2. Theory 

The use of characteristics in studying the scramjet is by no means original. 
Antonio Ferri in his classic paper [1] on axisymmetric scramjets used 
characteristics with a simple combustion model. In a paper by R.J. Stalker, 
et al. [2] the thrust production process for a two dimensional scramjet was 
described using diverging and converging Mach waves based on an approximate 
theory developed by Weinbaum [3]. Our approach to this problem employs the 
physical model used by Stalker, et al. with a nonisentropic (or more 
precisely a nonhomentropic ) characteristics method as shown in Vincetti and 
Kruger [4]. 

The fundamental equations are the steady, two dimensional, compressible Euler 
equations cast in Cartesian coordinates: 
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If we were to assume the flow was completely lsentroplc (a false assumption 
for the scramjet problem), then we could further simplify equations (4) & (5) 
Into the classical Method of Characteristics solution form: 


4 (W " 9) 
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Where u Is the Prandtl-Meyer function which Is defined as: 
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It should be mentioned there is a nonhomentroplc formulation which uses the 
Prandtl-Meyer function as a dependent variable. This form can be found In 
Llepmann and Roshko 15]. Using the Prandtl-Meyer function for the 
nonhomentroplc case Initially seemed advantageous since this formulation Is 
quite compact and accurate for the lsentroplc regions of the flow. 
Unfortunately the use of the Prandtl-Meyer function in nonhomentroplc flow 
requires carrying the temperature as a dependent variable In the vorticity 
terms. This created a host of added and unnecessary problems. Therefore It 
was concluded the pressure, theta technique shown In equations (4) & (5) was 
the superior approach. 


The classical Method of Characteristics method shown in equation (6) & (7) 
can be integrated exactly. However the nonhomentroplc formulation in 
equations (4) & (5) needs to be integrated numerically. If a first order 
finite difference technique Is used on equations (4) & (5) then the following 


Is derived: 
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The solution procedure assumes all conditions are known at points: 1 and j. 
The values for P^ and 0^ are found directly from equations (9) & (10). 
However In order to propagate the solution the Mach number for point k needs 
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3. Computer Program ^ ^ to be usable on a 

The computer program develope MS-DOS or UNIX operating 

personal computer or work station e “ P ^J^ wrltten in C) can compile first 
systems. The computer programs source ^ i>e . IBH-PC, HP-9000, 

time without error on a wid^ ^ ^ ^ for both speed 

Perkin-Elmer , Iris, etc. trm-PC/AT) The basic 

and minimum memory retirement (25 OU seconds on “ Flgure 2; 

scram iet geometry described by the computer program is depicted 
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The computation starts at the leading Mach wave which is the first Mach wave 
of the expansion fan radiating from the diverging section corner as shown in 
Figure 2. The fuel/air mixing layer along the leading Mach wave is divided 
into forty characteristic pairs. The air region above the fuel/air layer 
region is divided into another forty pairs up to the top of the duct. The 
scram jet is assumed to be an open duct type, viz. no upper surface. 

The expansion fan is divided into five regions where each region is 
subdivided into seven sectors. Six hundred and forty characteristic pairs 
are generated for each region. When the next region is calculated, all but 
the last sector of characteristic pairs from the previous region are thrown 
away resulting in the computer's data memory requirement being kept below 64 
kilobytes. This 64 kilobyte memory restriction is a consequence of using 
personal computers with the Intel 80x86 microprocessor. There are programing 
methods for getting around this problem, i.e. dynamic memory allocation with 
pointer arrays. However these alternative methods increase the program s 
complexity, and reduce portability. 

Once the expansion fan has been calculated, the program then calculates the 
characteristics down to the surface of the diverging section (thrust 
surface). The characteristics are each evaluated to determine if they cross 
prior to encountering the thrust surface. If the characteristics do cross, 
then it is assumed a shock wave has formed resulting in the computer program 
generating a warning message. Characteristics reflecting away from the 
thrust surface are not checked to see if they cross since it is not believed 
that reflected induced shock waves are significant to thrust production. 
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The characteristics leading to the thrust surface are divided Into ten 
regions each with six hundred and forty characteristic pairs. Calculation of 
total thrust is found by Integrating the local wall pressure over the thrust 
region using the trapezoidal rule. 


The scram jet Is assumed to always have central fuel Injection. Wall 
injection Is certainly worth Investigation. However to theoretically 
Investigate wall Injection, one must consider boundary layer effects and wall 
heat transfer, both of which would greatly complicate the theory and computer 

program. 

The fuel/air mixing layer as It encounters the leading Mach wave is described 
by Gaussian distributions for the Mach number. The initial velocity vector 
is assumed to be parallel to the wall. The Initial static pressure is 
assumed to be uniform across the duct. The Gaussian distributions Mach 
number are matched to the duct air layers on both boundaries of the fuel/alr 
mixing layer. The ratio of specific heats Is assumed to be uniform 
throughout the scramjet. It should be emphasized that these initial 
conditions are Input parameters based on experimental data and are not 

actually calculated. 

The program was tested against a nozzle designing program written by P. 
Jacobs 16] to establish it's basic validity. The nozzle designing program 
was developed to solve for homentroplc flow with Riemann invariant 
characteristics. It was found the two programs generated results with very 
good agreement even though the characteristic mesh resolution and solution 
procedures were different. 

4. Experimental Data 

The experimental data used for comparison with this computer program was 
taken from the R.J. Stalker, et al. paper. The data was generated In the T-3 
shock tunnel at the Australian National University, Canberra. Table 1 shows 
the test conditions for the three experiments used In calibrating the 
software. This data roughly corresponds to a flight Mach number of about 

eight. 
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Table 1. Experimental Test 

Conditions 

(from ref. 

(2] ) 

Mack Ko 

8 
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Test Case No. 

1 

2 

3 

Stagnation Enthalpy (MJ/kg)“ 

4.3 

6. 1 

8.7 

Precombustion Pressure (kPa) b 

140 

180 

160 

Precombustion Temperature (K) 

1300 

1900 

2500 

Precombustion Velocity (km/sec) 

2.42 

2.85 

3.26 

Precombustion Mach Number 

3.50 

3.35 

3.40 

Ratio of Specific Heats, y 

1.35 

1.32 

1.32 


^Experimental error ± 3%, Experimental error ± 5% 


ORIGJMAL F/’qp w 
» POOR QUALITY 



/ - distance from jet centreline 


Figure 3. Mach Number Profiles Prior to jthe Expansion Fan (from ref. [2]) 


The initial Mach number profiles in the duct are shown in Figure 3. These 
Mach number profiles are measured normal to the duct’s axis of symmetry in 
the fuel/air mixing layer Just prior to where the duct's flow encounters the 
leading Mach wave (see Figure 2). 


Figure 4 shows the pressure on the thrust surface corresponding to the three 
test conditions. The vertical axis in Figure 4 uses the pressure parameter 
AP/P. This pressure parameter is based on static pressures and is defined as: 


AP _ Pressure (Fuel On) - Pressure (Fuel Off) equ(10) 

~p * Pressure ( Precombus t i on ) 

The precombustion pressure was measured in the duct Just prior to the leading 
Mach wave of the expansion fan. The Fuel On and Fuel Off pressures were 
measured from a set of pressure taps on the thrust surface. The Fuel On 
value was the wall pressure measured when combustion had occurred. The Fuel 
Off value was the wall pressure measured when combustion had not occurred due 
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to the fuel (hydrogen) not being Injected. The horizontal axis In Figure 4 
to the luei vny « # diverging section 
is th. radial distance In .llllneter. downstream fro. 

romer as shown In Figure 2. 



x - distance from corner (mm) 

Figure 4. Open Duct Pressure Distributions (fro. ref. U) > 

. , Fieure 4 are double humped. The right most 

The pressure profiles shown in Figure pertains 

h „p is relevant to th. analysis of this paper -hll. the 
to other combustion processes within the duct. 

^ actual dimensions of th. experimental scra.J.1 — not provided In th. 
original paper. They are now provided In Figure 5. 



Figure 5. Physical Dimensions of the Experimental Scramjet 


* The Effect of Ramp Angle on Thrust Production 

of the T-3 experimental diverging section was fifteen degrees 
” 5 , A parametric study -as performed which Involved varyl^ the angle 

F 8 , rt i on The results of this studying are depicted in 

of th. diverging section. cases „„, lned cxperl.ent.Hy. 

Ft cures 6. 1-6. 3 which correspond to the three cas 
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The vertical and horizontal axis are the saae a. used in Figure 

.hows nine different pressure trace, for different diverging section ang 

starting -1th an angle of three degrees and going through a ran*, of angle. 

. to a final diverging section angle of nineteen 
at two degree Increments to a final ai e g ng 

degrees. It -as observed that the eaxl.ua pressure occurred In gu . 
a reap angle of 10.9 degrees. In Figure 6.2 the «t» Pressure occurred at 
a ranp angle of 9.7 degrees. In Figure 6.3 the -axl.u. pressure occurred at 

a ramp angle of 10. 6 degrees. 

Thrust -as calculate by Integrating th. pressure over the entire length of 
th. trace. Plots sho-lng thrust as a function of diverging section angle 
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Before this study was performed, it was speculated that the Euler equations 
might reveal an optimum diverging section angle which yielded maximum thrust. 
Figure 7 shows the total thrust increased monotonically as diverging section 
angle increased. Axial thrust is the component of total thrust which would 
overcome drag in a hypersonic vehicle. Figure 8 shows axial thrust 
increased as diverging section angle increased. Obviously the Euler equation 
model for thrust production would breakdown where the diverging section angle 
increased sufficiently to cause flow separation from the wall. Also our 
model of a scram jet assumed the thrust surface was infinitely long. Clearly 
in an actual scram jet the thrust surface must have a finite length. The 
relevance of thrust surface length was demonstrated in Figures 6.1 - 6. 3, 
where the pressure hump increased in length with an increase in diverging 
section angle. At some angle the pressure hump would spill off of the thrust 
surface. Although it was not revealed by our approach with the Euler 
equations, for a real scramjet there would be an optimal diverging section 
angle for maximum thrust production. 

6. Conclusions 

In this paper we have investigated the thrust production process in the two 
dimensional scramjet using a high speed computer program written for small 
personal computers. We have shown that for specific diverging section angles 
the peak pressure on the thrust surface achieves a maximum value. However we 
have also shown the axial thrust does not have a maximum value for a specific 
diverging section angle provided the thrust surface is infinitely long. 
Rather the axial thrust Increases monotonically with diverging section angle. 

Computer investigations of the scramjet are usually restricted to the 
province of super computers doing runs on the order of CPU hours. This study 
has demonstrated that useful information can be obtained from a small and 
fast computer program using the Method of Characteristics. 

References 

1. Ferri, A., Journal of the Royal Aero. Society 68:575-597 (1964). 

2. Stalker, R. J. , Morgan, R. G. , and Netterfield, M. P. , Combust . Flame 71: 
63-77 (1988). 

3. Weinbaum, S. , A. I. A. A. Journal 4:217-226 (1966). 

4. Vincenti, W.G. , and Kruger, C. . Introduction to Physical Gas Dynamics , 
John Wiley & Sons, New York, 1965. 



5. Llepmann, H. W. , and Roshko, A.. Elements of Gasdynamics, John Wiley ft 
Sons, New York, 1967, p. 295. 

6 Jacobs, P. A. , An Interactive Graphics Program for Computer- Assisted 
Calculation of Isentropic Supersonic Flows, Report 7/88 Department 
of Mechanical Engineering, University of Queensland, 1988. 


I 


126 


The Design of a Mass Spectrometer for use in Hypersonic Impulse Facilities 

K. Skinner 

The value of Mass Spectroscopy as an analytical tool to Chemists has 
long been established. The use of Shock Tunnels to produce high density 

particle flows (of the order of 0.1 atmospheres) with a stream energy of up 
to 10eV, although not as widely applied, has also become well established. 
Little work has been done utilising both techniques simultaneously and what 
has been done has been mostly concerned with the study of transient chemical 
species produced in the reflected shock region of the tunnel. Several 
researchers have expanded the test gas through nozzles and skimmers before 
sampling the chemistry of the gas. Crane and Stalker were the first to 

attempt Mass Spectroscopy in a high density hypersonic flow expanded from a 
reflected shock region through a contoured nozzle, whose purpose was to 

produce a uniform test flow, rather than to chemically freeze the expanding 

gas in its stagnation state. Their work, however, was never applied to 
hypersonic combustion. This work has been to investigate the use of a mass 
spectrometer to analyse the reacting gases in hypersonic combustion. 

While hypersonic combustion is uniquely suited to mass spectroscopic 
analysis, a rather severe set of constraints exist at present which act to 
limit and define the design and the performance of such an experiment. A 
mass spectrometer is uniquely suited analysis of hypersonic combustion 
because it provides the possibility of sampling free stream flow without 
causing the chemistry of the sample to be altered. Static combustion mass 
spectroscopy suffers from the handicap that the sample that is analysed is 
drawn f rom a region adjacent to the combustion chamber wall and the presence 
of the wall influences the chemistry. If the combustion is occuring in a 
flow at subsonic speeds then the presence of the instrument in the flow will 
affect the flow upstream of the instrument and hence perturb the chemistry. 
In supersonic flow it is possible to introduce a hollow conical skimmer which 
samples the flow at the tip of the cone without affecting the flow upstream. 
The cone angle on such a skimmer must be kept sufficiently low that the shock 
wave off the nose must remain attached. The cone angle of the hollow 
interior of the skimmer must be large enough to prevent those molecules which 
strike the inside walls of the skimmer from having a large effect on the 
chemistry of the axial flow. Only in hypersonic flow can these mutual 
conditions be satisfied. 

Due to the nature of hypersonic facilities which are currently 

127 


■ 



available, the design of a mass spirometer is limited In several critical 
aspects: 

j Hypersonic test facilities typically have very short test times. 

2 They deliver an impulsive load to test instruments. 

3. in all Mass Spectrometers a high vacuum must be maintained 
throughout the test. 

4 . For safety reasons combustion experiments must take place 
within a pressure vessel. 

1 . Test Time 

The -cry short test time immediately restricts the general type of Mass 
Spectrometer which may be used productively to a Time-of-Flight instrument. 
Magnetic Sector or Quadrupol. mass filters may be used, but If so, then only 
one mass ma, be sampled during each test time. A Time-of-Flight mass 
spectrometer, however, is capable of mass analysing a complete spectrum m 
on,, ,00 microseconds from a single pulse of ionised sample gas. Magnetic 
sector mass filtering or magnetic focussing are also undesirable on the 
grounds that they require precise mechanic* alignment and suffer from 

vibration problems. 

2. Impulsive Test Flow 

Regardless of the available volume of the test section of the flow 
facility, the cross-sectional area of the Instrument must be kept to 
minimum. This is to avoid the problems of choking and/or upstream influence 
which may occur when the instrument is moved close to a wall, strut or ue 
injector, and to reduce the amount of impulsive force which is applied during 
a test run. To attain any required Internal volume, therefore, the length o 
the instrument must be adjusted while the cross section is kept to the 
mimimum required to maintain sufficient expansion of the sample within « 

instrument. 

3. (a) Vacuum - Pump requirements 

The vacuum needed within the Mass Spectrometer is determined b, the 
length through which the ions must travel to the ion detector after leaving 
the ionisation region. The mean free path of particles must be longer than 
this distance, which w. ma, assume to be of the order of one meter, -nils 
means that the vacuum necessary is of the order of 0.001 Pascals. Pumps that 
are used to produce this order of vacuum are usually complex and sometimes 
delicate. Whatever type of pump is used, it should be built into the 
instrument so that the distance between pump ami the high vacuum regions .s 
kept to a minimum. Placing the pump at some distance from the instrument, 
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possibly even outside the test section vessel, limits the effective pumping 
speed to the diffusion rate along the connecting tube. 

3.(b) Vacuum - Skimmer Constraints 

As mentioned, a sample of gas from the combustion region may be expanded 
within a hollow conical skimmer with a small orifice at its tip. The stream 
tube, passing through this orifice, experiences close to a free jet expansion 
when the interior walls move away from the molecules faster than their 
thermal velocity. An important consideration in the design of these skimmers 
is the detrimental effect that the front orifice is going to have on the 

vacuum inside the instrument. The final vacuum attainable inside the 
instrument prior to a test, is determined when the leak rate through the 

front orifice is equal to the pump throughput. Small orificed hollow 
skimmers are hard to manufacture in practice, so a better ultimate vacuum may 
be attained by increasing the number of skimmers and introducing intermediate 
volumes which are kept at intermediate pressures. Each extra skimmer , 

however, introduces new possibilities for interference with the sample and 

the need for another vacuum pump. 

3.(c) Vacuum - During Test 

During a test, while gas is flowing into the instrument through the 

skimmer orifices, the vacuum inside the instrument will deteriorate unless 

the volume behind the skimmmers is large. In a time-of -flight instrument, 
the important regions requiring a high vacuum are in the ionisation region, 

the drift region for the ions, and the ion detection region. (See Figure 1. ) 
Since the molecules passing through the skimmer have highest mobility in the 
direction of the skimmer axis, the vacuum in these important regions are 

maintained longest if they are situated at an angle (optimally 90 degrees) to 
the direction of the molecular beam. Another consideration limits the extent 
to which this can be achieved, however. Ions which are formed from the 

molecular beam behind the skimmer have an initial velocity component in the 

direction of the neutral beam. Different masses, however, have different 

energies which correspond to an equivalent velocity. This means that as the 
ions are accelerated radially by an electrostatic field, and thus given 

different radial velocities, their resultant velocity vectors will make 
different angles to the axis of the skimmers. This fanning out of the ions 

can be minimised by reducing the angle at which the ions are accelerated and 
by increasing the accelerating voltage on the ions. A limit is placed on the 
accelerating voltage, however, as a high ion energy implies a faster ion 
transit time and hence a faster output spectrum. The speed of the available 



the value of the ion accelerating 


data recorder places an upper limit on 
voltage. 

4. Pressure Vessel 

Conducting experiments with combustion of flammable gases requires 
design which will »ot endanger lives and if possible, equipment. The test 
section must therefore also be a pressure vessel, and the experiment shoul 
be completely confined within this vessel. In the normal course of op.rat.on 
or an impulse test facility, this test section will hold any pressure from 
vacuum up to the maximum rating of the vessel. This means that Pumps must be 
able to operate within a fluctuating pressure environment and all electrical 

connections must be insulated against arcing. 

dearly the shape, size, skimmer arrangement, and the vacuum system for 
the mass spectrometer are strictly limited by the environment of the Impulse 
Hypersonic facility. Remaining design questions focus on the production. 

treatment and detection of the ions within the instrument. 

When deciding what method of ionisation to use. three things need to be 
considered: Will all the different species be ionised approximately equally 

in a predictable way? Will the ioniser interfere with the chemistry of the 
sample in an unpredictable way? And will the ioniser be cheap, dependable, 
small and robust? The first two requirements will only be satisfied by a 
beam style of ioniser, as all other methods require avalanches of electrons, 
accelerated in an external field. Of the possible beams, (ion. electron, or 
photon), only electron beams satisfy most of the third requirement. Ion 
beams are not simple to produce and will pollute the chemistry of the sample. 
Photo-ionisation occurs at wavelengths of 800 angstroms or less (far 
ultra-violet) for single photon ionisation, and even so, the ionisation 
cross-section is two or three orders or magnitude smaller than for an 

electron beam. 

The treatment of the ions after formation consists of passing them 
through an acceleration region and then focusing them onto an ion detector. 
The choice of ion detector is not constrained by this application. 
Constraint is placed on the means of data recording, however. It must be 
fast enough to maintain reasonable resolution over the period between the 
arrival of the first and last ions in the spectrum. With a sufficiently fast 
data recorder multiple spectra may be taken throughout the test time of the 
flow each spectra corresponding to a single pulse of ionising electrons. 

1. ’crane. K.C.A., and Stalker, R.J., "Mass Spectrometric analysis of 
hypersonic flows". J Phys. D: Appl. Phys.. Vol. 10, pp. 679, 1977. 
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Figure 1. The regions of critical high vacuum 
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SHOCK TUNNEL DRAG MEASUREMENT 


By S.L. Tuttle & J.M. Simmons 
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SHOCK TUNNEL DRAG MEASUREMENT - S.L. Tuttle & J.M. Simmons 


The develppment of the drag measurement technique Initially started by 
Sanderson and Simmons In 1988 was continued In 1989. The method was to 
allow a stress wave caused by the Impact of the flow on the model to 
propagate along a sting of sufficient length to enable the measurement 
of the resulting compressive strain before the return of any reflection. 
This requires the knowledge of the model’s Impulse response so that 
the Inverse force problem may be solved by means of deconvolution. 

The method was at a stage where the measured and deconvoluted drag would 
agree with Taylor-Maccoll theory to within 10*/.. The model geometry 
was such that skin friction could be entirely neglected ( the model was a 
cone of 15 degree semi-vertex angle ). 

The new model was a 5 degree semi-vertex angle cone of more than double 
the length of the previous one. This made the Inclusion of skin 
friction essential, and made the stress waves reflected within the model 
very apparent. 

Probably the single most significant difference was the signal-to-noise 
ratio. The strains measured were of the order of several microstrain 
only, the drag being a factor of between 4 and 10 less than for the 
15 degree cone. The noise on the signals obtained from the shots with the 
5 degree cone was in many cases the same size as the overall signal upon 
which It was superimposed. The noise In the case of the earlier 15 degree 
cone was on average only about 6Y. of the total signal, and could 
effectively be completely filtered out. 

The term ’noise* In the case of the 5 degree cone data refers not only to 
the random electrical noise, but also to the passage of small ’packets’ 
of reflected stress waves through the overall signal. These were most 
apparent In the dynamic calibrations where the loading was far more 
coherent than that occur ing in the shock tunnel. 


The process of filtering the data was thoroughly explored and it soon 
became apparent that no filter would remove sufficient noise, whilst 
leaving the overall signal acceptably unaltered. to enable an 
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Intelligible 


level of drag to be deduced after deconvolution. 


This then necessitated polynomial curve fitting to the experimenta 
in order to obtain a sufficiently clean signal, as noise Is greatly 
amplified In the deconvolution. This curve fitting approach gives quit, 
close agreement between the measured and predicted drag, although there 
are short comings with the method, again predominantly due to the 
excessive amount of noise which makes a very good curve fit difficult. 


Two methods are being used to predict the drag, each a variation on 
Newtonian theory. The first uses the measured pitot pressure, while the 
second replaces this with pv 2 calculated fro. T4 data. Both methods 
include a skin friction component which Is scaled to follow the T4 pitot 
pressure trace, starting from an Initial value calculated using simple 

boundary layer theory. 

At this stage the agreement varies between very good and about 10% 
difference at worst. Five shots were done covering four different 
tunnel conditions and examples of these results are attached, with the 
shock tunnel conditions included for each. For each shot there 
from wire resistance strain gauges and semiconductor strain gauges, both 
mounted at the same axial sting location. 


Another ares which needs to be addressed Is the curve fit which is again 
due to the poor signal-to-noise ratio and it Is believed that this can be 

improved. 

in conclusion, it can be seen that despite the severe difficulties 
Imposed by the excessive amount of noise on the signals, the drag .which 
may consist of a significant amount of skin friction, can be measured to 
a reasonable degree of accuracy for this stag. In the development of this 
drag measuring technique. 
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Shot # 1186 : Stagnation enthalpy : 
Stagnation pressure : 

Mach number 
Static pressure : 

Static temperature : 

Density 

Velocity : 


Shot # 1187 : Stagnation enthalpy : 
Stagnation pressure : 
Mach number 
Static pressure : 

Static temperature 
Density 
Velocity 


Shot # 1189 : Stagnation enthalpy : 
Stagnation pressure : 
Mach number 
Static pressure : 

Static temperature 
Density 
Velocity 
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18.67 Mj/kg 
48 MPa 
: 5.18 

24.31 kPa 
2516 K 

0.0303 kg/m 3 
5240 m/s 


10.60 Mj/kg 
49 MPa 
: 5.40 

22.00 kPa 
1500 K 
0.052 kg/m 3 
4050 m/s 


6.20 Mj/kg 
50 MPa 
: 5.81 

18.40 kPa 
785 K 

0.0811 kg/m 3 
3220 m/s 


I 






KEY TO DRAG PLOTS 


pV 2 Newtonian + skin friction 
pitot Newtonian + skin friction 
deconvoluted drag 
curve fit to measured drag 
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Force (N) 





FDrce (N) 


SHOT • 1187 - wire resistance strain guafes 


Force balance output 



SHOT • 1187 - semi-conductor strain guages 

Force balance output 





Force (N) 


SHOT • 1189 - wire resistance strain guages 



SHOT t 1189 - semi-conductor strain *“•*“ 



DEVELOPMENT OF A SKIN FRICTION GAUGE FOR USE IN AN IMPULSE FACILITY 

Gabrielle Kelly 

The skin friction gauge designed for use in T4, as described in the previous 
submission, was designed using piezoceramics. The ceramic PZT-7A Shear Plate 
was selected as the appropriate piezoceramic material due to its sensitivity 
In shear and theoretical zero cross axis sensitivity. A sensitivity to axial 
pressure was detected and in order to decouple the extranous pressure effect 
which distorts the true shear signal two identical sections of the ceramic 
were used. One was inverted with respect to the other in the axial plane. 
Thus the summation of the two separate signals cancels the pressure effect 
and doubles the shear stress signal. The design was tested in June 1989 in T4 
using an existing flat plate model: Since the signal generated was predicted 
to be small it was amplified as close as possible to the site of signal 
detection in order to prevent contamination of the signal by noise. As such a 
cavity had to be built into the model to house a suitable charge amplifier. 

In the first set of results we found the signals to be much larger than 
expected. This we summised was for two reasons. Firstly, despite the fact 
that heat transfer calculations had indicated the thermal cover thickness was 
sufficient to precludethe piezoceramic from "seeing" a temperature rise in 
the test time, thermal expansion of the cover material had not been taken 
into account. As such part of the signal generated was due to thermal 
expansion of the cap placing the piezoceramic under strain. This effect was 
in fact in the opposite direction on each of the ceramics so that summation 
of the signals canncelled this thermal effect to the first order. Thus the 
original design catered for this oversight.lt was however preferable to 
reduce this effect since the magnitude of the real shear. stress signal was 
swaped by the thermal signal and any real accuracy in determining shear 
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stress was dl.lnlshed. The other extraneous effeot was that resulting from 
the test gas coning Into contact with the electrical configuration of the 
charge amplifier. Since In this particular nodel It was Inposslbl. to 
completely seal off the cavity In which the charge amplifier was situated the 
effect of the gas contacting the electrical connections -as apparent. The 
signal was concluded to be an Ionization effect even though little Ionization 
was predicted at the test conditions used. In the presence of Argon the 
signal generated was as expected. In air It was swamped by the Ionization 
factor. Reduction In temperature reduced the effect but It was Impossible to 
remove It other than through the use of Argon. 

The next series of tests run, were using a » fUt plate mode, In which the 
cavity for the charge amplifier was able to be completely sealed. The results 
were vastly Improved. The new mode, also was designed with a thlch.r thermal 
cover which was capped off with bahellte disc .a material which exhibits a 
zero coefficient of thermal expansion). This series of tests brought the 
signals down to a sensible and expected level. Unfortunately the signal did 
not emulate the behaviour of the stagnation pressure. After an original rise 
the signal seemed to have a decaying sinusoidal wave superimposed on It. This 
behaviour was exhibited by both ceramic transucers - identical In size and 
nature, but opposite In sign. This would see. to suggest that the effect was 

. , a. Ko _ vial direction of the transducer and perhaps 
one experienced in the axial direction 

. thp , nresence of stress waves through the 

could be explained as the presence 

thermal cover to the ceramic. Unfortunately due to the limited time 
available In T4 the tests were concluded In October 1989. Due to the 
fact that I wish to bring these tests to a conclusion as soon as 
possible I have opted to test the gauges In TQ. This facility is 
less demand and can be operated by the user alone and has the added 
advantage of a short turn around time. This effectively means the gauge can 
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undergo more thorough testing and 1 can reach conclusions as to the last 
extraneous effect as soon as possible. I feel confident that this behaviour 
can be explained by either the presence of stress waves or a still 
unresolved thermal effect. The results of these tests will be finalised in 
the next two months. 
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HYPERVELOCITY FLOW IN AXISYMMETRIC NOZZLES 
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AKTIUCT. 

We describe a simple procedure for the design 
of axisymmetric supersonic nozzles for use In 
ref lec ted -mode shock tunnels. This method has 
been used to design a aoderate Mach number (M ■ 4) 
nozzle end a high Mach number (H • 10) nozzle. 
Both nozzles have been calibrated by measuring the 
pitot profiles near the nozzle exit planes and, 
although the H « 4 nozzle performs well, the 

M * 10 appears to have reached a Mach number - 
pressure limit in which the unsteady nozzle 
boundary layers significantly affect the test 
flow. 


INTRODUCTION. 

A new generation of aerospace planes Is 
currently being developed These vehicles will be 
powered by air-breathing engines and so spend a 
considerable fraction of their flight accelerating 
through the atmosphere This means that high 
speed (v - 3-Tkm/s ) aerodynamics will again become 
a focus for fluid dynamic researchers. 

Experimental facilities capable of providing 
test flows at these speeds supply a pulse of test 
gas that lasts only a few milliseconds. To make 
the best use of the high enthalpy test gas. It Is 
expanded to a high-speed uniform and parallel 
flow. For the shock tunnel facility T4 (Stalker A 
Morgan 1988), this is achieved by using an 
axisymmetric nozzle in which the gas is initially 
allowed to expand through a conical section and is 
then redirected by the contoured part of the 
nozzle wall to produce a uniform test flow at the 
nozzle exit plane. 

Impulse facilities such as T4 typically 
operate at total enthalpies (H^) of 10 - 30 MJ/kg 

with associated stagnation temperatures (T^) of 
5000 - 12000 K and stagnation pressures (P R ) of 

40 MPa. At these conditions there is a strong 
coupling between the chemical reactions of the 
dissociated air and the axisymmetric gas flow. 
This coupling complicates the design calculations 
significantly However axisymmetric nozzles for 
non- ref lecting shock tunnels have been designed 
using the method of characteristics (HOC) with 
chemical reactions by Hudford ef ai (1980). In 
their design they note that the contour shapes 
computed for a chemically reacting flow were 
similar to that computed for a perfect gas with 
a suitably chosen ratio of specific heats, t 

This leads us to approach the design process 
In a simple fashion and treat the total flow as 
two relatively simple flows patched together. We 


treat the early expansion of the test gas in the 
conical section as a quasi -one -dimensional flow 
with finite-rate chemical kinetics. At the exit 
of this cone we assume that the flow is a uniform 
source flow and that the chemical reactions are 
frozen. We then treat the flow in the contoured 
section as axisymmetric flow of a perfect gas 
where 7 has been chosen to approximate a 
chemically rmactlng flow over the same expansion 
The HOC calculation of the flow in the contoured 
section Is then the "standard* perfect gas 
procedure described in Liepmann 8 Roshko (1957). 

This approach has been used to design a 
nozzle with exit Mach number, M » 4 for use in 
supersonic combustion studies (see figure la) 
Its performance in terms of flow uniformity 
appears to be adequate A high Mach number 
(M ■ 10) nozzle was also constructed (see figure 
lb). It has a two-stage initial expansion to 
reduce the time required for the test flow to 
reach steady state. Unfortunately, for nozzle 
stagnation pressure P, < 40 MPa. the boundary 
layers that develop on the nozzle wall 
significantly affect the uniformity of the flow. 


MACH 4 NOZZLE. 

We will consider the design of a specific 
nozzle with M * 4 and a throat diameter of 25. 0mm. 
at nominal stagnation conditions P. * 30.4MPa. 
H, - 16. 2MJ/kg (T. » 8000k). The conceptual 

layout of the nozzle is shown in figure 2. For 
computational convenience, we set the origin of 
the axial coordinate, x. at the start of the 
conical expansion. 

For the Initial conical expansion, we choose 
a cone half-angle of 12° and we assume the flow to 
be a uniform source flow at the end of this 
section (x * 98mm). Zonars (1967) indicated that 
throats with a constant area section with length 
equal to the throat diameter produced "excellent 
high temperature source flow characteristics when 
measured in conical nozzles". Hence we do not 
analyse the transonic flow near the throat as done 
in a number of other studies ie.g. Sivells 1978). 

The flow in the conical section is strongly 
influenced by the chemical reactions associated 
with the dissociation and recombination of the 
molecules in test gas We treat this flow as a 
quasi -one-dimensional flow with finite rate 
chemical kinetics. Computations are performed 
with a FORTRAN program WENZF (Lordl tt ai 1964). 
The stagnation region (labelled 1 In figure 2) and 
the subsonic flow up to the throat is considered 
to be in chemical equilibrium but. once the 
program steps into the conical expansion, a 
nonequilibrium chemistry model Is used 
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figure 1. 



M • 4 and M - 10 axlsymmetrlc nozzles 
constructed for the T4 shock tunnel 
facility. 


M * 4 


100m 
l i 
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M » 10 


We use the flow properties (H and jr) at the 
end of the conical section as the starting point 
for the HOC calculation of the flow In region 3 
(source flow), region 4 (transition from source 
flow to parallel flow) and region 5 (parallel and 
uniform flow). The gas flow Is considered to be 
chemically frozen and a value of M at the end of 
the conical expansion Is calculated using an 
estimate of the speed of sound 
a » (y 5 Kq 1000 T) os , where S is the sum of the 
•P^cles concentrations (mole/gm-mlxture) 
*0 » 9.314 J/gm-mole/K and T Is the static 
temperature in K. Here a - 1380 m/s. The 
computation was performed with the aid of the 
program ’N0C“ documented In Jacobs (1988). The 
Inlet boundary is specified as a uniform source 
flow with M 4 • v/a • 2.804 and z • 1 334. Only 12 
points were used on the Inlet boundary as the 
program retains all of the mesh data In the memory 
of the microcomputer We believe this to be 
sufficient because Schurmeler (1959) indicates 
that a 10 point mesh produced a 0. IX error in exit 
height for a two- dimensional nozzle calculation. 
The Mach number on the axis at point C (In figure 
2) Is computed to be H c • 4. 12. We then step 
along the axis downstream of point C and compute 
the flow field In region 4 by proceeding upstream 
along characteristics such as CA. 

Once the characteristic mesh is generated, a 
*t r «**llne can be Interpolated through the mesh, 
starting at point A and finishing where It 
intersects the characteristic CD. The data points 
on the Interpolated streamline are then used In a 
spline fitting routine to define the nozzle wall 
as a cubic spline with eight knots and specified 
end slopes See table 1 for the wall coordinates. 
A boundary layer modification was added to give a 
total wall radius 



8 C Nozzle Axis 


Figure 2. Conceptual view of a hypersonic nozzle 

1. Stagnation region 

2. Transition to source flow 

3. Source flow region. 

4. Transition to parallel flow 

5. Uniform flow parallel to nozzle 
axis. 


x * 

r total r lnviscld * 511 *exlt 
where 8 was estimated to be 1 4mm 


Calibration of tha Mach 4 Nozzle. 

The performance of the nozzle was evaluated 

by measuring the pitot pressure. P at a 

pi tot 

plane normal to the nozzle axis and located a 
distance z • 120mm downstream of the nozzle exit 
plane. Each pitot probe was fitted with a PCB-112 
plezo-electrlc pressure transducer which measured 
the stagnation pressure behind a detached shock 
that formed over the upstream face of the probe 
Several probes were mounted In a rake and a number 
of shots of the shock tunnel were required to 
build up each pitot profile 

Figure 3 shows the time history of some of 

the pressures. The "raw- traces for both P and 

s 

%ltot show the impulsive start and 

subsequent decay associated with "under-tal lored" 
°P erm ^l° n of the shock-tube In this mode, the 
shock that compresses the test gas reflects In 
such a way as to allow the compressed test gas to 
expand back up the shock-tube. This reduces P 

s 

during the test flow time but delays contamination 
of the test gas (air) by the driver gas (helium). 
To eliminate the time variation from our pitot 
pressure measurements, we normalize P . by P 

pitot 7 s 

with a suitable time delay (here 0 2 msec) and 
filter (time constant • 0.05msec). The normalized 
traces for 4 pitot probes are shown In figure 3B 
We then discard the first 0.5msec of the trace 
which contains the starting pulse and measure the 
mean value over the next 0 5msec. This mean 
value, together with an estimate of Its variation 
over the test tier is plotted as a single point on 
the pitot preset.- profile 

Figure 4 shows pitot profiles for a nominal 
stagnation pressure P # - l3Mpa and three 

enthalpies H # * 16. 8 8 and 6.6MJ/kg. There is 

so ** fill-off at the edge of the core flow but 
this Is due to the .expansion fan from the trailing 
edge of the nozzle propagating into the test flow 
The performance of the nozzle appears to be 
satisfactory but we cannot measure how parallel 
the flow is because the pitot probes are 
Insensitive to small changes In flow angle. 
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Flgwe 4. pitot pr«**ur« profile* V* * “ 4 

^ , _ t> - iKna H ■ 120ea 


A -Raw" atagnation end pitot pressure 

B Bores 11 sed pitot pressure. 


MACH 10 MOZZLX. 

The atarting process for the flow in • 
hypersonic nozzle Involve the Kop^.Uonof 
prlsary and aecondary 

•xpanslon through the no2zl * S * lth the f low to 
«h£m that the tl. ( t^. required for the nou io 

„,„=h ..... ■;"-^;; r U , " e .T«'.2 

tlae taken to sweep the upetreaa nemo 

unsteady expansion out of the nozzle. 

nozzles with fixed divergence angle, t^ Increase 

rapidly -1th M. To reduce t # (and hence reduce 

the a-ount of gat consiaed In the 
process) for a nozzle -1th large * * 

Increase the divergence angle of the lnl 1 
conical section. However. If the angle Is 
large, the flow -111 separate. 

U. ulll no- describe a nozzle with M - 10 

-v .trjsjsr 1 tsss>» 

through * cone i ¥lth * h . wms the largest 
previous ^rlsnce. th r#msonably uni for. 
divergence that wouic g» half-engle i» 

a^.*rr=sj^srs 

straighten the flow. 

.T. - >■«*> - 


nozzle. P # ■ 13Hp* 

A. 

- 6.6HJ/kg 

B. 

• g.BMJ/k* 

c. 

- 16HJ/kg 

a. 

NEKZF was uitd 

1.41 

at the etart 


of the contoured 


section. True re *»*.*▼« *7 ***•*■ 7 7 ~ r * mr 

result of the large fraction of aonatoaic specie 
ft the end of the conical 
Coordinates for the complete contour are *‘ v "7 d 
table 2. Ho boundary layer modification -»»»dded 
to these invlscid coordinates. To reduce the 

length of the fabricated nozzle without effecting 
the usable test core, -e truncated the **>«** *1 
the x-st*tlon -here the characteristic CD 
Intersects the esUsat^ posltlon £/* r £». xlt 


I i* 


the boundary layer. Because of the large exit 
Itod sabtr. the fabricated nozzle -as such 
m^te^Ti . 1036.) th. the full Invlscid 
contour (J “ 1T77.). 


Calibration of the Hsd» 10 hozsle. 

to for the H • 4 nozzle, the ** rf or **"7 
the H - 10 nozzle -as evaluated by -assuring the 
pitot pressure do-nstresn of the no**l« ** 
plane Figure 5 shows pitot profil es, for 
conditions* Profile A (P. - ZOKPs. ", * 

indicates that the nozzle la not performing -ell 
at^ .v even have a conical shock present In the 
YL even » ^ ^ flnU .. difference 

program and obtain*. ^l-Uar *. 

^r^promT. possibly because the ^ 

nsssed out of the test flow core st this axial 
position. The s*p..lon PropM*^ ^ £ 
tralllna sdge of the nozzle ha* reduced me 
dlaneter of the teat core at this axial position. 
He note also that the normalized pitot traces for 
this condition changed continuously through the 
Yl., ttas This indict., that the boundary 
££r. on the nozzle -all did not reach steady 


146 


I 







Itlt*. 


0 002 


Operating the nozzle at higher pressure (see 
profile C. • 40HPa) significantly laproves the 

flov uniformity but there Is still a depression In 
the centre of the pitot profile. We suspect that 
this improvement In nozzle performance Is a result 
of reducing the boundary layer that grows along 
the nozzle wall For profiles A and C. we 
estimate the total boundary layer thickness at the 
exit plane to be 4 • S2mm and 6 - S8m. These 

- 1/2 C 

values scale as Re * which indicates that the 
boundary layer Is laglnar. Assuming that 
displacement thickness 4 • 0. 33 4, the effective 

•rea ratio of the nozzle is reduced from the 
design value 2290 to 1450 and 1680 for profiles A 
and C respectively. Values of pitot pressure 
computed with WDCF for these effective area 
ratios are shown as dashed lines In figure 5. 

In conclusion, it appears that our design 
approach Is satisfactory for hypervelocity nozzlds 
with moderate exit Hach number. However, nozzles 
with high values of exit Hach number need to be 
operated at sufficiently high pressure to ensure 
that their boundary layers do not disturb the core 
flow. 
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Pitot pressure profiles for the M ■ 10 
nozzle. The dashed lines Indicate the 

NENZF computed level with 4 /4 » 0.33. 

A P # • 2 OH pa. - 30HJ/kg. z - 139m 

B. P % • 2 OH pa. - 30HJ/kg. z - 369am 

C. « 4 0M pa, - 25HJ/kg. z - 139mm 
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Table 1. Coordinates for the M « 4 inviscid 
contour. Points with x a 98mm are 
knots on a cubic spline. 


x (mm) 


r lnvlscld ** 


comment (slope) 


0.0 

12 5 

start of 12* cone 

98 43 

33.43 

start of contour (0.212< 

157. 48 

44.45 


216.54 

52.79 


275 59 

58 56 


334 65 

62.30 


393 70 

64.66 


452 76 

65.88 


511.81 

66.07 

end of contour (-0.0048 

Table 2. 

Coordinates for the K • 10 Inviscid 


contour. 1 

Points with x a 47ma are 


knots on a 

cubic spline. 

x (mm) 

r lnvlscld 

(am) comment (slope) 

0 0 

3.00 

start of 15* cone 

12 0 

6 24 

end of 15* cone 

15.0 

7.09 


18 0 

8.03 


21 0 

9.00 


24 0 

10.03 

start of 20* cone 

47. 70 

18 66 

start of contour (0.3640 

294 78 

74.51 


541 85 

104.66 


788 93 

121. 57 


1036 00 

132 76 


1283 08 

139 13 


1530. 16 

142 44 


1777 23 

143 55 

end of contour 
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REAL gas effects m hypervelocity FLOWS over AN INCLINE nc 

Kw 


by 

R.M.Krek 


1. INTRODUCTION 

. he „ measurements obtained from experiments on a 

Pressure and heat transier 

, 0 , 15 ° at angles of Incidence p - 0 and P - 

cone with a half angle Three 

. hypervelocity Hot. with a nomina, Mach number of S are P— ™ 

H - 6 28. and 36 MJ/kg. The pressure 

stagnation enthalpies were used. . 

results agree rm.son.bl, well with Ne«tonl» theory. Taylor- acoo 

and with detailed three-dimensional calculations when the flow » roz 

.assure measurements at high entlmlP, — — — J 

Heat transfer results are also given and are 
enthalpy on the leeward flow. Heat — ^ ^ from 

compared with theoretical values for zero inci en 

I 30- incidence exp., men, s. show the - Trends as - Pressure 
distribution around the cone. 

The aim of this project is understand the basic dissociation pheno 

which occur in three dimensional high enthalpy flows. A simplif 
consisting of an inclined con. In a h,P«rv«locit, flow of pure Nitrogen is 

a f this flow at sufficiently high enthalpy should 

used Experimental study of this flow a 

ton between nonequilibrium dissociation 
provide insight into the mterac i evaluate the 

chemistry and hyperve.ocity aerodynamics. and » 

. ™ effects on the pressure and the heat transf 

influence of real gas 

distributions on the cone. 

r hvnervelocity flow about an inclined cone is 
The expected structure of hyperveloc 

• , shock is of greatest strength at the 

i The conical bow shoot is a 

shown in figure l. 1UC . . If 

i hut weakens with increasing azimuthal ang e 4>. 
windward smgnation plane normal 

Urn angle of incidence A. is sufficiently Urge, the 
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to the cone exit is supe™onic. As the flow expends around the cone surface, 
the fluid may be expected to first accelerate, and then decelerate as It 
approaches the leeward plane of symmetry. The pressure recovery required at 
the symmetry pUne requires a prfr of cross-flow shocks. Owing to the strong 
entropy gradients along these shocks the eortlcity oT the shock-processed 
flow Is greatly enhanced, with the resultant formation of a complex 
three-dimensional flow containing one or more pairs of streamwise vortices, 
and possible shock induced separation from the cone surface. This flow 
structure Is essentially Inviscid (see Marconi 1989) but will be modified by 
the presence of shock induced boundary layer separation. We note that for 
frozen flow at low enthalpy, lie. constant r). the approximately conical flow 
will contain no dominant length scale. 


2. EXPERIMENTAL PROCEDURE 

•me experiments were carried out in The University of Queenslands T4 
free piston shock tunnel which is capable of achieving the high temperatures 
and thus the enthalpy conditions which are experienced during reentry. A 
Nitrogen test gas was used to further isolate the basic dissociation chemical 
reactions which occur at these high temperatures. The expected freestream 
conditions in the test flow for the three enthalpy conditions are presented 

in Table 1 below. 


H 0 

MJ/kg 

T 

CD 

K 

P 

CO 

kPa 

kg/m 3 

u 

CD 

km/s 

Mach 

N° 

7 

a 

00 

6.42 

775 

5.88 

.0256 

3.35 

5 .84 

1.4 

~ 0 

28. 0 1 

2560 

9.92 

.0099 

5. 67 

4.79 

1.38 

.140 

36.0 

3518 

10.05 

.0076 

6. 17 

4.49 

1.44 

.268 


Table 1 Freestream Conditions in Test Section 
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The cone used for the experiment, had . heir angle • - »* 
lOTgth l80 nun. The Instrumentation on the cone consisted orlgln.ll, ot ten 
Eh™ pressure transducers and «n Pt-PtRh thwmocouple transfer 

gauges. BOO, the pressure transducers and the re mounted ,n 

COS,, of rise gauges each. -1th the rows ^spaced around the cone. The 

arrangement of the cone and it's Instrumentation Is shown in figure 1. Due 

. „ s,. ms we W ere only able to obtain five gauges that 

to pressure transducer problems, we were « , 

operated satisfactorily. 

3. PRESSURE RESULTS 

The results from the calibration experiments at « * 0 are displayed in 
figure 2 which show plots of P,/Ppitot versus cone posit 
conditions with H.- 6.« Mt/V* and H.- g7.PS Ul/kg. ™e experiments « 

compared with T.,lor-Maccoll -.lotions, for which an appropriate ~l»e of » 

* * .... it was found that the flow is 

is required. For the equilibrium case, it was 10 

, . 1 4 For H =28 MJ/kg. the 
vibrationally frozen, and that we can take r 

, a,u. „ s 4+ot/3. assuming ideal 

flow is relaxing vibrationally and we tak 7f 

dissociating gas behaviour. Figure » show. Urn. the measure 

agrees quite -ell with the T.,lor-Maccoll solution when H„- 6 UJ/kg. The 

hv 57 from the Taylor-Maccoll theory. For 
average of the experiments differs by 57. rrom xne y 

the non-equilibrium case, figure 2b. the average value of the pressure 
signals is 28 heiow that of the Taylor-Maccoi. vaiue. Ap^t from the 

discrepancy in figure 2b. it can be seen that the pressure on the 

surface during the experiment at aero incidence, for both enthalpies is 

approximately constant along a ray. 

NOW consider the results obtained from aliments for a - 30". Figures 

3-S show plots of pressure in two forms. 1) P^ma, •“<* 21 1 “ (P v /P . ) 
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verses f at each of the five measuring stations, for H c - 6, 28 and 36 MJ/kg, 
respectively. The experiments are compared with the Newtonian estimate for 
pressure, which in the hypersonic limit and assuming perfect gas behaviour. 

is given by, 

P ,/P - >M*(sin 0 cos p ♦ cos 0 sin P cos $) 2 ♦ 1. ^ 

Figure 3 also shows results obtained from detailed three-dimensional 
calculations performed by Macrossan et al (1989), and these compare 
favourably with experiments. For all enthalpy cases it can be seen that the 
Newtonian value of pressure for the first 60° of the windward surface is 
higher than the values derived from experiments (figures 3b-5b). 

The main aim of this project is to investigate the effects of real gas 
chemistry on the leeward surface of the cone. From figures 3a-5a it can be 
seen that the data provides some evidence that cross-flow shocks are 

occurring on the leeward surface, as suggested in figure 1. The cross-flow 

shocks occur at about * - 150°. but the resolution of the plots and the 
reliability of the pressure gauges make any conclusive quantitative 

assessment of the results difficult. It is proposed that further experiments 
on the leeward surface of the cone be conducted using more sensitive pressure 
gauges, and that a finer survey of the leeward surface be taken. This should 
allow accurate measurement of the azimuthal position of cross-flow shocks. 

4. HEAT TRANSFER RESULTS 

As said above heat transfer was measured via the use of Pt-PtRh 
thermocouples. Calibration of the thermocouples independent of the 
experiments proved impractical, thus the thermocouples were calibrated by 
scaling the gauge measurements for heat transfer to the cone surface at zero 
incidence, and in low enthalpy flow, to the predications of Eber (1952). 
From Eber’s laminar boundary laminar analysis, the heat transfer at a 

151 


I 



, . alon£ the cone surface may be expressed as 

distance x from the apex along 


_ -in _-i/3 
St, ■ 0.575 Re* Pr 


( 2 ) 


. St_ Re, the Stanton number. Reynolds nun, bar. and Brandt, 

where St*. Ke,, 

number restively and me bamd on -.stream cation.- 
Knowing tbe tbeoretioa, values tor tbe beat 
c^ibrat. the experimental data. Figure 6 shows the heat transfer » 
plotted in terms of St verms Re. for the three enthalpy oondtt.ons and are 

compared with their respective theoretical values. Figure * 

. a Ar *c ran be determined for 
f „ _ W hich the calibration factors can oe 
equilibrium case from whicn in 

end individual thermocouple. Figures *b and Sc are the two noncuii.br, 
conditions, from which it cm, be seen that calibrated experiment., data me 

higher than the theoretical values. 

The dad obtained from experiments performed at B * 30 me nornrnl.se 
wlth the zero incidence data for each enthalpy condition. Figures VS show 

st/st verses ♦. for H„- 0. «. - * 

dtw to the heat transfer until there to a .oca, minimum - ~ 

, . with a nse to the heat transfer “ *** 

pressure distribution mound toe cone. Thto -gain surest toe existence of 
separation and counter-rotating vortices occurring on the leeward mnfac. 
as with toe pressure data toe resolution of heat transfer is not suffeen 
to establish an, quantitative arguments of to. leeward surface. 

5. CONCLUSIONS 

w. pressure and heat transfer results presented compare reasonably with 

their respective theories. Both toe pressure and heat transfer suggest - 

, thf. leeward surface, but more detailed 

existence of crossflow shocks on toe leeward 

* ,K,t will be conducted will have 

experiments me required. The experiments that «,U be 


152 


more sensitive pressure transducers and thinfilm heat transfer gauges to 
investigate a finer survey of the leeward surface of the cone. From these 
experiments we shall be able to pinpoint the crossflow shocks and determine 
the effects of varying enthalpy on the crossflow shocks. A further use for 
the thinfilm gauges is to compare actual experimental data with theory at 
zero incidence and to calibrate the windward thermocouples. 
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Figure 2: Peoo^W « r »« cone P 08 * 1 * 011 fcr ^ “ 0 *‘ 

(•) Ho “ 6 MJ/kg, r » 1.4. M * 5.84. ♦ *500. • *501. * *502. 

(b) Ho « 28 MJ/kg. r - 1-38. M - 4.79. «■ *503, • *511. x *512 

Taylor-Maccoll Theory. 
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•2 

•3 



Figure ±: fi « 30°. H. « 22. MJ/kg 

(a) lnfPj/P^) verses * 

(b) P c /P PH ot »«rses 4 > 

157 . Position 1, ♦ Position 2, x Position 3, • Position 4. • Position 5 


Newtonian. 




pitot 
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BUB S ' - 30*. H. - * « U/ *« 

(a) lntPe^J »er*es * 

(b) P^pitot Tcrtes ♦ 

. position I. ♦ Po*Won 2. x Position 3. 

Newtonian. 


Position 4, • Position 5 


ft 0500 before 
0 0501 Hftn coltbroUou 
4 0500 etUfcrftU* 
x x501 ®*Ubrote4 


► 0500 before 
» Til before eellbroliou 
4 0003 c*Ubroie4 
i ifill eellbrete4 


JUjDol4f Number 


Jtejuoldj Number 


0551 before coMbreUo® 
s562 before wlibntloa 
0561 eaUbr»te4 
b 562 celibrote4 
theory 


t*pot4i Number 


rimre 6: Sunton Number eeree. RejnoWe Number for t - O* 

(,) H. - 6 UJ/k*. (b) H. - 28 NJ/k*. tc) H. - 36 NJ/k* 








. on° H - 6 UJA« 
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SHOCK TUNNEL DEVELOPMENT (Supplementary Project) 


R.J. Stalker and R.G. Morgan 


Raising Operating Pressures 

A major target of the 1989 program was to raise the operating pressure levels 
of The University of Queensland shock tunnel T4. Specifically, this meant 
raising the pressure after shock reflection at the downstream end of the 
shock tube - i.e. at the nozzle entrance. With a free piston shock tunnel, 
this can be done in two ways. Either, by operating at a high driver gas 
volumetric compression ratio, a high main diaphragm burst pressure. The 
Page-Stalker effect (i.e. the loss of pressure at shock reflection in the 
driver gas at the downstream end of the shock tube), then causes the pressure 
at the nozzle entrance (nozzle stagnation pressure) to be reduced. Or, by 
operating at a lower driver gas volumetric compression ratio, and a lower 
diaphragm burst pressure. The reduction in the Page-Stalker effect then can 
allow similar nozzle stagnation pressure levels to be reached. 

In both cases, the pressure driving the piston remains the same, ensuring 
that the same energy is imparted to the driver gas during the compression 
process. The first method has the advantage that longer test times are 
obtained at a given stagnation enthalpy. The second method has the advantage 
that lower peak pressures are produced in the apparatus. 

Tests were made to compare the operation of T4 in the two modes, and 
resulting nozzle stagnation pressure records are shown in Fig. 1. Record A 
was obtained with a main diaphragm burst pressure of 170 MPa, a driver 
volumetric compression ratio of 120 and a piston driver reservoir pressure of 
8.3 MPa. Record B was obtained with a main diaphragm burst pressure of 
86 MPa, a driver volumetric compression ratio of 60, and a piston driver 
reservoir pressure of 8.9 MPa. 

It can be seen that the nozzle stagnation pressure for B is a little higher 
than for A, even though the main diaphragm burst pressure is only one half of 
that for A. Thus a given nozzle stagnation pressure is produced with much 
lower peak pressures in the apparatus. 

The number of flow passes over an 0. 5 m long model also is shown in the 
figure, together with the estimated time at which contamination by driver gas 
occurs. The stagnation enthalpy for A was 27 MJ/kg, whilst that for B was 
13 MJ/kg (the test gas for A was nitrogen, and that for B was air, but this 
does not influence the operating characteristics under discussion) and it can 
be seen that approximately the same number of flow passes takes place before 
contamination in each case, in spite of the fact that increasing stagnation 
enthalpy with given driver conditions reduces the time for contamination. 
Thus, changing the driver conditions by raising the diaphragm burst pressure 
and the driver volumetric compression allows operation at high enthalpies 
without reducing the effective test time. 

It might be noted that the higher diaphragm burst pressures are associated 
with considerable technical problems. For example, the high pressure 
stainless steel diaphragms which are used do not "petal" cleanly, and 
fragments of the order of 10 to 20 mm in diameter tend break away, to impact 
the downstream end of the shock tube, or to pass into the nozzle. The steel 
high pressure plate of the entrance to the shock tube tends to ablate, and 
the high stresses induced in the nozzle by high shock reflection pressures 
cause deformation of the surface contour. Although all of these are problems 
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can be solved by the application of time and technical effort, they are 
nevertheless a practical discouragement to high pressure operation. 

Hypersonic Combustion 


The purpose of raising the tunnel operating pressure levels was to make It 
posslbl^ to produce reasonably vigorous, hypersonic combustion of hydrogen 
f^el With the pressures which were achieved, it was considered that It may 
be necessary to Increase the combustion duct length in order to allow an 
adequate streamwlse distance for the combustion to develop. Accordingly, the 
long 11 duct designed for the scaling experiments see Press ^*:^^ 
Correlations In Supersonic Combustion") was used. The Intake was modified to 
allow effective "direct connect" operation, rather than the expansion after 
tie intake which is shown In Fig. 2 of the above report, and the same central 

Injection strut was used. 

Static pressure distributions along the duct are shown In Fig. 2 Tests 
conducted not only with hydrogen, but also with ethane, as fuel. Although 
the wave structure resulting from the modified Inlet somewhat confuses the 
display of results, the clear cut Increase In static pressure, amountin f *° * 
factor of two , which is achieved when air replaces nitrogen as test gas Is 
evidence that combustion Is taking place. It Is Interesting to ^ note that 
ethane produces combustion pressure Increases which are only a little belo 
those experienced with hydrogen. 
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ST/>6. MtESSUhL US Tlttf. 



H. - NO. or FLOW MASSES (C> ,S 
c - 3>e«vrc-« Gas Cont^mimaTioN 


STAG. P1ESSUU US TIHE. 



FIG. T DRIVFR OPTIONS 
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PRESSURE VS POSITION, 
absolute tiie = 8.988 isc 
flow velocity = 3708.888 «d/»sc 

5km HiU steel rupture? 71 Up# He lira driver coipression ratio 5S 

Nowinal enthalpy 8.4 Hj/Rg 
* HZ/air phi 1.3, ru* 1858 
filled square C2H6/air pki 1.2, ruu 1874 
hollow square HZ/NZ equivaleut phi 1.3, run 1859 



( M — 5*5 , Stag. Press. - 67-5 MPq. ) 
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INVESTIGATION OF FLOW CHARACTERISTICS JN 
TQ EXPANSION TUBE 


NASA Report 1990 


Andrew Nedy 
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In the eariy part of 1989 the work on the analytical predictions of test gas 
flow conditions for an ideal dissociating gas and an ionizing gas m an expansion 
tube was continued the work continued to the point where the simulation was 
shown to worit for low enthalpy flows but was still having problems with the high 

enthalpy flows. 

The program Is driven by the Input of the initial fill pressures in the shock 
tube and the acceleration tube and the measured shock speeds in the shock and 
acceleration tubes. The calculation process is set out in flow chart form in 
FIGURES 1 Si 2 . The governing equations for the flow processes were set out in 

last year's report. 

Initial calculations with the program using the expected end velocity 
from 1 - D perfect gas calculations returned a negative value of a final indicating 
that fun recombination had occurred during the expansion process. Substitution of 
a = 0 into the velocity expression confirms this. That is it confirms that ionisation 
will fall to approximately 0 before the expansion process is complete and the end 
conditions are reached. The point where this occurs is designated the point of 

transition. 

It is necessary to find this transition point, but as the exact isentropic 
relation fails for a zero level of ionisation it is required to find the finite but very 
small level of ionisation which when reached during the expansion can nolonger 

contribute to the flow energy. 
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This transition level is found by varying a final and observing its effect on 
the final Test Gas pressure Pj\ This pressure can be found by relating the gas flow 
conditions through the unsteady expansion from the secondary diaphragm but 
also by relating conditions back across the secondary shock. This is possible because 
by physical constraints the pressures on either side of the acceleration gas/test gas 
interface must be equal. The test gas velocity is similarly constrained to the velocity 
of the acceleration gas behind the secondary shock. 


The relations used to find P4 and U4 are 
2 y M 2 Y-M 


P 4 = P3 


Y+l Y+T 


u < * $ * 


where 


U4 

As 


M 4 = 



y-\)\ 
2 >) 


T 4 

Ts 


YMs : 


Y-lV, Y - 1 t)* 


1 + 


2 M 3 ■ 


4M3‘ 


The expansion process is thus divided into two stages. A real gas model is 
used to relate from the beginning of the acceleration tube through the unsteady 
expansion to the transition point at which the level of ionisation has reached a 
final negligible level, then a simple perfect gas model is used to relate through to 
the end of the tube Since the gas at this stage can be considered to be behaving in 
an ideal manner, using the shock speed and the quiescent acceleration gas pressure 
the conditions behind the secondary shock can be found using 1-D shock relations. 
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TRANSITION POINT 


SECONDARY 



Real Gas Perfect Gas 

behaviour behaviour 


These relations give values of velocity and pressure In the region behind 
the secondary shock. i.e U 4 and Pa By physical constraint these are known to be 
equal to the velocity and pressure of the test gas that lies behind the secondary 

contact surface. 

In the case of Dissociating gas the same technique that was used in the 
ionisation case is followed. This is because once again it is observed that when the 
level of dissociation is set to zero the resulting velocity does not correspond with 
the known test conditions, it is too low. Tbe flow process is divided into real and 
ideal gas stages, with the change occuning at the transition point. 
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The equations are set up and the value of the final level of dissociation is 
adjusted until the test gas pressure Pr matches the acceleration gas pressure P 4 at 
the gas interface, giving the end flow conditions. 

These two models allow the prediction of the flow conditions of the test gas 
at the end of the acceleration tube given the initial fill pressures and the observed 
shock velocities The models assume that an equilibrium flow process is occuring. 

A comparison of the measured and predicted flow conditions for a low 
enthalpy dissociation shot using Argon driver gas are as follow 


AIRTESTGAS 

MEASURED 

PREDICTED 

ACCELERATIOirUBE 



Shock Velocity (m/s) 

4000 


TEST FLOW 



Dissociation Fraction 



Static Pressure (Pa) 

5158 

5170 

Pitot Pressure (kPa) 

230 

194 

Density (kg/ nr*) 


0.0197 

Temperature (K) 


912 


While the agreement has been good for low enthalpy test conditions , as 
yet the program has not been able to match static pressures for the high enthalpy 
shots The program is being modified to account for the Kkelyhood that total 
recombination does not occur for these conditions. 





To assist in the calibration of the simulation programs an adapter section 
with a stauc pressure probe was added to the end of the acceleration tube. This 
transducer allowed the measurement of the static pressure for the duration of the 
flow process providing further evidence of the steady state of the test flow as 

shown by the flat traces. (FIGURE 3) 

It wasalsoproposedto use the laser interferometer originally set up by Dr 
T. Posillico at the beginning of 1989, to measure the density changes in the flow 

and thus determine the free stream density of the test gas. 

The laser interferometer (FIGURE 5) passes two parallel beams from the 
same source through the test section. One beam passes through the test flow 
exiting the end of the acceleration tube while the other does not. The beams are 
optically recombined and made to interfere on a light detecting diode. The change 
in intensity of the combined beam can then be related to the change in density of 

the test flow. 


The system works on an infinite fringe pattern where 


Fringe shift 



P2-1V 
PS ; 


(Liepman & Roshko) 


where L = 

X = 

PS = 
Pi = 
P2 = 


length of flight across the test section which in this case 
is the width of the jet of gas exiting the tube, approx 32 mm 
wavelength of the light source i.e. the He - Ne laser 6318 nm 
density of gas at standard temperature and pressure, 
initial density of acceleration gas 
density of gas behind the shock 
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The laser interferometer was shown to be able to detect the arrival of a 
shock and qualitatively show the density change in the flow (FIGURE 4). The 
problem came when it was attempted to quantify the density change observed. 
The diode tight detector outputs a voltage proportional to the intensity of the 
interfering beams. For an infinite fringe pattern the intensity varies sinusoidally 
with fringe shift. 

Le. Intensity = a + b sin ( N + 0 ) 

Hence to calibrate the interferometer the starting point must be known and 
the maximum full scale change in intensity must be determined. The only way to 
determine these constants is to calibrate the system with a known density change. 
The electronics will only respond to a very high speed change of the order of time 
scale associated with flow in the tunnel. 

The tunnel was fired in shock tube mode in a series of tests to provide 
known high speed density changes , but it still proved inadequate to determine 
the required constants. It was decided that calibration of the system is impractical 
in terms of time and effort at this time but if an accurate and simple method of 
calibration could be found then the laser interferometer would be a useful n on- 
intrusive diagnostic tool. 
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FIGURE-1 


SHOCKTUBESECTION 


User inputs observed shock velocity in the shock tube 

and initial fin pressure in the shock tube 


▼ 

User inputs initial estimate of ot2 (or X2> 


Solve the quadratic for density ratio E 

„ i 

Calculate p- and thus P2 using e 

i 

»r 

Calculate £ and thus T2 using e 
M 1 


Substitute values into law of mass action 
and find error 


Adjust value of (*2 to reduce error until it is acceptable 


Using c*2 calculate other conditions at the end of shock tube. 
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FIGURE 2 


ACCELERATIONTUBESECTION 


User inputs observed shock velocity in acceleration tube 
and initial acceleration tube fill pressure 


1 


User inputs initial estimate of a final (Xf j na i ) 


-► Calculate U4 and P4 using ideal 1 - D shock relations 
set Uj = U4 


I 


Calculate the gas velocity, sonic speed and gas temperature at 
transition point 


1 


Calculate pressures using law of mass action 


Calculate Aj, Tj, Pr using ideal unsteady expansion relations 


Compare Pr and P4 


Adjust a final and repeat until Pr and P4 match 




Have end flow conditions 


175 



CH2 


CHANNEL Z* 

200 MILLIVOLTS/DIV 
TIMEBASE=25 . 6MICROSEC 
CAPTION ONLY*l 


/ DIV, 



FIGURE 3 : Static Pressure Trace of Flow At Exit of Acceleration Tube 


CHANNEL 3 
50 MILLIVOLTS/DIV 
TIMEBASE=51 . 2MIC-R0SEC / DIV, 
CAPTION ONLY® 1 


CH3 



FIGURE 4 : Laser Interferometer Trace 
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Diagnostic Beam 
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DISTURBANCES IN THE DRIVER GAS 
OF A SHOCK TUBE 


A. Pauli and R.J. Stalker 

Mechanical Engineering Department 
University of Queensland 


1. Introduction 

n has bean observed (Stalker, that expansion tube test tines are i» 
at low enthalpies by large regular disturbances which disrupt the centre ne 
pitot pressure. Previously, these disturbances have been attribute 
bubbles which are mixtures of driver and test gases and which are '*** ““ 
than the test gas. However, laser lnt.rforo.etry experiments (Possl 
which measured test gas density at the expansion tube exit detected a much 
smaller density change than was expected for a large bubble. This Indies « 
that either, the bubbles are elongated In the flow direction or. mo 
fundamentally, the bubbles do not cause these disturbances. In this report 
lt am be assumed that the latter is correct and an alternative theory which 
is based on longitudinal and lateral waves is developed which is expected o 
explain the structure of the pitot pressure disturbances. . 

It has been observed that similar disturbances exist in the tes an nv 
gases following a shock in a shock tube. It is believed that a model or 
these disturbances could be used to determine test time limitations in the 
more complex flows of an expansion tube. There are advantages of analysing 
shock tube flows rather than the flow fro. an expansion tube. The primary 
advantage is that shock tube flows are not accelerated by the unst 
expansion centred at the secondary diaphragm of an expansion tube . Th 
reduces the complexity of the flow history and thus reduces the possibility 
of additional unknowns distorting the fundamental problem. 

This report describes the disturbances which are observed in the driver gas 
of a shock tube. 


2 Experimental Description 

This report centres on a series of tests in which the exit wall pressure and 
centreline pitot pressure where measured in a shock tube (see figure . 
shock tube was 5.435 « long and has an inside diameter of 38 ... free 
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piston driver was used to heat the driver gas. 

Throughout the experiments the driver conditions remained the same. The 
filling pressure of the reservoir behind the piston was 3.4 MPa and the 
compression tube was filled with helium to 130 kPa. The helium driver gas 
was compressed to rupture a 1mm cold rolled mild steel diaphragm at nominally 

34.5 Mpa. 

Shock speeds were varied between 2.4 km/s and 6 km/s by altering the shock 
tube filling pressure between 320 mm^and 5 mm ^respectively. The shock tube 
and dump tank were filled with undried air extruded from the atmosphere. 
Shock speeds were measured with ionization gauges stationed at five locations 
as shown in figure 1. 

Wall pressure was measured 15 mm from the shock tube exit and the pitot probe 
leading edge was positioned centrally approximately 5mm from the shock tube 
exit. The pitot design is shown in figure 2. For all experiments each analog 
signal was sampled every microsecond. 


3. Results and Observations 

Figure 3 displays the pitot pressure as a function of time measured for 
different shock speeds. Figure 4 displays the corresponding wall pressure 
records and table 1 displays the shock speeds calculated from the time delay 
between excitation of neighbouring ionization gauges. A typical record from 
the ionization gauges is shown in figure 7. The predicted driver-test gas 
interface location is displayed in figures 3 and 4. This prediction was made 
using turbulent boundary layer entrainment theory (Mlrel). The shock speeds 
where assumed to be that observed between the two ionization gauges closest 
to the diaphragm. 

Verification of the repetitiveness of these results is given in figures 5 and 
6 where the the shock tube filling pressures of 5, 10 and 20 mm were 
repeated. In these experiments a more sensitive wall pressure gauge was used. 
In can be seen that the results from these series of experiments are in 
reasonable agreement with each other. 

Three features of the unsteady nature of each pair of wall and pitot traces 
should be observed. 

(i) A sharp dip is sometimes observed in the pitot pressure in the vicinity 
of the interface. 

(ii) At the higher shock speeds the pitot pressure mean does not follow the 
wall pressure mean. 

(iii) The pitot pressure has well defined oscillations superimposed on its 
mean, but the oscillations in the wall pressure are not as pronounced. 

179 


I 



third feature Is believed to be a fluid mechanical feature and Is not the 
"u of stress saves In the pitot probe. This -as verified si. an 
experiment in shieh the tapping to the pitot pressure transducer was - 
Virtually no vibrational noise Is detected from the pitot pressure 

transducer. 


3 ’ 1 ’ ^dieted 1 * interface lies close to the dip or a point where the nature of 
p'tot pressure traces changes. It sill be seen from the theoretical 
resuns discussed later that in addition to a mean pitot pressure change^, 
to the difference betseen the the 

r^r :: z. used a. most 

interfa . „ , Ideally to the left and right of the 

positions are shown in figures 3-6. Ideally to 

interface the gases are air and helium respectively. 


prom m be seen for a shock tube filling pressure (P.) of S mm 

"" after the observed driver-test gas Interface the sail pressure rise 
Idily. The gradient of this rise becomes less as P is ™ ” 

r . gomm shere there is no rise. If P. U greater than SO mm 

hinrl aslt of this trend is seen In the pitot pressure record. For 

P 5., the pitot pressure rises rapidly after the Interface and then rises 

sloulT (If at all, for the remainder of the trace. This trend is repeated 

1 ... the rise m pitot pressure being less dramatic as 

for P = 10mm and 20mm with the ris p 

the shock tube filling pressured is increased. An interesting result oc 
It P.= 40mm. At this filling pressure, after the interface, the pitot 

Dressure falls as the wall pressure rises. 

It can 1 seen that the long term overall trend is for the pitot pressure to 
have a lesser gradient than the wall pressure. 


3 3 Oscillation features 

Osc nations can not be seen in the -all pressure traces of figure 4 a h 

h gher shock speeds. It is only -hen the disturbances are -treme a t 
' shock speeds that oscillations appear in the -all pressure. This is 
t a y ue to the insensitivity of the gauge used. In figure 6 -here the 
sensitivity of the -all pressure gauge -as approximately t-enty five greater 
.h„„ that used for figure 4 it can be seen that high frequency oscillations 
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do exist at the higher shock speeds, however, these oscillations are still 
small compared with those observed in the pitot pressure. 

The oscillation frequency increases with shock speed and at the three highest 
shock speeds the high frequency oscillations are superimposed on a much lower 
frequency disturbance which has a period of approximately 250 ps. 

4. Theoretical Results 

In this section a model based on longitudinal and lateral acoustic waves is 
developed which is consistent with the oscillatory nature of the flow. In 
addition, it is shown that identification of the compression following the 
dip with an unsteady expansion travelling upstream can explain the different 
trends observed between the pitot and wall pressure records. 

4.1. Acoustic waves 

Acoustic disturbances are those in which the sound speed remains constant to 
first order. If the velocity and pressure are perturbed about a mean value 
then it can be shown (Whitham) that to first order there exists a potential 0 
for which 


u*= 70 (41) 

P * p o- Po (42) 

P = Po- Po c ~ 2 It <43) 

c 2 = P ! (P 0 ) = r <«- 4 > 

Po 

where y,u*,P,p c and t are the specific heat ratio, velocity, pressure, 
density, sound speed and time respectively and subscripts "o" refer to the 
condition which is being perturbed. The spatial origin is stationary with 
respect to the mean flow so that the unperturbed velocity is zero. From 
conservation of mass and momentum it can be shown that 0 obeys the wave 
equation 


i+ = c 2 v 2 0 

at 2 


If a change of origin is made and is fixed with respect to the shock tube 


l 
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then (4.1) Is rewritten as 


U = Uq* 


( 4 . 6 ) 


c ,. .. 4) and (4.6) and the 1 ' 1 « h PUot P" ss ' ire f ° r “ la “ 

hi 0 - u first order in the derivatives of ♦ that for supersonic flows the 

pitot pressure Is 


9 = Tr 


j _ L — ♦ 
1 c 2 St 


2y(2Mn-l) \a± 


luo(2yMo-r + l)J 


( 4 . 7 ) 


. * d <(> a 

k the oosltlve z-dlrection is In the unperturbed flow direction u - 3l • 
where the * „ lth respect to ti.e Is the contribution 

I 7pr" fluctuations and the tern contain^ the derivative with 

respect to z is the contribution to to velocity fluctuations. 


«, 2 . Wav. elation solutions ^ uhlch have th e properties 

In this section solutions btalned It will be assumed that 

observed in the experimental results are obtained It will 

- — - *» - 
i“ — - -» « - — « — 

is slip but no penetration so that 


^ (z.a.e.t) = 0. 
or 


( 4 . 8 ) 


solutions of the wave equation obtained by separation ^ 

variables and which are single valued, periodic in time and 

have the form 


l(ut±p n Z±me) 

<t> = W) e 


( 4 . 9 ) 


where m and n are integers. 


■!(?)•- 4 


( 4 . 10 ) 


a „d J is an . th. order Bessel function of the first kind, 
condition (4.8) is imposed on (4.9) then this requires that 


If the boundary 
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( 4 . 11 ) 


jl(X n a) = 0. 

Hence, A n a' is the n th. root of (4.11). 

Any mode or any combination of modes given by different values of m and n may 
exist, however the properties of most of these modes will not fit the 
properties observed in figures 3-6. For example, if m*0 then the m th. order 
Bessel is zero at r=0, however, experimental results show that the centreline 
pitot pressure is not less, but considerably larger than the wall pressure. 
The only solutions of the form (4.9) with this property are the zero th order 
(m=0) Bessel function solutions. The zero th order Bessel function has its 
largest maxima at r=0 (figure 8) and oscillates to zero as r approaches 
infinity. If it is therefore assumed that the required solutions have m-0 

then as -Jj is equal to the derivative of J 0 it follows that A n = Z n /a, where 

is the n th zero of which are available from tables (Abramowitz) . 

If n-0 then A is zero (as Zq is zero) and solutions are strictly dependent on 
z and time. Solutions of this form will be referred to as longitudinal waves. 
In contrast to the zero th order mode all other solutions have a radial 

dependence and thus will be referred to as lateral waves. A first order 

lateral wave has n=l and A n =3.83/a. Pressure contours of this mode in the 
(z.r) space are plotted in figure 9. Higher order lateral waves may also 
occur, however, at this stage studies will be restricted to the first order 
lateral waves and the longitudinal waves on the grounds that more energy will 
be needed to excite higher order modes and therefore the lower order modes 
would be excited first. 

4.3 Pitot pressures of longitudinal and lateral waves. 

From (4.7) it can be seen that if a disturbance is travelling in the flow 
direction (^(z-ct)) the velocity fluctuations enhance the pressure 
fluctuations to give a larger percentage fluctuation in pitot pressure than 
would be seen in the wall pressure. If the disturbance travels upstream the 
effect on pitot pressure due to pressure fluctuations is decreased by the 
velocity fluctuations. 

If the peak centreline pressure fluctuation is P r then from (4.2) and (4.9) 
along the centreline 

, _ _i _Pj_ (wt±0 n Z) ,4.12) 

9 Po u 

Hence, the magnitude of the centreline pitot pressure fluctuation is 
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( 4 . 13 ) 


9 , 


0 Pr 


_ 2{2M 2 - 1) R 

[ T M(2rM 2 -T+l) 


where the Mach number M = U 0 /c and 


B„= 


1 - (¥) ■ 

^ J 


( 4 . 14 ) 


If y=1.667 then the coefficient 


2(2M 2 - 1 ) 
M(2yM 2 -7+l) 


( 4 . 15 ) 


, M - 1 8 In the case of a longitudinal wave B 0 is equal 

is equai ;:;:M. i i — — *- — - 

l ° r is less than three the contribution fro. longitudinal wave velocity 
number Is less th fluctu ,tions Is comparable to that from 

fluctuations to the pitot pressure travelling In th. 

the static pressure variations. If a W and tt> . pltot 

fl n„ direction .then the ^ terms br-e-J ^ ^ ^ pressure 

pressure fluctuation will b PP^ ^ ^ lhen the bracked 

” ni subtract, hence. It will be difficult to detect in the pitot 

. ....... .... .... 

:r:rjssr::s~. - «. •- 

if a iongi lateral wave travelling in the flow 

nressure fluctuation observed for a lateral wa 

d rectlon would be less than that of the longitudinal wave because B . 

P rthermore a lateral wave travelling upstream would produce a greater pitot 
pressure fluctuation than the longitudinal wav. travelling upstream. 

can also be seen for constant Hash number that as the sound s eed 
lateral wave will produce a greater fluctuation in the pitot 
decreases irop ortant in the application of acoustic wave theory 

pressure. This may be i p ^ ln the sound speed across 

to the flow in an expansion tube sin lateral 

an expansion will tend to increase the relative importance of the 

waves. 
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4.4. Doppler shifts 
If the change of origin 

X = ut + Z (4.16) 

is made so that the origin of the (x.r.t) space is fixed with respect to the 
shock tube then the coefficient of the exponent in (4.9) with m=0 is 

f 

i + u£ n jt ± x0 n 

Hence, a Doppler shift occurs in the frequency due to the velocity of the 
gas. If a wave Is travelling in the direction of the flow (the positive 
x-direction) then the frequency is increased to w + u0 n . However, if waves 
are travelling upstream then the frequency decreases to u - u0 n . It should 
be seen that if the velocity is sufficiently large so that w - u0 n < 0 then 
this wave would also represent a wave travelling in the positive x-direction. 
In the simple case of a longitudinal wave this would occur if u>c. 

It should also be observed from (4.17) that the frequency shift is dependent 
upon sound speed. For a longitudinal wave the frequency shift is w(l±u/c). 
Thus across an interface where the sound speed changes there would be a 
change in the period of the disturbance which is dependent upon the Mach 
number either side of the interface. A similar result is true for lateral 

waves, however, the frequency shift is dependent upon u,c,u and a 

independently. 

4.5. Reflection and transmition coefficients 

jf a wave travelling in the positive 2 — direction strikes an interface at 

which the sound speed and the specific heat ratio change then the reflection 

and transmition coefficients are obtained by requiring that the pressure and 
velocity be continuous across the interface. If an incoming longitudinal 
wave has unit amplitude then the transmitted and reflected waves have 
amplitudes of 

2*1 Pnl*-^ 

(4.18) 

*1 @n2^2 + *20nl^l 

and 

*2^111^1 ~ *lPn2^2 

(4.19) 

*lPn2^2 + *20nlCl 
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respectively (n = 0), where 


Pn. 


■ l »,)•• < f 


( 4 , 20 ) 


which the incident and 


and subscripts 1 and 2 refer to the mediums 
transmitted waves respectively traverse. 

:f the incident wave Is a lateral wave then It Is not pass lb . to 
collectively satisfy continuity in the pressure, the longitudinal veloc y 
component (*/»z> and the lateral velocity component 0«/ar) across the 
interface. One of these restrictions has to be Ignored. To be consisted 
with the wall boundary condition, slip at the interface will be tolerate . 
With this assumption the reflection and transition coefficients -ill be the 
same as that for a longitudinal wave. 

There is an additional point to consider for lateral wave transmitlon whlc 
does not arise with longitudinal waves. There is the possibility that the 
transmitted lateral wave can not be supported in the new medium. If reg on 

supports a lateral wave then o > C.V however If the sound speed in region 

.. . v r \ then Bn2 will be imaginary and thus 

is sufficiently large so that u < C 2 X n then pm win 

the transmitted wave will decay exponentially. 

4.6. Receding unsteady expansions 

this section the pitot pressure properties of a receding unsteady 
expansion are determined as a function of Mach number. Fro. these properties 
it will be concluded (section S) that the compression observed following the 
driver-test gas Interface can be modeled as an unsteady expansion travelling 
upstream. The unsteady expansion Is assumed to be a disturbance for 
is too large to ignore sound speed variations. Thus, the propert es o 
one-dimensional isentroplc Inviscld flow are used. 

Across an unsteady adiabatic expansion travelling upstream into an ideal 
inviscld isentroplc gas the velocity and sound speed are related to the 
steady upstream condition (subscript o) by 


u + 


FI c - J=I °° 


( 4 . 21 ) 


The wall pressure and sound speed 




_ j. 1 


The Rayleigh pitot pressure formula written in terms of wall pressure and 
Mach number (M=u/c) is 


¥ = 



(4.23) 


From (4.21) the sound speed across an unsteady expansion may be written in 
terms of Mach number and thus from (4.22) and (4.23) the wall and thus the 
pitot pressure can also be written as functions of the Mach number. 

By examining the derivative of the pitot pressure with respect to the Mach 
number it can be shown that 

— < 0 ; M > 1 * ^ a 1.62 (4.24) 

dM 2 


and 


d¥ 

dM 


> 0 


1 < M < 


1 + 5 ‘ 


(4. 25) 


independent of specific heat ratio. 

If an unsteady expansion is travelling upstream in a flow with Mach number 
greater than one then the wall pressure will rise and the Mach number will 
fall as the expansion is convected downstream past a fixed point on the tube 
and thus the expansion will appear as a compression. If the flow Mach number 
at the beginning of the compression is greater than 1.62 then since the Mach 
number is decreasing with time the pitot pressure will increase until the 
Mach number is 1.62. If the Mach number decreases further the wall pressure 
shall continue to increase but the pitot pressure wall fall. 

5. Receding Unsteady Expansion Results 

The following three observations show that the compression following the 
interface observed in the wall pressure measurements can be modeled as a 
receding unsteady expansion. 

1. The expansion has a decreasing gradient as the flow Mach number upstream 
of the expansion decreases. This is consistent with an expansion travelling 
upstream. The expansion head travels upstream at a speed c 0 with respect to 
the flow and is convected downstream with a velocity Uq ' and as the shock 
speed decreases Uq decreases and c 0 increases thus the separation between the 
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expansion head and the Interface will Increase. If “o' c 0 then the expansio 
dead is stationary with respect to the tube, so the observed compression 

would be very slight. 

2 Table 2 Hats the pitot pressure (>-.) and wall pressure IP.) which are 
measured at the expansion head at the four highest shoe* speeds Fro. these 
two quantities the flow Mach number (M.) at the expahsion head Is 
using Rayleigh pitot pressure formula (4.23). The fifth column of Table 2 
lists the wall pressure at the driver-test gas interface, which is assume 
be the trailing edge of the expansion. It Is then assumed that the expans on 
is unsteady and the Mach number behind the Interface (M.) Is determine u 

the relation 



(5.1) 


which is obtained fro. the unsteady and adiabatic relationships (4.25) and 

(4.22) respectively. . ... 

. . i j ska interface is determined with 

Knowing M, and P x the pitot pressure behind the interlace 

i *7 Tt ran be seen that these values agree 

(4.23) and is tabulated in column 7. It can be seen 

very well with those observed (column 4). 

At lower shock speeds than those tabulated the agreement is not as good. 

i- a tho wall trace for P« = 80mm has only a 
However, the compression observed , 

slight gradient which would indicate that the upstream Mach number 
a preaching one (and the Mach number is less for P.oSOmm). thus the ito 

probe may be outside the Mach cohe In which case the above analysts is 

longer valid. Theoretical Ideal gas calculations indicate that the driver 
gas is at Mach 1 when the shock speed Is approximately 3.3 km/s which s 
consistent with this deduction as this theoretical shock speed ts on y 
slightly less that that measured (see Table 1). 

3 It was seen In the previous section that an unsteady expansion vou 
produce a rise In both the pitot and wall pressures provided the Mach num er 
1 . greater than 1.62 and the pitot pressure would fall as the wall pressure 
rose if the Mach number was less than 1.62 (and greater than one). It can be 
seen from figures 3 and 4 that this trend is observed, however, the turning 
point would appear to be approximately Mach 2 rather than Mach 1.62. 

6 Spectral Distributions of Measured Disturbances 

Fourier analysis have been used to study the frequency distributions of the 
pitot pressure disturbances in the driver gas. The frequency spectrums 
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figure 10 where made using fast Fourier transform techniques over 256 points 
from data sampled at a rate of 1 MHz. 

It can be seen that there is a definite band of frequencies different from 
zero at which the disturbance is dominate. The band width is a function of 
filling pressure. It can be seen that the Fourier decomposition of the 
repeat pitot pressures (figure 5) have similar forms. The decomposition has 
been taken over the time period in which the unsteady expansion was detected. 
It can be seen from 4. 10 that for a given sound speed that the fundamental 
frequency of a lateral wave must be greater than cA n . All lateral waves with 
fundamental frequencies less than this will decay. It should be observed 
that there Is not a non-zero lower bound on the fundamental frequency of a 
longitudinal wave since A 0 =0. 

From Table 2 the Mach number immediately behind the interface (M 1 ) for P s =5mm 
is 2.57. The shock speed measured (Table 1) was 5.82 km/s. Hence, for an 
ideal gas (y=1.4) the gas speed would be 4.85 km/s. Therefore, the sound 
speed of the driver gas behind the interface would be 1.89 km/s. For a first 
order lateral wave (n=l) A 1 =201.6, hence, the fundamental frequency of a 
first order lateral wave must be greater than 380450/2w=60550 Hz. From 4. 17 
because B n >0 it then follows that in the laboratory frame of reference the 
frequency of a first order lateral wave travelling in the direction of the 
flow would also be greater than this frequency. Following the same arguments 
it can be shown that the frequency of a second order lateral wave travelling 
in the flow direction must be greater than 110980 Hz. It can be seen from 
the Fourier decompositions of figure 10 that a large number of the dominant 
frequency components are between the lowest observable frequencies of the 
first and second order lateral waves for P B =5mm (Similar statements can be 
made for P B =10mm and P B =20mm). From this alone it can not be concluded that 
these disturbances are lateral waves because longitudinal waves can exist at 
any frequency. 

However, the dominance of the lateral wave would be consistent with the 
reduced amplitudes observed in the wall pressure disturbances. The Bessel 
function J 0 (A 1 a) s 0. 4 J o (0), thus the static pressure along the wall is 0.4 
that of the centreline static pressure. This enhances the increase in pitot 
pressure fluctuations which results from the contributing velocity component 
from a disturbance travelling in the flow direction (section 4). 

Furthermore the dominance of lateral waves is also consistent with 
Possilico’s interf orometry results which did not detect a large variation in 
the integrated density across the test gas. In a lateral wave the density 
variation is both positive and negative across the test section and thus the 



Integrated density would be reduced. It can be shown that 

J 0 (\r)dr 


L 


f a 


1-08 1 
3.83 4 


( 6 . 1 ) 


so that the fringe shift which would be observed In a lateral wave would be 
approximately one duarter of that for a ^ dl5turbanC es 

r/n cld lt be explained through the existence of longitudinal waves can 
w. understood if the disturbances are lateral waves. 


7. Hypulse Results which It was a 

The ,,,L Hypulse facility hellu ,. the filling 

shock tube. , 8 5 nun respectively, 

pressures of the driver and shocX ^ „ shocX 

The inside diameter of the Hypulse _ flgura „ t hat 

ZtrZ'TZZ to those observed In the University of Queensland T Q 
facility are also observed in the Hypulse facility. 

A furth er o rr . n r — :irr — 

— ’.r=,r . - — •• - 

centre which is consistent with lateral waves. 


8. conclusions uhlch slarts mediately 

1 An compression Is observed in the drl g . , an 

behind the driver-test gas interface. This compression can be modelled 

unsteady expansion travelling upstream. * „ the dominate 

2. First order lateral waves (n-1) wouxa pp 
nc=ri llitory disturbance which exist in the driver g 
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Table 1. 

SHOCK TUBE FILLING 


PRESSURE 

(mm) 

320 

160 

80 

40 

20 

10 

5 

SHOCK SPEED BETWEEN 
GAUGES (km/s) 








1 AND 

2 

_ 

_ 

4.65 

5.06 

5. 44 

5.72 

6. 19 

2 AND 

3 

- 

3.48 

4. 15 

4.69 

5. 17 

5.76 

6. 13 

3 AND 

4 

2. 39 

3. 10 

3. 60 

4. 19 

4.82 

5. 50 

5. 97 

4 AND 

5 

2. 12 

2.69 

3. 18 

3.79 

4.41 

5. 27 

5. 82 

Table 2. 

p, 

mm 

*0 

kPa 

P 0 

kPa 

Mo 

^l»ea« 

kPa 


Pi 

kPa 

Mi 

^lcalc 1 

kPa 

5 

3350 

500 

2. 05 

2950 


300 

2.57 

3070 

10 

3150 

550 

1.89 

2880 


350 

2. 35 

2980 

20 

3350 

670 

1.75 

3690 


460 

2. 13 

3240 

40 

3880 

830 

1.69 

4550 


560 

2.07 

3800 
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Figure 1. Shock tube and test section. 
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HI Host likely interface position 
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Figure 3. Figure 4. 

Pitot pressure as a function of time for Wall pressure as a function of time for 


DIFFERENT SHOCK TUBE FILLING PRESSURES. DIFFERENT SHOCK TUBE FILLING PRESSURES. 
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Figures 5 6. Com>ar,son of pitot a® wall presses between 
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Figure 6 1 
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Figure 11. Hypulse pitot pressure as a function of time at different 

RADIAL LOCATIONS. 
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